^D-A14i  210 

UNCLflSSIFIE 

KINETIC-ELASTIC  APPROACH  FOR  TIME-DEPENDENT  RHEOLOGICAt 
DATA  ON  SLURRV  FU.  .  (U)  OHIO  STATE  UNIV  RESEARCH 
FOUNDATION  COLUMBUS  R  S  BRODKEV  ET  AL.  JAN  84 
OSURF-762084/712747  AFWAL-TR-82-2085  F/G  20/4 

1/3 

NL 

1 

■ 

■ 

■ 

■ 

1 

J 

■ 

ft 

K 

IB 

^B 

^B 

^B 

^B 

^B 

^B 

^B 

Bi 

AD-A141  210 


A  FirJAL-  TR-G3-  2UHS 


NOTICE 


when  Government  drawings,  specifications ,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications ,  or  other  data,  is  not  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS ,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


_ 1  ^  ^ l  _ 

TERESA  A.  BOOS 

Fuels  Branch,  Fuels  and  Lubrication  Division 
Aero  Propulsion  Laboratory 


ARTHUR  V.  CHURCHILL 

Chief,  Fuels  Branch 

Fuels  and  Lubrication  Division 

Aero  Propulsion  Laboratory 


FOR  THE  COMMANDER 


R03ERT  D.  SHERRILL,  Chief 
Fuels  and  Lubrication  Division 
Aero  Prooulsion  Laboratory 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or 
if  the  addressee  is  no  longer  employed  by  your  organization  please  notify  AFWAL/POSF 
w-PAFB,  OH  45433  to  help  us  maintain  a  current  mailing  list". 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations ,  contractual  obligations,  or  notice  on  a  specific  document. 


Unclassif ied 


SECuR'Ty  CLASSIFICATION  of  this  PAGE  (TWioi  Dmto  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

t.  REPORT  NUMBER  |2.  GOVT  ACCESSION  NO. 

A!  1A;  -  '.R- 8  3-2085  if>j\  .  Al  v  /  7^1  O 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

KINETIC-ELASTIC  APPROACH  FOR  TIME-DEPENDENT 
RHEOLOGICAL  DATA  ON  SLURRY  FUELS  AND  POLYMERS 

S.  TYPE  bf  REPORT  A  PERIOD  COVERED 

Final 'Report  for  Period 

Mar  80  -  Sep  82 

6.  PERFORMING  ORG.  REPORT  NUMBER 

762084/712747 

7.  AUTHOR/*) 

Robert  S.  Brodkey 

(K.  Park,  S.  Weng,  A.  Ouibrahim,  S.  F.  Lin) 

8.  CONTRACT  OR  GRANT  NUMBER/*; 

F33  61 5-80-C-2021 

9.  performing  organization  name  and  aooress 

The  Ohio  State  University 

Research  Foundation,  1314  Kinnear  Road 

Columbus,  Ohio  43212 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  8  WORK  UNIT  NUMBERS 

P.EU  62203F 

3048-05-91 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Aero  Propulsion  Laboratory  (ATVAL/POSF) 

12.  REPORT  DATE 

January  1964 

Air  Force  Wright  Aeronautical  Laboratories  (/FSC) 
Wright -Patterson  Air  Force  Base,  Ohio  45433 

13.  NUMBER  OF  PAGES 

212 

14.  MONITORING  AGENCY  name  4  ADDRESS (It  dlllorent  /ran  Controlling  Olltcm) 

15.  SECURITY  CLASS,  (ol  Nil*  report) 

Unclassified 

is*.  declassification/ downgrading 

SCHEDULE 

16.  01 S  ’  RIBUTION  STATEMENT  f  •  thlo  Report) 

Approved  for  public  r  lease;  distribution  uni  mited 

• 

17  DISTRIBUTION  STATEMEN  >t  the  obotroet  ontorvi  in  Block  10,  II  d  eronl  Irom  Report) 

ta.  st.  'Plementapy  not*;s 

1 

19  KEY  WORDS  rcontlnuo  on  r evoroo  *  do  It  nocootory  ond  Idontlly  by  bio  number; 

rheology,  slurry  fuels,  JP-10,  jet  fuel,  carb  n  black  slurries,  sheer  thinning, 
sheer  thickening,  thixotropic,  dilation,  polyn er  solution,  lubricants, 

r __  . 

20.  \a^oTRAC  T  (Contlnuo  on  reverse  tldo  II  neceeeory  ond  Idontlly  by  bloc!  nutnbor) 

'’’At  The  Ohio  State  University,  we  have  been  developing  our  kinetic  interpre¬ 
tation  of  non-Newtonian  fluid  behavior  under  support  from  tbe  Aero  Propulsion 
Laboratory  with  materials  of  interest  to  the  Air  Force.  Specifically,  we  have 
considered  slurry  fuel  systems  and  a  reference,  high-viscosity  lubricant  (5P4E) . 

In  addition,  we  continued  to  obtain  data  on  a  polymeric  system  previously 
studied  by  us  so  as  to  ascertain  the  adequacy  of  our  measurements  and  of  our 
theory, 

(continued  on  reverse  side) 

DO  ,  r*,  1473  EOITi  >N  '  -  1  NOV  «S  IS  OSS  7LETE 


_ _ L'nclat,  if  i  ;d _ 

5ECUR  TV  CL  ASSIFICATIO'  OF  T  US  PAGE  (When  Dot'  -ed) 


f 


SlftiniTY  c  Av.iFir.AriON  OF  This  P AGEfHTiwi  Data  Entered) 


Block  20;  Abstract  (continued) 

\  i  -  - — — 

^The  kinetic  theory  is  a  phenomenological  rate  concept  used  to  describe  | 
the  time  rate  of  change  of  a  material  structure.  We  call  this  part  of  the 
development  "the  kinetic  model  for  thixotropic  change."  Corresponding  to  the 
thixotropic  fluid  structure  is  a  thixotropic  structural  viscosity  and  a  cor¬ 
responding  thixotropic  stress.  The  real  stress  in  the  system  is  different 
from  the  thixotropic  stress  because  of  the  elastic  properties  of  the  material. 
Our  "elastic  model  for  viscoelasticity"  utilizes  a  modified  version  of 
Oldroyd's  development  in  convected  coordinates.  A  combination  of  the  two 
models  allows  the  representation  of  rheological  data  with  what  we  call  "the 
kinetic-elastic  model." 

To  date  we  have  shown  conclusively  that  much  of  the  older  data  could  not 
be  used  because  of  instrumentation  problems.  We  have  obtained  new  polymer 
data  (polymethylmethacrylate  in  DEP,  PMMA/DEP)  wi  th  our  improved  instrumenta¬ 
tion  and  data  acquisition  system.  Using  this  new  data  has  allowed  for  the 
first  time  a  critical  test  of  tK&  kinetic-elastic  theory.  The  comparison 
between  the  polymer  data  and  the  theory  are  very  good.  We  wish  to  study 
further  the  variation  of  parameters  and  a  possible  modification  of  the  elastic 
part  of  the  theory.  We  have  completed  our  measurements  of  a  number  of  low 
viscosity  slurry  fuels  and  have  completed  an  extensive  data  gathering  program 
on  the  reference  lubricant.  The  slurry  fuels  are  clearly  non-Newtonian,  and 
thixotropic,  and  some  have  a  low  critical  yield  stress  that  must  be  exceeded 
before  they  will  flow.  These  materials  are  Theologically  complex  and  turn 
out  to  be  similar  to  another  colloidial  material  we  studied  some  years  ago. 

We  have  also  developed  the  kinetic-elastic  theory  to  account  for  the  yield 
stress.  The  comparison  of  slurry  data  to  the  theory  is  adequate  but  the  logic 
use  could  be  improved.  The  reference  lubricant  was  found  to  be  Newtonian  at 
room  temperature.  As  the  temperature  was  reduced  towards  the  pour  point 
(-5  C),  the  material's  viscosity  Increased  rapidly  and  the  material  became 
non-Newtonian  and  viscoelastic.  We  have  formulated  and  applied  to  a  test  set 
of  this  data  a  method  to  remove  the  effect  of  the  upper  cone  movement  during 
transient  testing.  Since  this  material  is  so  much  like  the  polymeric  material 
we  plan  to  eventually  reduce  the  data  and  Interpret  the  results  in  terms  of 
the  kinetic-elastic  theory  as  used  for  the  polymeric  data.  I 
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The  pre  ent  work  has  been  carried  out  under  contract  number  F33615- 
80-C-2021  w  th  Wright-I’atterson  Air  Force  Base  over  the  period  from 
March  24,  1('80  to  September  30,  1982. 

Many  individuals  have  contributed  to  this  work.  In  particular, 
graduate  students  Drs.  Kyon  Park  and  Shew-Fen  Lin  spent  many  hours 
working  on  rheology  in  the  course  of  obtaining  their  degrees  in  Chemical 
Engineering.  Drs.  Sue  Weng  and  A.  Ouibrahim  likewise  spent  many  hours 
on  post-doctoral  appointments,  the  former  working  both  on  the  polymer 
solutions  and  slurry  aspects  and  the  latter  on  lubricants. 

In  addition,  Mr.  Chung  Pong  Chung  is  currently  investigating  pol /mer 
solutions  in  even  more  detail.  The  time  and  suggestions  offered  by  Dr. 
H.  Lander  aid  T.  Boos  of  the  Fuels  Branch  at  Wright-Patterson  are 
i  greatly  appreciated  as  well  as  the  loan  of  key  equipment  for  the  low 

temperature  work.  Finally,  The  Ohio  State  University  was  also  a  major 
contributor  through  additional  student  support  and  by  providing  the 
funds  neede>  for  the  data  acquisition  system. 

Of  key  nporcance  to  the  slurry  st  idles  were  the  samples  provide  I  to 
us  by  the  F>  ’Is  Branch  that  were  obtai  led  from  Sun  Techand  from  Exxon. 

There  a  J  available  through  t’nivei  ;ity  Microfilms,  Ph.D.  theses  ’ 
Park  (1981)  and  'in  (1983)  that  provie  •  more  details  on  some  aspects  of 
this  study. 
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OBJECTIVES 


Our  charge  in  the  study  covered  by  this  final  report  was  threefold  : 

1)  To  establish  the  adequacy  of  the  time  dependent  polymer  data  previously 
obtained;  if  inadequate,  to  obtain  a  new  set  of  data  and  to  evaluate  the 
kinetic-elastic  approach  for  the  representation  of  this  data. 

2)  To  obtain  rheological  data  on  a  variety  of  potential  carbon  black  - 
JP-10  slurry  fuel  formulations  in  order  to  help  establish  their  general 
rheological  characteristics  as  a  function  of  shear  level  and  temperature 
and  to  adapt  the  kinetic-elastic  approach  for  the  description  of  such 
materials. 

3)  To  obtain,  if  time  permitted,  a  data  set  for  the  high  viscosity 
reference  lubricant  5P4E. 

Our  objectives  have  been  accomplished  and  this  report  provides  the  detailed 

summary  of  our  efforts. 


SIGNIFICANCE 


Systems  that  involve  rheological  fluids,  such  as  those  found  in  polymer 
processing  and  in  the  fields  of  fuels  and  lubricants,  often  involve 
rapid-response,  unsteady-state  conditions.  A  rheological  model  must  be  able  to 
adequately  represent  data  in  the  time  domain.  The  model  should  also  be  able  to 
account  for  changes  in  concentration  and  temperature.  As  important  is  that  the 
constants  and  parameters  of  the  model  have  basic  meaning  that  can  be  interpreted 
in  terms  of  the  fluid's  composition  and  can  be  related  to  the  performance  of  the 
fluid  in  its  application.  Ideally,  one  would  like  to  know  enough  about  the 
fluid's  rheology  so  that  modifications  of  the  fluid  can  be  made  in  response  to  a 
change  in  the  fluid's  performance  in  a  system.  Such  change  would  be  caused,  for 
example,  by  more  severe  conditions  or  by  a  change  in  the  physical  system  itself. 
The  "kinetic-elastic  model"  developed  at  Ohio  State  has  the  potential  of 
representing  complex  rheological  data.  The  general  approach  of  our  effort  is 
one  of  a  complete  characterization  of  any  given  non-Newtonian  material  so  that 
one  can  apply  this  information  to  evaluate  the  flow  of  the  material. 

The  work  has  significance  with  respect  to  the  flow  of  various  exotic- 
materials  such  as  thick  slurry  fuels,  shale  oil  products,  coal-oil  mixtures,  and 
modern  lubricants.  'I here  is,  for  example,  amazingly  little  known  about  the  flow 
properties  of  slurry  materials  with  yields.  There  is  no  adequate  means  of 
representing  such  data  when  the  material  has  viscoelastic  characteristics.  It 
is  a  simple  matter  to  extend  our  approach  to  allow  description  of  materials  of 
this  nature.  Such  an  analysis  would  apply  to  slurry  fuels  and  coal-oil 
mixtures.  ft  has  only  been  in  recent  years  that  workers  involved  in  lubrication 
have  come  to  realize  that  materials  they  are  dealing  with  are  highly 
non-Newtonian  and  viscoelastic  in  their  chai  icteristics.  Indeed,  modern 
lubricants  contain  a  high  amount  of  polymeric  additives  and  these  are 
notoriously  viscoelastic.  It  is  interesting  to  speculate  that  the  elastic 
characteristics  involved  in  sudden  deformation  could  be  an  explanation  as  to  why 
some  lubricants  fail  at  high  stress  conditions  at  low  temperatures. 


SUMMARY  OF  RESEARCH  ENDEAVORS 


In  our  original  proposal  we  provided  a  detailed  summary  of  the  field,  our 
earlier  endeavors,  and  a  detailed  development  of  the  theory.  We  repeat 
this  material  here  as  Appendix  D.  Included  in  the  review  is  a  literature 
survey  of  the  limited  work  done  by  others  on  lubricants  and  on  slurry  fuels. 

The  various  approaches  towards  the  elucidation  of  the  rheological 
characteristics  of  materials  generally  fall  into  one  of  five  classifications: 
empirical,  phenomenological  (mostly  rate  processes),  linear  viscoelastic, 
nonlinear  viscoelastic,  and  microrheological  analyses.  The  empirical  methods 
correlate  data  by  curve-fitting  techniques.  The  rate  theories  have  as  their 
basis  the  assumption  that  the  nonlinear  characteristics  can  be  associated  with 
some  structural  change  of  the  material  whether  it  involves  particle  associations 
and  dissociations,  link  formations  and  ruptures,  or  molecular  entanglements  and 
disentanglements.  The  linear  viscoelastic  models  are  based  upon  linear 
combinations  of  hook's  law  of  elasticity  and  Mewton's  law  of  viscosity.  The 
nonlinear  viscoelastic  models  are  based  on  cc  itinuum  mechanics  and  nonlinear 
combination  mechanical  models.  Finally,  mice  orheological  analysis  starts  with 
the  basic  moLecular  or  microscopic  variables  such  as  particle  sizes,  molecuiar 
interactions,  and  chain  lengths.  \  simplified  mechanism  is  proposed,  and  then 
it  is  mathematically  represented  a. id  solved. 

Each  approach  has  made  contributions  to  the  field,  but  with  the  presenL 
state-of-the-art,  we  cannot  completely  describe  the  non-Newtonian  behavior  even 
in  simple  geometries.  Steady-state  flow  behavior  has  been  extensively 
investigated,  and  the  representations  available  are  adequate.  However,  little 
is  understood  about  the  unsteady-state  flow,  such  as  shear  stress  and  normal 
stress  growth  after  the  onset  of  a  sudden  shear  rate.  The  theories  in  the 
literature  are  not  always  able  to  predict,  for  example,  stress  growth  at  a  given 
constant  shear  rate  from  data  obtained  on  other  experiments. 

Our  efforts  over  the  past  two  years  were  both  experimental  and  theoretical. 
The  experimental  work  involved  a  complete  ana  Lysis  of  our  older  data  on  the 
PMMA/DEP  polymeric  system  so  as  to  determine  the  adequacy  of  the  data.  We  made 
the  necessary  repairs  to  the  R-16  instrument  to  improve  it  so  that  adequate  data 
for  the  polymer  system  could  be  obtained.  This  involved  instrument  repair  and 
installing  a  new  data  acquisition  system.  The  R-16  instrument  was  modified 
further  to  allow  accurate  measurements  to  be  made  on  low  viscosity  materials 
such  as  the  slurry-fuel  and  typical  lubricants.  In  addition,  the  limited 
temperature  control  system  was  changed  so  that  measurements  could  be  made  at 
reduced  temperatures.  Measurements  were  then  made  on  both  the  Haake  system  and 
the  R-16  instrument  over  a  range  oi  temperatures  for  the  slurry-fuel  and 
lubricant  materials.  Finally,  the  theory  was  tested  and  adapted  to  fit  the 
data.  In  the  next  section,  a  detailed  review  is  given  of  the  accomplishments 
outlined  above.  Here,  these  efforts  are  briefly  reviewed  in  less  detail  to 
give  the  reader  a  feel  for  our  efforts. 
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We  were  unable  to  find  any  data  in  the  literature  that  was  taken  with  high 
enough  time  resolution  to  allow  evaluation  of  the  parameters  of  our 
kinetic-elastic  model.  Thus,  we  turned  to  the  evaluation  of  our  own  data  on 
PMMA/DEP.  The  sets  of  data  for  data  evaluation  were  very  carefully  reduced  to 
remove  noise.  Although  the  result  of  our  analysis  raised  some  questions  as  to 
the  adequacy  of  this  data,  we  felt  that  the  transient  sets  of  data  were 
satisfactory.  The  measurement  of  the  exact  time  delay  from  the  moment  of 
activation  of  the  brake-drive  switch  to  the  shear  stress  signal  is  important  in 
obtaining  the  true  initial  slope  at  zero  time  in  a  stress  build-up  experiment. 
The  response  time  obtained  is  better  than  0.01  second. 

With  the  simplest  form  of  the  kinetic  equation  of  A  *■  B,  no  combination  of 
constants  would  provide  an  adequate  fit.  In  order  to  improve  the  model  so  that 
it  would  fit  all  types  of  data,  two  steps  were  found  necessary.  As  was  done 
earlier  and  described  in  the  review  presented  in  Appendix  D,  the  first 
step  was  to  use  a  more  empirical  form  for  the  elastic  part  of  the  theory. 

Rather  than  picturing  the  problem  as  a  coupled  elasto-viscous  phenomena,  the 
older  view  pictured  the  elastic  contribution  as  an  elastic  response  to  a 
changing  viscosity  basis  as  caused  by  the  kinetic  change.  The  second  step^ 
required  was  to  modify  the  simple  one-step  kinetic  model  of  A+-BtoA+-B<-C. 
In  the  initial  tests  of  this  idea,  quite  satisfactory  fits  to  the  transient  data 
were  obtained. 

A  detailed  reevaluation  of  all  the  equilibrium  data  was  made  as  well  as  an 
in-depth  evaluation  of  one  constant  shear  rate  experiment  at  a  shear  rate  of 
10.8  sec“l. 

A  preliminary  estimate  of  the  kinetic  constants  for  the  experimental  run 
had  to  be  made  so  that  there  would  be  a  good  starting  point  for  the  optimization 
of  the  constants.  A  study  of  the  effect  of  the  variation  of  the  various 
constants  of  the  kinetic-elastic  approach  on  the  final  response  curve  was 
probably  one  of  the  most  interesting  results  of  the  analysis.  These  results 
provided  the  guidelines  necessary  for  the  optimization  by  iteration.  The 
optimization  of  constants  was  not  done  by  computer  optimization  procedures 
because  of  the  large  number  of  constants  involved.  Rather  the  results  of  the 
previous  parametric  study  were  used  to  adjust  the  theoretical  curve  to  obtain 
the  best  fit  to  the  data.  This  was  accomplished  by  human  iteration  using  the 
computer  graphics  system.  The  predicted  data  from  the  kinetic-elastic  theory 
deviated  approximately  3.95%  from  the  transient  experimental  data  for  the 
constant  shear  rate  experiments.  The  value  of  the  parameters  that  were  varied 
generally  increased  over  the  range  from  3.4  to  21.5  sec-^. 

For  the  stress  relaxation  experiments,  various  reasonable  assumptions  for 
the  initial  conditions  of  the  shear  rate  did  not  affect  the  elastic  equation 
prediction.  The  normal  assumption  of  no  shear  rate  during  the  entire  stress 
relaxation  was  adequate.  However,  a  combination  of  thixotropic  viscosity  and 
the  elastic  equation  was  necessary  to  account  for  the  elastic  contribution  to 
the  stress  relaxation.  An  improved  prediction  of  stress  decay  was  obtained  by 
modest  changes  of  the  kinetic  constants.  It  was  concluded  either  that  the  data 
available  were  not  good  or  that  the  elastic  equation  of  the  theory  for  stress 
relaxation  was  not  exactly  correct. 
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The  initial  shear  rate  for  the  constant  shear  stress  experiments  should  be 
less  than  0.001  sec~l  and  the  value  of  the  second  modulus  elastic  parameter  is 
much  smaller  than  the  first.  The  modified  kinetie'-elast ic  theory  shows 
excellent  agreement  between  the  experimental  and  predicted  shear  strain  data  for 
the  constant  shear  stress  experiments.  The  average  percent  standard  deviation 
from  the  experimental  strain  data  was  approximately  1.80%  from  the  constant 
shear  stress  experiments. 

The  first  elastic  modulus  parameter  cannot  be  considered  a  constant  because 
it  is  affected  by  the  level  of  shear  rate.  The  modulus  parameter  for  stress 
relaxation  is  different  from  that  for  stress  growth.  The  second  parameter  has 
no  significant  dependency  on  shear  rate  for  the  constant  shear  stress 
experiment. 

In  general,  we  can  conclude  that  a  two-step  reaction  concept  of  thixotropy 
coupled  with  a  changing  viscosity  basis  for  the  elastic  response  shows  excellent 
agreement  with  steady  shear  rate,  constant  shear  rate,  and  constant  shear  stress 
i xperiments.  The  tl eory  predicts  well  the  stress  overshoot  at  high  shear  rates 
and  unsteady  shear  i ate  variation  at  constant  shear  stress.  The  present  theory 
cannot  use  a  constant  vaLue  of  viscosity  (q  0)  in  the  elastic  equation  to 
describe  stress  rel.  xation  experiments.  However,  the  value  of  0t  in  t he 
elastic  equation  with  slightly  modified  values  of  the  kinetic  parameters 
predicts  well  the  shear  stress  r  Taxation  results. 

As  a  result  of  further  analysis  of  the  kinetic  theory,  we  discovered  that 
there  coul  i  be  a  mote  va , Id  int.  rpretat ion  of  thi  upper  Newtonian  viscosity  than 
uu  have  male  it  the  past.  In  g*  ieral,  we  had  previously  taken  the  upper 
Newtot  ian  visit  sity  as  t  .e  solve  it  viscosity  (in  this  case  0.09  poise). 

However,  nulecular  heor  of  vis  ous  polymer  solutions  shows  that  the  upper 
Newtonian  visco  ity  is  ptobahly  inch  higher  man  this  (in  the  range  of  btiO  to 
12,0Un  poi  -.e  for  ou  solutions).  Even  with  this  major  change,  the  data  cannot 
be  fitted  vith  the  irapl  one  si  p  kinetics;  thus,  the  more  tomplicated  two  step 
kinetics  are  still  'equiied. 

It  may  be  that  the  adequate  fit  obtained  from  the  two-step  model  is  a 
consequence  of  the  model  having  more  adjustable  constants  rather  than  being  a 
hetter  representat i(  n  of  the  kinetic  change  mechanism.  Then-  is,  for  example, 
an  inconsistency  In  that  we  use  nQ  In  the  stress  growth,  but  must  use  nt 
in  the  stress  relaxation.  We  were  not  sure  at  this  point  if  the  theory  or  the 
polymer  data  was  the  problem.  Tims,  we  undertook  two  steps.  First,  since  the 
polymer  data  being  used  was  the  noisy  data  obtained  several  years  ago  by  Song, 
we  obtained  a  completely  new  set  ot  polymer  data  as  described  in  the  next 
section.  The  new  data  was  quite  good  and  no  noise  removal  (filtering)  was 
needed.  The  second  step  was  a  soul-searching  analysis  of  the  kinetic-elastic 
theory. 

Why  should  the  two-step  model  be  necessary  when  the  one-step  concept  is 
simpler  and  has  fewer  adjustable  constants?  At  this  point  we  made  a  concession 
to  history  and  our  own  work.  We  had  been  forcing  the  first  elasticity  parameter 
to  be  a  function  of  shear  rate  level  only.  Historical ly  it  has  been  taken  as  a 
spectrum  of  values  at  any  one  shear  rate,  i.e.,  a  range  ot  values  that  come  into 
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play  at  various  times.  This  is  in  agreement  with  our  earlier  work  which  showed 
that  the  parameter  determined  from  oscillatory  data  is  both  a  function  of  shear 
level  and  frequency  of  oscillation.  We  had  not  used  this  concept  because  we 
thought  the  parameter  could  be  a  constant  and  that  the  variation  observed  could 
be  a  result  of  the  variation  of  Hc.  In  this,  we  were  incorrect.  The 
modification  to  our  one-step  theory  is  almost  trivial;  i.e.,  the  modulus 
parameters  are  now  real  parameters  and  as  such  are  determined  so  as  to  give  a 
perfect  fit  for  the  one-step  model.  The  test  comes  from  the  simplicity  of 
representation  of  the  parameters  and  from  comparison  of  values  obtained  from 
other  independent  experiments.  The  results  of  this  new  approach  are  very 
promising.  The  values  of  the  constants  for  the  representation  of  the  parameters 
were  determined  independently  foe  each  shear  rate.  Therefore,  if  any  of  them 
varied  with  the  shear  rate  that  variation  would  be  clearly  observed.  Only  three 
of  the  constants  showed  a  systematic  variation:  the  forward  specific  rate 
constant,  the  initial  (or  primary)  value  of  the  elastic  parameter,  and  the  final 
(or  asymptotic)  value  of  the  elastic  parameter.  Further  investigation  of  the 
effecc  of  variations  in  the  elastic  part  of  the  theory  must  await  future  work 
under  the  renewal  of  the  contract.  The  details  of  the  polymer  data  acquisit  •'  on, 
analysis,  and  the  new  developments  are  included  in  the  next  section.  In 
addition,  suggestions  for  furthe  •  work  are  included. 


The  flow  behavi jr  of  suspen  ions  encompasses  the  fulL  range  jf  material 
characteristics.  They  can  vary  in  consistency  from  fluids  to  almost  solids,  and 
they  can  be  Newtonia 1  or  noa-New  onian.  Dilute  suspensions  often  behave  as 
simple  Newtonian  liquids.  As  th.  concentration  o:  particles  increases, 
suspensions  become  r  m-Newtonian .  Thixotropy,  sh> ar-thinning,  rheopexy,  and 
e/e n  viscoelasticity  can  be  seen  in  such  systems.  At  certain  consent  rat  ions , 
suspensions  may  even  function  as  solids.  In  order  to  handle  suspensions 
properly  during  processing  or  pui  ping,  their  flow  behavior  should  be  well 
understood  and  ciiara  terized. 


Carbon-black  has  been  widel.  used  as  a  filler  or  reinforcing  agent  in 
rubbers,  paints,  and  other  polymers.  Though  carbon  black  can  modify  the 
properties  of  polymers,  like  smooth  low-swell  extrusion,  it  also  changes  the 
rheological  properties  of  the  system.  The  change  of  rheological  properties 
depends  on  the  fundanental  properties  of  the  carbon  additive:  particle  size, 
porosity,  and  structure.  Besides  being  used  as  a  filler  in  polymers,  to 
increase  the  available  energy  per  unit  mass  as  a  solid  fuel,  carbon  black,  a 
combustible  solid,  has  been  added  to  a  liquid  fuel,  JP-10,  to  form  a  higher 
energy  slurry  fuel. 


Though  the  presence  of  solids  increases  the  available  energy,  the  flow 
behavior  of  the  slurry  is  made  complicated  by  the  solids.  The  slurry  is  no 
longer  a  simple  Newtonian  fluid  in  spite  of  having  a  Newtonian  solvent,  JP-lU. 
Its  rheological  properties  become  difficult  to  characterize  as  a  result  of  the 
combination  of  solid  and  liquid  properties.  In  order  to  add  to  our  knowledge 
about  the  slurry  fuels,  their  rheological  properties  (viscosity,  elasticity,  and 
plasticity)  were  studied  and  analyzed  in  detail.  These  details  can  be  found  in 
the  next  section  and  are  briefly  cited  here. 


Slurry  fuels  show  non-Newtonian  fluid  behavior  under  imposed  forces.  The 
shear  stress-shear  rate  relationships  were  studied  as  a  function  of  time  over  a 
wide  range  of  shear  rates  and  temperatures.  Since  the  deviation  from  Newtonian 
behavior  is  due  to  the  presence  of  solids,  the  non-Newtonian  characteristics  are 
highly  dependent  on  the  solid  properties:  type  of  solid,  particle  size, 
concentration,  shape,  and  distribution.  No  effort  was  made  to  relate  the 
rheological  characteristics  to  the  properties  of  solids  in  this  preliminary 
study.  Instead,  slurries  with  different  particle  sizes  and  stabilizers  (as 
supplied)  were  studied  to  provide  general  insights  into  the  rheological  nature 
of  such  materials.  These  should  be  useful  to  the  researcher  to  help  design  the 
proper  formulation. 

In  order  to  have  a  better  understanding  of  the  flow  mechanism  of  the 
slurries,  both  transient  and  steady  state  conditions  were  measured.  Transient 
data  was  obtained  b  '  suddenly  imposing  a  step  change  of  shear  rate  on  the  system 
and  measuring  the  suear  stress  as  a  function  of  time.  It  was  found  that  slurry 
fuels  possess  a  time-dependent  behavior  called  thixotropy.  Elasticity,  another 
property  that  can  be  seen  in  suspensions  due  to  tne  interactions  of  particles, 
was  estimated  from  -he  transient  data. 

YieLd  stress  i;  still  another  pronounced  property  of  suspensions.  The 
yield  stress  must  b.  overcome  b;  any  pumping  system  and  will  dictate  the  maximum 
lower  required.  Th  •  yield  stre  ses  of  the  slurry  fuels  were  determined  and 
studied  to  see  how  hey  affect  the  flow. 

.he  change  of  the  v  scosity  of  uels  with  temperature  is  important,  especially 
. ; t  low  temperatures  Viscosity  is  due  to  the  internal  friction  between 
molecules  in  liquid;.  Since  temperature  is  a  measurement  of  molecular  motion, 
temperature  is  the  most  important  variable  affecting  the  viscosity  of 
suspensions.  As  th  :  temperature  drops,  the  solid  properties  of  the  system  could 
dominate  and  be  cri  ical  for  pumping.  In  practice,  the  viscosity  of  fuels 
should  be  as  little  temperature-dependent  as  possible  so  as  to  minimize  the 
required  power  at  low  temperatures.  The  viscosity  of  the  slurries  was  studied 
over  a  wide  range  of  temperatures. 

The  final  objective  was  to  construct  a  rheological  equation  of  state  that 
could  be  used  with  the  equations  of  motion  to  predict  complicated  flow  behavior 
of  slurries  in  practical  situations.  To  accomplish  this,  the  kinetic-elastic 
model  was  modified  to  describe  the  complex  behavior  of  suspensions.  The  fit  of 
the  model  to  the  equilibrium  dai.a  is  excellent.  The  fit  in  the  time  domain  is 
not  as  precise,  but  is  adequate.  Certainly  it  could  be  improved  with  increased 
knowledge  of  the  mechanism  of  the  initial  break-down  of  the  solid  structure.  It 
is  this  study  that  would  be  a  part  of  our  renewal  contract  and  is  given  in  more 
detail  in  the  next  section. 

Besides  the  collection  and  analysis  of  data  for  the  polymer  system  and  the 
in-depth  data  on  the  slurry  fuels,  additional  data  was  obtained  on  the  reference 
Lubricant  5P4E.  This  data  has  been  obtained  over  a  wide  range  of  temperatures 


and  because  of  Its  polymer  like  nature  should  be  interpreted  similarly.  The 
data  set  as  well  as  the  future  work  necessary  to  analyze  this  data  is  detailed 
in  the  next  section. 
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DETAILS  OF  RESEARCH  RESULTS 


I.  POLYMERIC  MATERIAL 


We  were  unable  to  find  any  data  in  the  literature  that  was  taken  with  high 
enough  time  resolution  to  allow  evaluation  of  the  parameters  of  our 
Kinetic-Elastic  model.  Thus,  we  turned  to  the  evaluation  of  our  own  data  on 
PMMA/DEP. 

A.  Instrument  Response 

Since  the  elastic  modulus  parameter  of  the  Kinetic-Elastic  model  is  related 
with  the  initial  slope  at  zero  time,  it  is  important  to  establish  the  exact  time 
delay  from  the  moment  of  activation  of  the  break-drive  unit  of  the  Weissenberg 
rheogoniometer  to  the  appearance  of  the  shear  stress  signal.  Thus,  the 
measurement  of  the  mechanical  response  of  the  Model  R-16  Weissenberg 
rheogoniometer  was  undertaken  for  analyzing  the  transient  data  during  the 
initial  period  of  stress  growth.  The  resulting  data  and  analysis  allowed  us  to 
obtain  the  time-delay  to  the  true  beginning  in  a  stress-buildup  experiment.  The 
response  time  obtained  was  much  faster  than  that  suggested  in  the  instruction 
manual  (about  0.04  seconds).  The  response  times  of  the  viscometer  at  constant 
shear  rates  of  3.408  sec-*-,  10.8  sec”'  and  21.48  sec~^  are  0.009  sec, 

0.007  sec,  and  0.008  sec  respectively. 

B.  Noise 

The  sets  of  data  for  data  evaluation  were  very  carefully  reduced  so  as  to 
remove  any  extraneou ;  noise  introduced  from  the  gearbox  system.  We  feel  i hat 
the  transient  sets  of  data  are  satisfactory. 

C.  Theory  for  One-Step  Reaction 

Appendix  B  of  the  proposal  gives  a  complete  tabulation  of  the  equations  of 
the  modified  Kinetic-Elastic  theory  used  in  the  initial  theory  evaluation.  The 
necessary  equations  from  that  Appendix  are  reproduced  here  as  Appendix  B. 
Equations  1,  6,  and  7  are  the  kinetic  rate  expressions.  Equation  11  is  the 
elastic  contribution.  Equations  12  to  14  are  used  to  evaluate  the  equilibrium 
constants.  Equations  4,  5,  and  8  are  used  for  fitting  of  the  stress  growth 
data. 


Equilibrium  Evaluation  for  One-Step  Reaction 


The  evaluation  of  the  equilibrium  kinetic  constants  (K  and  p)  of  the 
Kinetic-Elastic  model  has  been  done  at  the  equilibrium  state.  These  and  other 
constants  of  the  theory  are  noted  in  Appendix  B  of  this  report  and  in  more 
detail  in  Appendix  B  of  the  proposal.  In  this  analysis,  4  combinations  of  m,  n, 
and  a  were  used  (m  ■  1,  n  “  2  or  3,  a  ■  1  or  3.5).  The  results  show  that  all 
combinations  of  m,  n,  and  a  give  reasonable  prediction  of  the  shear  rate  and 
equilibrium  shear  stress.  The  values  obtained  for  p  and  K  are  listed  in  Park's 
thesis  (75)*.  The  maximum  range  of  %  deviation  of  the  predicted  shear  rate  and 
equilibrium  stress  are  19%  and  16%,  respectively.  The  equilibrium  constants  (p 
and  K)  show  a  weak  dependency  on  concentration.  The  p  values  are  slightly  more 
sensitive  to  concentration  change  than  the  K  values. 


*  Specific  thesis  c  ted  herein  are  given  in  the  proposal  reference  list 
reproduced  here  as  Appendix  A. 
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tit.  the  transient  dit.»  was  then  ased  t  i  tit  t  ne  entire  mrvt.  1  ne  initial 
slope  was  used  to  establist  t  tie  elastic  parameter.  various  parameter  search 
techniques  were  used  t.  establish  t  tie  remaining  constants.  A  typical  example  is 
shown  in  Figure  1. 

With  the  simplest  iorm  oi  the  kinetic  equation  of  As-  B  as  described  in  the 
original  proposal,  no  combination  of  constants  would  provide  an  adequate  fit. 

If  the  time  to  the  peak  in  the  response  curve  were  matched,  the  degree  of 
overshoot  was  excessive.  It  the  amplitude  of  the  over-shoot  was  matched,  the 
time  to  the  peak  was  tor  short.  The  poor  initial  fit  is  a  result  of  the 
non-flexible  nature  of  the  model,  i.e.  ,  logically  sa“ isfactory  but  not 
satisfactory  for  an  empirical  fit. 


D.  Two-Step  Reaction 

In  order  to  improve  the  model  so  that  it  would  fit  all  types  of  data,  two 
steps  were  necessary.  The  first  step  was  to  use  a  aore  empirical  form  for  the 
elastic  part  of  tht  theory  as  was  ione  earlier  and  is  described  in  the  review  of 
our  proposal.  Rather  than  picturing  the  problem  as  a  coupled  elasto-viscous 
phenomena,  we  used  the  older  view  which  pictured  the  elastic  contribution  as  an 
elastic  response  t<  a  changing  viscosity  basis  as  caused  by  a  kinetic  change. 

The  se-ond  step  re<  aired  was  to  modify  the  simple  one-step  kinetic  model  of 
A  t  B  to  something  more  complex  lice  A  i  B  Z  C.  This  change  makes  evaluation  of 
the  constants  more  difficult.  All  of  the  variables,  constants,  and  parameters 
used  in  the  sequel  are  given  in  Tible  I  and  have  the  same  parallel  meaning  as 
the  corresponding  terms  have  in  the  one  step  reaction;  e.g.,  KA  and  Kg  are 
the  two  equilibriur  constants  that  are  now  needed  to  replace  K;  m^,  m2, 
n^,  n2  are  reactioi  orders  needed  to  replace  m  and  n,  etc. 


Equilibrium  Evaluation  for  Two-Step  Reaction 


The  value  of  a  was  selected  as  3.5  based  on  previous  work  and  not 
modified  or  adjusted  in  any  way.  The  value  of  nw  was  taken  as  0.09  poise  bused 
on  experimental  analysis  and  not  modified  or  adjusted.  Six  combinations  of 
ra^ ,  m2,  n^,  and  n2  were  selected  for  testing.  These  six  were  selected 
based  on  logical  arguments  for  the  orders  of  the  forward  and  reverse  reactions. 
The  variation  of  the  parameter  S  was  studied  for  each  of  the  six  cases.  (The 
parameter  S  is  defined  in  equation  (L)  and  (M)  in  Table  I).  The  specific 
evaluation  of  the  constants  was  done  by  first  using  equation  (E)  to  obtain  KA 
and  pA  from  a  log-log  plot  of  the  LHS  (left-hand-side)  against  the  equilibrium 
shear  stress  Teg.  The  constants  Kg  and  pB  were  then  evaluated  in  a 
similar  manner  using  equation  (F)  and  the  values  of  KA  and  pA  evaluated  in 
the  previous  step.  Values  of  S  below  2  resulted  in  F  being  negative,  which  is 
unrealistic.  Too  high  values  of  S  reverted  the  two  step  kinetics  to  the  older 
single  step  system.  Since  the  selection  of  S  was  not  critical  to  fit  the  data, 
a  value  of  2.5  was  selected  for  further  evaluation  of  the  data.  All  six  cases 
for  S  =  2.5  were  quite  satisfactory  and  within  the  scatter  of  the  experimental 
equilibrium  data. 


From  past  wor*  on  equilibrium  data,  we  have  learned  that  data  at  low  rites 
of  shear  should  be  weighted  less  than  that  at  higher  rates.  This  is  a  result  of 
Fa  being  near  unity.  When  the  ter  a  that  involves  1  -  Fa  ts  evaluated,  we 
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TABLE  I 
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Figure  2.  Evaluation  of  K.  and 


introduce  considerable  error  due  to  the  taking  of  differences  of  large  numbers. 
Thus,  the  equilibrium  results  were  evaluated  by  using  a  variable  weighting 
factor  to  improve  fit.  To  illustrate  the  adequacy  of  the  data  fit.  Figure  2 
gives  the  evaluation  of  KA  and  pA  from  the  plot  according  to  equation  (E) 
and  Figure  3  gives  the  evaluation  of  Kg  and  Pg  from  the  plot  according  to 
equation  (F).  These  are  presented  only  for  case  1,  but  all  other  cases  were 
nearly  as  good.  Finally,  Figure  4  gives  the  final  backcalculated  result  of  the 
experimental  and  predicted  basic  shear  diagram  equilibrium  results  for  this  same 
case.  Clearly  these  results  are  more  than  adequate  to  represent  the  data  over  a 
range  of  concentrations. 

D.2  Transient  Evaluation  of  Two-Step  Reaction  for  Stress  Growth 

As  was  seen  in  the  previous  section,  and  as  was  true  with  the  one-step 
model,  no  problems  were  experienced  in  fitting  the  equilibrium  data.  Completely 
adequate  fits  were  obtained  in  all  cases.  However,  the  fits  to  the  complete 
transient  curves  were  only  somewhat  improved  over  those  shown  in  Figure  1  (see 
Figure  5).  This  result  strongly  suggested  chat  the  kinetic  aspect  of  the  theory 
was  not  the  main  caufe  of  the  equation's  failure  to  provide  an  adequate  fit  for 
the  data. 

In  addition  to  tne  kinetic  c  lange,  a  change  in  the  elastic  response 
ormulation  of  the  theory  was  mad*.  Rather  than  picturing  the  problem  as  a 
■  oupled  elasto-vi scou  ;  phenoaena,  an  older  view  of  the  elattic  contributicn 
ieir  an  elastic  resp >nse  to  a  ch  inging  vis  osity  basis  as  a  ised  by  the  kinei  Lc 
'■han  ;e  was  used.  In  the  initial  ests  of  t  lis  idea,  quite  sai isfactory  fits  to 
the  transient  dati  were  obtained  ^see  Figure  6). 

Preliminary  values  for  kinetic  constants.  A  preliminary  estimate  of  all 
the  kinetic  constants  was  made  so  that  there  would  be  a  good  starting  point  for 
the  optimization  of  constants. 

The  procedure  for  evaluation  of  k^  and  P ^  is  first  to  calculate  FA 
from  the  data  using  equations  (G)  and  (I)  and  then  to  plot  this  FA  as  a 
function  of  time.  The  calculation  can  be  done  since  Gj  will  have  been 
previously  estimated  from  the  initial  slope  of  the  data  and  n0  is  known. 

From  the  plot,  the  initial  rate  cm  be  obtained  from  the  slope  at  t  =  0,  i.e., 
an  evaluation  of  the  term  -  (dFA/lt) t„Q.  Equation  (A)  at  time  zero  reduces 
to  equation  (J),  whic  \  in  turn  ca  i  be  rearranged  to  equation  (K)  by  using 
initial  conditions.  A  plot  of  th  ;  LHS  of  equation  (K)  versus  Y  will  give  a 
sio^e  of  PJ,  and  from  the  intercept,  the  concentration,  and  H(),  kg  can  be 
calculated.  Since  P2  and  k£  are  uniquely  related  to  pj,  kj,  KA,  aid 
PA,  k2  an<*  P2  are  als  i  known  (for  k2  evaluation  it  is  assumed  that 
kg  =  0).  The  constant  kg  appears  in  the  rate  equations  in  (A)  and  (ti)  of 
Table  I.  It  describe;  the  Browni  ig  motion  buiLd-up  of  structure  as  is  explained 
in  Appendix  B  of  the  proposal.  Tie  term,  kg,  is  small  compared  to  the  term  to 
which  it  must  be  adde  1.  Since  da  a  at  various  shear  rates  were  used  in  this 
evaluation,  values  were  not  obtal  led  that  are  a  function  ot  shear  rate.  The 
value  for  pi  was  0.38^9  and  that  "or  k^  was  0.2204.  These  values  will  serve 
as  the  initial  estima.es  for  optimization  at.  all  shear  rates. 

Evaluation  of  th *  preliminary  values  for  P4  and  kA  was  more  difficult 
because  data  near  equ Llibrium  had  to  be  used.  Equation  (D),  with  i  ficed  value 
of  1>  of  2.3,  was  ased  to  obtain  d  g/dt  from  the  values  of  dr'-,  dt  previously 
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Figure  6.  Comparison  of  experimental  data  to  complex  reaction  and 
elast ic— r esponse  to  a  changing  viscosity  of  kinetic- 

elastic  model  for  one  set  of  constants. 


calculated  from  equations  (G)  and  (I).  After  some  rearrangement  of  equations 
(A),  (B),  and  (D)  as  well  as  the  introduction  of  the  relation  between  X4, 
k.5,  and  Kg  and  of  that  between  P4,  P5,  and  pg,  equation  (L)  can  be 
obtained. 

After  further  multiplication,  this  can  be  rearranged  into  a  more  convenient  form 
given  by  equation  (M).  A  plot  of  the  LHS  of  equation  (M)  versus  the  shear 
s'tress  provides  for  the  evaluation  of  P4  and  £4,  from  which  P5  and  k.5 
can  be  readily  obtained  with  the  assumption  of  kg  =  0.  This  is  clearly 
parallel  to  taking  k3  ■  0  in  the  same  equation.  All  cases  were  evaluated. 

Parametric  analysis  of  equations.  A  study  of  the  effect  of  the  variation 
of  the  various  constants  of  the  kinetic-elastic  approach  on  the  final  response 
curve  produced  one  of  the  most  interesting  results  of  the  analysis.  A  figure  is 
needed  to  show  each  constant.  Figure  7  shows  the  effect  of  KA;  the  data  curve 
for  run  7/322  is  shown  for  comparison.  Figure  8  shows  the  effect  of  that 
changes  both  the  initial  slope  and  the  peak  overshoot  value.  The  figures  for  the 
effect  of  the  other  parameters  can  be  found  in  reference  74.  These  results 
provide  the  guideline  .  necessary  lor  the  optimization  by  iteration  to  be 
described  next. 

Optimization  of  tonstants.  "he  optimization  of  constants  was  not  done  by 
computer  optimization  procedures  because  of  the  large  number  of  constants 
involved.  Father  the  results  of  the  previous  parametric  study  were  used  to 
adjust  the  theoretica  curve  so  a:  to  obtain  the  best  fit  to  toe  data.  Th  _s  was 
accomplished  by  human  iteration  u.-  ing  the  computer  graphics  system.  The 
specific  procedure  was  first  to  plot  the  experimental  data.  The  equilibrium 
constants  were  then  ai justed  to  provide  the  best  fit  for  the  specific  run.  This 
involved  only  a  small  change  in  pA;  KA,  Kg,  and  pg  remained  as 
previously  estimated.  Next  G^  was  adjusted,  if  necessary,  to  give  the  best 
fit  for  the  initial  slope.  Finally  pj  and  P4  were  slightly  adjusted  to 
provide  the  best  final  fit  for  the  system.  The  constants  k^  and  k4  were  not 
changed  from  their  ini tial  estimates.  Recall  that  the  trend  in  changing  them 
was  iiKactly  the  same  as  changing  p^  and  P4.  As  a  test  on  the  adequacy  of 
the  fit,  a  mean  squared  deviation  given  as  equation  (N)  in  Table  I  was  used. 

The  result  is  shown  in  Figure  9.  The  final  values  for  all  constants  and  all  of 
the  plotted  results  can  be  found  in  reference  74.  The  standard  deviation  of  all 
of  the  data  was  less  t  han  4%. 

0.3  Transient  Evaluation  of  Two-Step  Reaction  for  Stress  Relaxation 

The  constants  from  the  stress  growth  experiment  were  used  to  predict 
stress  relaxation.  A  constant  viscosity  (ne,j),  which  poorly  predicted 
results,  was  used  in  the  elastic  equation,  Much  better  results  were  obtained 
when  we  replaced  nec.  with  a  variable  viscosity  (n^).  The  equations  used  are 
gLven  in  Table  I.  We  found  that  using  a  constant  viscosity  was  better  for 
stress  growth  and  a  variable  viscosity  was  better  for  stress  relaxation.  We  are 
uncomfortable  with  th: s,  but  have  no  explana!  Lon  for  the  results. 

The  results  of  t!  e  parametric  study  show  that  the  stress  decay  increases  .it 
a  faster  rate  with  an  increase  in  Gi  (see  Figure  10).  pA  has  m  effect  on 
the  initial  slope.  Gj  controls  the  initial  slope.  Pi  has  the  same  effect 
as  pA.  P4  has  an  effect  on  the  curve  near  the  equilibrium  state.  The 
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higher  P4,  the  higher  the  tail  part  of  the  curve.  However,  P4  has  no 
significant  effect  within  a  range  of  -0.29  to  -0.14.  K4  has  the  same  effect 
as  p4,  but  no  significant  effect  within  the  above  range. 

Use  of  the  maximum  values  for  the  parameters  that  were  observed  for  stress 
growth  provided  a  reasonable  prediction  for  the  stress  decay.  Therefore,  the 
equilibrium  and  kinetic  constants,  except  p^,  were  retained  and  the  value  of 
Gj  was  used  as  0.5  x  10^  dyne/cm^  for  all  the  data  sets.  Here,  p^  was 
changed  within  a  range  of  0.26  to  0.30  for  the  best  fit. 

The  standard  deviation  was  used  as  a  fitting  criteria;  the  average  value 
was  0.05.  Comparison  of  one  of  the  predicted  results  to  the  experimental  data 
is  illustrated  in  Figure  11. 

Effect  of  Initial  Conditions.  Previously  we  assumed  that  the  shear  rate 
dropped  instantaneously  to  zero  at  the  initial  time  of  the  stress  relaxation 

experiment.  Since  this  is  not  exactly  true,  we  realized  that  the  initial  shape 

of  the  relaxation  curve  could  be  affected.  Thus,  we  modified  the  analysis  to 
allow  a  nonzero  initial  shear  rate  condition.  Due  to  a  response  delay  in  the 
Weissenberg  rheogoniometer,  the  lower  plate  of  the  rheogoniometer  continues  to 
rotate  for  a  very  short  time  after  the  brake  switch  for  the  lower  plate  is 
activated.  To  approximate  this  alternate  initial  condition,  the  shear  rate  was 
expressed  as  a  function  of  time  rither  than  being  set  to  zero.  The  resulting 
set  of  equations  that  must  be  simultaneously  integrated  is  the  same  with  only 
the  first  equation  ft-  the  viscoelastic  part  of  the  theory  requiring 
modification. 

Since  the  brake- drive  system  response  time  of  the  rheogoniometer  was  within 
0.01  seconds,  the  value  of  8  should  be  between  100  and  150.  The  term  $  is 

defined  in  equation  (Q)  in  Table  i.  The  result  of  the  parametric  study  of  the 

8  -effect  showed  that  a  value  of  3  between  25  to  200  had  no  significant  effect 
during  stress  relaxation. 

The  normally  assumed  initial  condition  of  zero  shear  rate  is  clearly 
satisfactory  for  the  stress  relaxation  experiment. 

D.4  Transient  Evaluation  for  Two-Step  Reaction  for  Shear  Change  at  Constant 
Stress 


The  shear  strain  data  during  the  constant  shear  stress  experiments  is 
not  contaminated  by  noise,  since  the  motor  drive  system  is  disengaged.  Thus, 
the  raw  data  obtained  by  Song  was  directly  analyzed. 


Since  the  consta 
(G2),  but  no  new  para 
concerned,  the  equili 
constant  shear  rate  ( 
obtained  at  CSR  for  a 
Song  was  shear  strain 
numerically  integrate 
this  procedure  is  tha 
numerical  differentia 
comparison  of  shear  r 


it  stress  experiments  introduce  one  new  elastic  parameter 
meters  as  far  as  the  kinetic  part  of  the  theory  is 
orium  shear  stress  was  carefully  compared  to  those  at 
CSR)  in  order  to  apply  values  of  the  kinetic  parameters 
nalysis  of  shear  strain  data.  Since  the  raw  data  taken  by 
instead  of  shear  rate,  the  elastic,  equation  was 
1  again  to  obtain  the  shear  strain.  Another  reason  for 
t  the  numerical  integration  involves  Les:  error  than 
tion,  which  tends  to  amplify  small  variat ions.  Therefore, 
nte  was  made  indirectly  by  strain  curves. 
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of  the  shear  stress  relaxation 
;  predicted  results  (SHI.AK  012) 


The  values  01  the  parameters  obtained  from  t  lie  cons  taut  -.hear  rate 
experiments  (having  the  same  equilibrium  shear  stress  as  those  at  I  l.e  constant 
shear  stress  experiments)  were  directly  used  for  the  jn.i  Lysis  of  transient  shear 
strain  and  shear  rate  during  the  first  5  seconds  of  the  experiment. 

The  results  of  a  parametric  study  show  that  the  initial  condition  of  the 
shear  rate  (  Yn  at  t  -  0)  has  a  shifting  effect  on  the  origin  (ttru  time)  and 
that  y/y  n  t_Q  should  be  less  than  0.01  sec-'.  The  elastic  parameter 
has  an  effect  on  curvature  of  shear  strain.  The  larger  the  value  of  the 

more  convex  the  strait  curve.  Values  of  G2  less  than  luo  have  no  significant 
effect . 

The  elastic  parameter  (Gj)  is  the  only  parameter  that  /.as  ur.ee.  to  control 
the  strain  curve  and  vis  found  to  be  approxit  atelv  l  lJ  dyne  for  the  best 

fit.  This  value  of  G^  is  much  smaller  than  loose  obtained  oy  other 
researc  hers.  Even  though  the  decreased  vaiut  of  1-2  gets  closer  to  tiie 
experimental  shear  st  ain  curve,  the  result  snows  that  the  i'i.nil  value  of  Ga 
obtained  for  each  stress  level  has  no  depend*  icy  on  the  shear  stress. 

The  prediction  ol  the  shear  strain  with  the  model  was  1  oir.,>areU  to  the 
experimental  data  and  the  percent  standard  dt  /iation  (PSTD)  tor  the  entire  data 
set  was  about  1.2/%.  in  order  to  show  the  e>  elleiit  agreement  between  the 
experimental  and  predicted  shear  strain  data,  a  typical  plot.  i.  shewn  in  Figure 
12. 


E.  Conparison  of  Kinc  tic-Elasticy  Two-Step  Ri  iction  Results  with  Earlier  Efforts 


The  theory  has  beta  analysed  lor  stress  growth,  stress  relaxation,  constant 
shear  stress  experimei  ts  and  oscillation  experiments  by  earlier  researchers. 
These  previous  researchers  have  used  the  simple  one-step  reaction  concept  of 
A  ^  B  for  the  kinetic  part  of  the  theory.  However,  in  this  work,  the  old  theory 
had  to  be  modified  by  introducing  3  two-step  reversible  reaction  (A  t-  b  t  C) 
concept.  Therefore,  the  values  of  the  constants  and  parameters  from  the  present 
work  will  be  different  from  those  obtained  by  others. 


The  previous  researchers  provided  values,  of  an  elastic  modulus  parameter, 
Gj,  for  their  models.  Song  (10)  calculated  Gj  from  the  initial  slope,  but 
he  had  difficulty  in  obtaining  the  true  initial  slope  due  to  the  noise  and  not 
knowing  the  exact  zero  in  time.  Thus,  he  obtained  his  values  of  G j  from  the 
approximated  slope  by  the  simple  window  averaging  technique.  The  values  of  Gj 
obtained  by  Song  were  within  the  range  of  66,000  to  110,000  dyne/cm*,  which 
were  determined  at  the  shear  rates  of  1.08  sec-'  and  21.48  sec-', 
respectively.  His  results  showed  that  Gj  ini  reased  with  an  increase  in  shear 
rate.  This  trend  of  G^  is  in  agreement  with  the  results  of  the  present  work. 

Pandalai  (12)  obtained  the  values  of  G^  by  a  linear  interpolation  for  his 
oscillation  experimen  s  based  on  the  fact  th.  t  the  elastic  parameters  (Gj  and 
G2)  are  functions  of  oncentratiou.  In  his  cork,  the  value  of  Gj  increased 
with  an  increase  in  s dear  rate.  However,  thi  value  of  G^  was  iependent;  on  the 
value  of  k2  used  to  fit  the  initial  slope  in  his  work.  Another  interesting 
observation  by  Pandalai  was  that  a  much  high'  r  value  of  Gj  had  to  be  used  for 
the  stress  relaxation  experiments  to  obtain  a  adequate  fit.  Pandalai's  results 


Tigure  12.  Comparison  of  the  shear  strain  data 
and  predicted  results  (STOSS  052)  . 


differ  from  the  preseat  work  which  shows  a  much  smaller  value  of  Cj  for  the 
stress  decay  compared  with  the  value  for  stress  overshoot. 


Lander  (13)  extensively  studied  the  elastic  modulus  parameters  (G^  and 
G2)  by  oscillation  methods.  The  effects  of  oscillatory  frequency,  shear  rate 
and  temperature  were  all  considered.  He  calculated  Gj  for  five  different 
oscillatory  frequencies  within  the  range  of  0.03  to  3.0  Hz  at  four  different 
temperatures.  In  Lander's  work,  the  values  of  G^  and  G2  were  obtained  by 
extrapolations  to  zero  amplitude  of  oscillation.  He  concluded  that  the  values 
of  increased  with  oscillatory  frequency  and  decreased  with  tempetature. 

The  viscosity  effect  on  Gj  in  his  work  is  similar  to  the  effect  of  the 
oscillatory  frequency.  For  the  shear  rate  effects  on  Gj,  the  value  of  Gj 
goes  to  zero  as  the  upper  Newtonian  region  is  approached.  However,  the  effect 
of  shear  rate  on  Gj  from  Lander's  work  is  contrary  to  those  obtained  from  the 
present  Kinetic-Elastic  model.  Also,  the  value  of  Gj  obtained  by  Lander  in 
the  near  non-Newtonian  transition  area  should  be  reexamined  at  higher  shear 
rates.  Lander  had  difficulty  establishing  the  value  of  Gj  from  his 
oscillatory  tests  at  shear  rates  over  4  sec-*. 

In  the  present  work,  the  values  of  G|  tor  the  constant  shear  rate 
experiments  are  obtained  from  the  initial  slope  of  the  filtered  data.  The  value 
of  G2  is  obtained  from  the  constant  shear  stress  experiments  and  is  the  only 
parameter  to  control  the  shear  strain  curve. 

Tables  II  and  III  summarize  the  values  of  the  elastic  modulus  parameters 
obtained  in  this  work  and  by  other  researchers. 

A  conclusion  from  this  work  and  that  done  by  others  is  that  the  elastic 
modulus  cannot  be  considered  constant,  but  is  a  function  of  shear  rate. 

F .  Conclusions  for  Kinetic-Elastic  Theory  -  Phase  I 


In  general,  we  can  conclude  that  a  two-step  reaction  concept  of  thixotropy 
coupled  with  a  changing  viscosity  basis  for  the  elastic  response  shows  excellent 
agreement  with  steady,  constant  shear  rate,  and  constant  shear  stress 
experiments.  The  theory  predicts  well  the  stress  overshoot  at  high  shear  rates 
and  unsteady  shear  rate  variation  at  constant  shear  stress.  The  present  theory 
cannot  use  a  constant  value  of  viscosity  (0o)  in  the  elastic  equation  to 
describe  stress  relaxation  experiments.  However,  the  value  of  H  t  in  the 
elastic  equation  with  slightly  modified  values  of  the  kinetic  parameters 
predicts  well  the  shear  stress  relaxation  results.  It  is  recognized  that  the 
Kinetic-Elastic  approach  has  great  potential  for  the  representation  of 
rheological  data  for  both  transient  and  equilibrium  states  but  requires  further 
development. 

More  specific  conclusions  are 

1.  The  measui ement  of  the  exact  time  delay  from  the  moment  of 

activatior  of  the  brake-drive  switch  of  the  shear  stress  signal 
is  important  in  obtaining  the  tru  ;  initial  slope  at  zero  time  in 
a  stress  build-up  experiment.  The  response  time  obtained  is  better 
than  0.01  second. 
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Table  TI  Comparison  of  Gj  obtained  by  various  investigators  using 
PMMA  in  DEP. 


Investigators 

Range  of  Gj 
(dyne/cm2 ) 

Remarks 

Present  work 
by  Park 

102,270^272,270 

obtained  with  shear  rates 
of  3.408  to  21.48  sec-1. 

35%  PMMA  in  DE?  only. 

Song 

66,000^110,000 

obtained  with  shear  rates 
of  1.08  to  21.48  sec-1. 

30$,  35$  and  40$  PMMA  in 
DEP. 

Pandalai 

30,000^180,000 

obtained  with  oscillatory 
frequencies  of  0.1563  Hz 
to  1.19  Hz. 

38.5$  PMMA  in  DEP. 

Lander 

13,660^133,800 

obtained  at  zero  shear  rate 
with  the  oscillatory  fre¬ 
quencies  of  0.3  Hz  to  3.0 

Hz. 

30$,  35$  and  40$  PMMA  in 

DEP. 

Table  II.T  Comparison  of  G2  obtained  :>y  various  investigators 
using  PMMA  in  DEP. 


Investigator 

Range  of  G2 
(dyne/ cm2) 

Remarks 

Present  work 
by  Park 

100^1,000 

obtained  with  shear  stress  of 
85,810  to  136,510  dyne/cm2. 

35$  PMMA  in  DEP  only. 

Song 

40,000^700,000 

obtained  with  shear  stress  of 
85,810  to  136,510  dyne/cm2. 

30$,  35$  and  40$  PMMA  in  DEP. 

Pandalai 

12,000^30,000 

obtained  with  oscillatory 
frequencies  of  0.1563  Hz  to 

1.19  Hz.  38.5$  PMMA  in  DEP. 

Lander 

36,700^34,000 

obtained  at  zero  shear  rate 
with  the  oscillatory  frequencies 
of  0.3  Hz  to  3-0  Ha. 

30$,  35$  and  40$  PMMA  in  DEP. 

2.  Separation  of  noise  from  the  fast  response  signal  is  essential  for 
the  analysis  of  the  transient  data  at  constant  shear  rates.  A  fast 
Fourier  analysis  has  been  effectively  used  for  the  noise  removal 
process.  The  net  result  obtained  was  a  clean  response  curve.  The 
contamination  noise  was  approximately  30  Hz  and  probably  is  from 
the  1800  rpm  (30  Hz)  motor. 

3.  With  the  simplest  form  of  the  kinetic  equation  of  A  t  B,  no 
combination  of  constants  provided  an  adequate  fit. 

4.  In  order  to  improve  the  model,  three  alternate  forms  were  introduced: 
1)  The  simple  one-step  kinetics  of  A  J  B,  coupled  with  the  lower 
Newtonian  viscosity  (  n0)  in  the  elastic  response  fails  to  improve 
the  fit  significantly,  2)  A  two-step  reaction  kinetic  concept  of 

A  £  B  £  C,  coupled  with  the  normal  viscoelastic  response  was 
introduced  and  was  only  a  slight  improvement  over^the  simpler  one- 
step  model.  Thus,  it  was  concluded  that  the  kinetic  aspect  of  the 
theory  could  not  be  the  main  cause  of  the  failure  to  provide  an 
adequate  fit  to  the  data.  3)  The  two-step  reaction  kinetics  of 
A  ?  B  ?  C,  coupled  with  the  time  constant  of  for  the 

elastic  response,  was  used  to  obtain  quite  satisfactory  fits  to  the 
transient  data. 

3.  Among  6  combinations  ol  mi,  m2,  n;  and  n2  tested,  m ^  =  1 , 

m2  =  1,  nj  ■  1  and  n2  *  2  (both  fi  cward  rate  orders  are  unity 
as  well  as  the  reverse  rate  for  tl  s  first  reaction,  with  the  reverse 
rate  for  the  second  reaction  bei^  two)  produce  the  best  fit. 

6.  A  new  parameter  S,  which  is  only  valid  at  the  equilibrium  state,  wa » 
introduced  for  the  evaluation  of  the  value  of  Fg  of  the 
intermediate.  A  value  of  2.5  for  5  was  selected  for  further 
evaluation  of  the  data,  since  the  value  of  S  was  not  critical  to  the 
fit  of  the  data. 

7.  It  was  found  that  a  weighting  factor  to  compensate  for  small  values 
of  log  (  Teq)  improved  the  fit.  he  results  were  more  than 
adequate  to  represent  the  data  over  a  range  of  concentrations. 

8.  A  study  of  the  effect  cf  parameters  of  the  Kinetic-Elastic  approach 
for  the  stress  overshoot  experiments  is  an  essential  part  of  the 
optimization  of  the  various  constants.  The  results  of  a  parameter 
analysis  provided  the  guidelines  necessary  for  the  optimization  using 
an  interactive  computer  graphics  system. 

9.  The  present  Kinetic-Elastic  theory  shows  excellent  agreement  with 
both  steady  state  and  transient  data  for  the  constant  shear  rate 
experiments.  For  the  constant  shear  rate  experiments,  the  percent 
standard  deviation  from  the  experimental  data  was  approximately  6.13% 
for  the  equilibrium  state  and  1.58%  for  the  transient  experimental 
data.  The  prediction  of  stress  overshoot  at  the  shear  rate  of 
21.48  sec-*  was  excellent  with  a  percent  standard  deviation  of  only 
1.09%. 


10.  For  the  constant  shear  rate  experiments,  k^  and  k^  have  no 
significant  effects  within  the  range  of  0  to  10  and  0  to  18, 
respectively. 

11.  The  same  elastic  parameters  and  kinetic  constants  (except  k^)  can¬ 
not  be  used  to  predict  the  results  of  the  stress  growth  for  all  the 
constant  shear  rate  experiments.  The  parameters  and  constants  were  a 
function  of  shear  rate.  The  values  of  the  parameters  and  constants 
were  found  to  generally  increase  over  the  range  from  3.4  to  21.5 
sec“l. 

12.  For  the  stress  relaxation  experiments,  various  reasonable  assumptions 
for  the  initial  conditions  of  the  shear  rate  do  not  affect  the 
elastic  equation  prediction.  The  normal  assumption  of  no  shear  rate 
during  the  entire  stress  relaxation  is  adequate. 

13.  A  combination  of  thixotropic  viscosity  and  the  elastic  equation  is 
necessary  to  account  for  the  elastic  contribution  to  the  stress 
relaxation.  An  improved  prediction  of  stress  decay  was  obtained 
by  modest  changes  of  the  kinetic  constants.  It  was  concluded  that 
either  the  data  available  are  not  good  or  that  the  elastic  equation 
of  the  theory  for  stress  relaxation  is  not  exactly  correct. 

14.  The  initial  shear  rate  for  the  constant  shear  stress  experiments 
should  be  less  than  0.  J01  sec”*.  The  value  of  &2  muc^ 
smaller  than  and  is  about  100  hyne/cm^. 

15.  Cone  angle  and  plate  diameter  have  no  significant  effect  on 
constant  shear  stress  »xperiments. 

16.  The  elastic  modulus  parameter  cannot  be  considered  a  constant. 

is  affected  by  the  level  of  shear  rate.  Gj  for  stress  relaxation 
is  not  the  same  as  Gj  cor  stress  growth.  G£  has  no  sig¬ 
nificant  dependency  on  shear  stress  for  the  constant  shear  stress 
experiments. 

17.  The  present  Kinetic-El sstic  theory  shows  excellent  agreement  between 
the  experimental  and  predicted  shear  strain  data  for  the  constant 
shear  stress  experiments.  For  the  constant  shear  stress  experiments 
the  average  percent  standard  deviation  from  the  experimental  strain 
data  was  approximately  1.80%. 

Some  General  Conclusions  are 

1.  A  two-step  reaction  concept  of  thixotropy  coupled  with  a  changing 
viscosity  basis  for  the  elastic  response  shows  excellent  agreement 
with  stiady,  constant  shear  rate,  and  constant  shear  stress  ex¬ 
periments.  The  theory  predicts  well  the  stress  overshoot  at  high 
shear  rites  and  the  unsteady  shear  rate  variation  at  constant  shear 
stress. 
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2.  The  present  theory  cannot  use  a  constant  value  of  viscosity  (H0) 
in  the  elastic  equation  to  describe  stress  relaxation  experiments. 
However,  the  value  of  rit  the  elastic  equation  with  slightly 
modified  values  of  the  kinetic  parameters  predicts  well  the  shear 
stress  relaxation  results. 

G.  Re-evaluation  of  the  Basic  Kinetic-Elastic  Concepts 


Although  the  data  can  be  fitted  with  the  two-step  model  there  are  too  many 
adjustable  constants.  The  need  for  two-step  versus  one-step  kinetics  is  not 
satisfying  and  there  are  some  inconsistencies  in  the  theory.  Thus,  even  though 
the  fit  is  quite  good,  we  made  the  decision  to  initiate  a  second  phase  of  our 
data  analysis  in  which  we  would  retain  the  one-step  model  and  try  to  determine 
why  it  does  not  fit  the  data  adequately.  First,  we  recognized  the 
upper-Newtonian  limit  was  not  the  solvent  viscosity  but  rather  something  much 
larger.  Second,  onLy  one  term  in  one  equation  needed  revision  to  satisfy  a 
known  limit.  Final Ly,  the  viscoelistic  constant,  Gj  was  probably  not  a 
constant  but  a  function  of  shear  rite  which  viries  with  stress  level.  We  also 
recognized  that  the  data  base  upon  which  all  our  conclusions  were  being  based 
was  not  the  best. 

G. 1  Re-evaluation  of  the  Upper-Newtonian  Viscosity  Limit 


As  a  result  of  further  analysis  of  the  kinetic  theory,  we  discovered 
that  there  could  be  a  more  valid  moans  of  interpreting  the  upper  Newtonian 
viscosity  than  has  leen  used  in  the  past.  In  general,  we  have  taken  the  upper 
Newtonian  viscosity  as  the  solvent  viscosity  (in  this  case  0.09  poise). 

However,  from  molec  ilar  theory  of  viscous  polymer  solutions  the  upper  Newtonian 
viscosity  is  probably  much  higher  than  this. 

From  the  analysis  by  F.  Bueche*,  the  upp.r  Newtonian  viscosity  could  be 
interpreted  as  shiwn  in  Figure  1J.  The  upper  Newtonian  viscosity  lies  at  the 
same  molecular  weig  it  as  the  lower  Newtonian  viscosity,  but  on  a  line  that  is 
the  extrapolation  o:  the  viscosity-molecular  weight.  This  line  is  marked  "A"  in 
Figure  13. 

For  PMMA  in  DE 1 : 

Me,  solution  -  (1-1»/C)  10,500  (1) 

where  C  is  the  concentration  in  gm/cm3.  At  M  , 

=-  Kj  Mc3*4  -  K2  M,.1*4  (2) 

The  values  of  Mc  and  (upper  Newtonian  viscosity)  for  the  three  solutions 
studied  by  Song  are  tabulated  in  Table  IV.  Note  how  different  the  values  are 
from  the  solvent  viscosity  of  0.09  poise. 


*  Bueche,  F.,  J.  Appl.  Phys.  26167,  738  (1955). 
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Table  IV  Molecular  Weight  Parameters  for  PMMA 

C(g/cm3)  (poise)  Mj.  K2  iv  (poise) 


0.33 

5000 

3.76 

X 

104 

4.38 

X 

10"5 

3.11 

X 

10-14 

533 

0.387 

30132 

3.20 

X 

104 

1.92 

X 

10"4 

1.87 

X 

10~13  2335 

0.443 

200000 

2.80 

X 

104 

9.73 

X 

10"4 

1.24 

X 

10-1 2 

1.18  x  10' 

In  spite  of  the  vast  difference  in  in  the  old  and  new  evaluation,  there 
was  no  difference  in  our  ability  to  fit  the  data.  At  this  point  we  simply  did 
not  know  if  the  problem  was  the  data  (Song's)  or  the  theory;  thus,  we  decided 
to  also  obtain  a  new  set  of  polymer  data. 

G.2  Re-evaluation  of  the  Theory  for  the  One-Step  Reaction 

If  one  solves  Eq.  ((J)  of  Appendix  C  of  the  proposal  for  rit,  one 

obtains 


nt  -  T/[y  -  (t/Gl) ]  (3) 

v/hich  in  the  limit  of  t  *  0  is  indeterminate  and  must  be  considered  wrong.  The 
viscosity  is  certainly  not  0,  1,  or  <»  .  If  Eq.  (9)  of  Appendix  C  Is  modified  to 
Eq.  (C)  of  Table  I  (with  G2  =  °°  ) 

t  +  (n0/c1)  t  =  nj  (4> 


as  used  by  us  earlier  and  used  by  Park  (74),  then, 

nt  *  fT  +  (n0/Gj)  x]/  y 


(5) 


which  in  the  limit  of  t  *  0  given  ■  hQ,  is  logically  correct.  We  know 

that  this  equation  cannot  fit  stress  relaxation  with  G^  constant.  It  is 
obvious  that  using  as  a  constant  was  an  error  and  that  Gj  is  actually  a 
function  of  T  or  F  (i.e.,  time  in  our  experiment). 

Two  approaches  have  been  used  for  the  systematic  evaluation  of  the 
constants  for  this  new  version  of  the  kinetic-elastic  theory.  Most  of  the 
evaluation  is  exactly  the  same  as  our  previous  endeavors  but  there  are  some 
differences  associited  with  obtaining  some  of  the  constants.  A  series  of 
programs  was  developed  in  order  t>  obtain  the  best  values  of  the  constants. 
Parallel  computatijns  were  made  it  order  to  ascertain  that  the  programs  and 
techniques  were  valid.  Both  the  main  University  computer  and  a  simulator 


package  available  on  our  departmental  VAX  li/780  system  were  used  for  this.  The 
latter  programs  are  quite  short.  The  first  program  is  used  to  adjust  the  p  and 
K.  values  from  the  equilibrium  evaluation.  This  is  exactly  as  was  done 
previously  (see  Figure  14).  The  same  program  is  further  used  to  adjust  the 
initial  slope  of  the  curve,  G^  (0),  which  gives  the  initial  values  for  the 
elasticity  parameter.  Again,  this  is  as  previously  done,  except  previously 
was  assumed  constant  (see  Figure  15).  The  values  of  p^  and  k^  are 
roughly  set  by  past  experience.  The  second  program  uses  initial  rate  theory  for 
a  more  accurate  evaluation  of  the  constant  p^.  With  this  new  value  of  pi, 
the  values  selected  in  the  first  program  are  checked.  If  they  are  inadequate, 
slight  variations  are  made.  With  the  constants  now  selected,  the  first  program 
is  used  to  establish  k^  more  accurately  so  that  the  thixotropic  stress 
intersects  the  experimental  stress  at  the  maximum  point.  This  is  a  new  fitting 
technique  (see  Figure  16).  The  same  program  is  then  used  to  see  if  k3  (taken 
as  zero)  can  help  the  shape  of  a  curve  near  equilibrium.  So  far  it  has  been 
zero.  Having  now  established  all  of  the  constants,  the  third  program  can 
calculate  G^  as  a  function  of  time  so  as  to  provide  the  final  fit.  Once  G^ 
has  been  evaluated  for  all  of  the  runs,  the  values  can  be  plotted  as  a  function 
of  shear  stress  or  F  (conversion  cf  the  thixotrophic  structure)  to  establish  the 
universality  of  the  parameter.  It  is  specifically  this  procedure,  just 
outlined,  that  will  be  described  in  more  detail  after  the  obtaining  of  n>  w 
polymer  data  is  described. 

H.  New  Polymer  Data 

H. 1  Problems 

As  previously  reported,  the  R16  unit  had  considerable  noise  that  made 
it  difficult  (impossible?)  to  establish  the  initial  slope  for  stress  growth. 

This  slope  provides  the  value  of  G^  from  the  kinetic-elastic  theory.  Most  of 
the  noise  was  removed  through  instrument  repair  and  adjustment  by  Mr.  Spooner  of 
Sangamo.  Recent  tests  on  polymethylmethacrylate  show  that  the  noise  has  been 
considerably  reduced  but  is  still  large  enough  to  introduce  more  error  than 
desired  in  obtaining  G^.  Running  the  instrument  in  reverse  and  using  the 
piezoelectric  crystal  in  tension  eliminated  the  noise  to  a  satisfactory  level. 

To  obtain  improved  data  for  testing  of  our  viscoelastic  model  for  polymer 
solutions,  it  was  necessary  to  obtain  a  complete  set  of  equilibrium  and 
transient  data  using  the  Weissenberg  Rheogoniometer.  Prior  to  the  experiment,  a 
calibration  curve  of  torque  versus  voltage  was  found  and  was  linear.  The  gap 
size  transducer  and  the  analog  recorder  were  also  calibrated.  None  of  the 
calibrations  are  presented  here  but  can  be  found  in  our  monthly  reports. 

Problems  were  encountered  with  the  range  switch  of  the  charge  amplifier;  thus, 
the  units  were  returned  to  the  manufacturer  fur  reconditioning  and  ■  alibration. 

In  spite  of  all  our  efforts  we  still  were  not  satisfied  with  the  output. 

New  data  sets  were  taken  using  the  LSI-1 1/VAX  11-780  data  acquisition  system 
(see  Figure  17)  and  as  noted  by  the  arrow  the  noise  problem  still  existed  but 
was  not  extreme.  Thus,  the  noise  was  still  coming  from  somewhere  in  the  gear 
system. 
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Experimentally  measured  stress  (TAUE) ,  calculated  Thixotropic  Stress 
(TAUT),  and  Stress  (TAU) .  Using  first  guess  values. 


Figure  15.  Experimentally  measured  stress  (TAUE) ,  calculated  Thixotropic  Stress 
(TAUT),  and  Stress  (TAU) .  Using  adjusted  G.  value. 


;ure  i  -  .  i  imentally  measure  stres:  (TAUE) ,  calculated  Thixotropic  Stress 

^  .  and  Stress  (TAT).  Usinc  adiusted  G,  and  k.  values. 


When  the  R-16  is  operated  in  the  forward  rotation  mode,  the  time  transient 
curves  are  not  smooth  as  shown  in  Figures  17  and  18.  The  cause  was  traced  to 
worm  gear  movement  and  could  be  completely  eliminated  by  operating  in  the 
reverse  rotation  mode  as  shown  in  Figures  19  and  20. 

To  test  the  effect  of  the  torsion  bar  on  noise,  the  3/8“  and  1/4”  torsion 
bars  were  used.  The  heavier  bar  reduced  the  signal  and  had  no  effect  on  the 
noise.  The  ratio  of  reduction  in  signal  is  almost  the  same  as  the  ratio  of  the 
calibration  of  k^'s  of  the  two  torsion  bars,  as  provided  by  the  manufacture. 

In  order  to  obtain  the  larger  signal,  the  1/4"  torsion  bar  was  used  for  further 
data  gathering. 

H. 2  New  Polymer  Data 


New  equilibrium  data  were  taken  at  various  shear  rates.  The  double 
logarithm  plot  of  viscosity  (poise)  and  shear  rate  (sec~^)  is  shown  in  Figure 
21.  The  behavior  of  the  curve  is  consistent  with  that  obtained  in  previous 
studies.  The  viscosity  at  zero  shear  rate  is  estimated  to  be  1.80  x  10^ 
poise. 


New  sets  of  time  dependent  shear  stress  data  for  the  38.5%  solution  of  PMMA 
in  OEP  were  obtained  using  the  Weissenberg  R-16  and  the  VAX  11/780  data 
gathering  systems.  The  data  include  shear  stress  growth  and  relaxation  at 
constant  shear  rates  over  the  range  of  shear  rates  from  0.11  sec-^  to  69.3 
sec-l.  The  shear  stress  curves  for  the  shear  rates  of  1.103  and  11.03  sec“l 
are  shown  in  Figures  22  through  25.  The  adequacy  of  these  data  should  be  clear. 
Furthermore,  no  filtering  has  been  used  on  any  of  the  data. 


1.  Detailed  Comparisons  of  Modified  Theory  and  New  Polymer  Data 


In  order  to  present  the  results  in  their  proper  context,  we  will  review  the 
theory  as  it  is  now  being  used.  The  kinetic-elastic  model  developed  in  this 
laboratory  involves  two  facets.  One  is  the  kinetic  model  for  thixotropic  change 
of  the  viscoelastic  material,  where  the  rate  of  change  of  a  material  structure 
is  described  by  a  rate  concept.  The  other  is  the  elastic  model  and  can  be 
des -ribed  by  a  modified  version  of  Oldroyd's  viscoelastic  equation. 


The  thixotropic  kinetic  theory  can  be  expressed  by  the  following  equation: 


A 


2B 


(6) 


Utilizing  the  homogeneous  kinetic  reaction  concept,  a  simple  kinetic  equation 
can  be  obtained  (as  shown  in  Appendix  C  of  the  proposal) 


d(CF) 


dt 


k^(CF)1  -  kJ[C(l  -  F)]2 


(7) 


I  f 

where  F  is  the  fraction  of  the  thixotropic  fluid  structure  unchanged,  k^,  k2 
are  forward  and  reverse  reaction  constants,  and  C  is  the  concentration  of  the 
polymer  solution. 
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Figure  18.  Stress  relaxation  with  dip  in  relaxation  curve  (after  forward  rotation). 


VitlV! 


Stress  growth  using  reverse  rotation. 
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Stress  relaxation  after  reverse  rotation  using  1/4"  torsion  bar. 
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Figure  25.  Shear  Stress  Relaxation  (11.03  sec 


The  function  of  the  fraction  of  the  thixotropic  fluid,  F,  is  further 
defined  as  _ 

nt  "  n» 

F  -  -E - =  , 

na  -  na  ( 

O  00 


the  viscosity  term;  nt»  is  the  thixotropic  structural  viscosity.  r\Q  and  hc» 
are  the  viscosities  at  zero  shear  rate  and  infinite  shear  rate,  respectively.  H0 
and  Hoo were  79,950  poise  and  4000  poise,  respectively,  for  this  specific  polymer 
solution.  The  constant  a  was  selected  to  be  (1/3.4)  from  the  evidence  (Fox  & 
Allen)  that  the  zero  shear  rate  viscosity , 0 Q,  is  proportional  to  the  3.4 
power  of  the  molecular  weight. 

The  rate  constants  k^  and  k2  in  Equation  (7)  are  functions  of  shear 
stress,  T  ,  and  temperature.  The  functional  forms  can  be  expressed  by  Equation 
(9)  and  (10)  as 


k3  +  k2T 


For  simplification,  kj  in  Equation  (10)  is  taken  as  ze  o  in  this  model 
analysis.  At  equilibrium  (or  steady  state),  the  forward  ra  e  and  backward  rate 
are  equaL,  i.e.,  d  ;CF)/dt  *  0  from  Equation  (7).  The  equilibrium  constant  can 
b‘  •  xpressed  in  the  following  form: 


Cd  -  F  ) 


Pi  -  P2 


K  TP 
eq  eq 


where  Feq  and  T  eq  are  the  F  and  t  at  equilibrium  conditions.  Based  on  the 
equilibrium  experimental  data,  and  P  were  found  to  be  1.313  x  10“^  and 

1.864,  respectively.  The  values  of  kj  and  in  Equation  (9)  can  be 

calculated  from  the  initial  rate  change,  i.e.,  the  rate  at  zero  time,  p^  was 
found  to  be  6  x  10-^  for  all  of  the  experimental  data. 

The  modified  Oldroyd  viscoelastic  equation  can  be  expressed  by 

T  +  (n  /G.)  X  »  o  Y 
01  t 


as  was  shown  earlier  in  section  G  2  and  where  x  and  t  are  the  shear  stress  and 
its  first  derivative,  Y  is  the  shear  rate,  and  G^  is  the  elastic  modulus 
proposed  to  be  a  function  of  T  or  F  (i.e.,  time  in  our  experimental)  for  our 
kinetic-elastic  model. 


Figure  26  shows  the  calculated  Gj^  values  as  a  function  of  time  from  one 
set  of  experimental  data.  It  may  be  that  Gj  is  a  simple  function  of  F,  as  is 
Ht  [Eq.  8],  Note  the  similarity  of  G^  versus  t  to  such  plots  as  H  or  G' 
versus  Y  .  The  G^  function  can  be  expressed  mathematically  by  Equation  (13) 

G^t)  =  (Ast)  +  G1(0)  exp  [-G^TC)  (t  -  t  )]  (13) 

where  G^(Ast)  is  the  asymptometrie  value  of  G^  function, 

G^CO)  is  the  G^  value  at  zero  time, 

tout  cfme  which  Gj  starts  to  decay,  and 

G^(TC)  is  the  decay  constant  started  at  tCUf 

The  best  fitting  kinetic  and  elastic  constants  for  the  kinetic-elastic  model  at 
various  shear  rates  are  listed  in  Table  V.  Examples  of  the  calculated  curves 
and  experimental  curves  for  the  individual  shear  rates  are  shown  in  Figures  27 
through  29.  Each  curve  in  this  initial  fit  was  optimized  individually  so  that 
t  he  various  constants  would  be  representative  and  any  dependency  on  shear  rate 
level  could  be  determined. 

As  can  be  see  i  from  the  examoles  given,  the  comparisons  are  excellent. 
However,  some  of  tie  constants  ta  lulatecl  in  Table  V  do  vary  with  the  shear  rat'. 
Hotl  G|(o)  and  Gj(4st)  are  functions  of  (log  Y  )  as  expectel,  with  Gj(0) 
increasing  and  G[(  \st)  decreasing  with  an  increase  in  Y.  i'lt  •  kj  value 
increases  linearly  with  Y«  P  anl  Fj  are  constant.  The  ot  ver  terns  [K.eq, 
t cut  >  G] ( TC) ]  vary,  but  do  not  co ‘relate  with  y  . 

J.  Implications  fir  Further  Work  on  Polymer  Systems 

There  is  an  alternate  evaluation  that  should  be  made.  This  is  to  use  those 
terms  that  did  not  correlate  with  shear  rate  [Keq,  tcut,  Gj(TC)]  as 
constants  at  their  average  values  and  re-evaluate  the  terms  that  are  a  function 
of  shear  rate.  This  should  give  a  clearer  and  more  consistent  correlation  for 
(^(0),  Gj(Ast),  and  k^. 

The  second  investigation  involves  a  modification  of  the  theory  to  ascertain 
whether  or  not  using  keq,  tcut,  and  Gj(TC)  as  constants  will  have  any 
effect.  The  modification  is  associated  with  a  change  in  the  viscoelastic  part 
of  the  theory.  Currently,  although  G^  is  computed  as  a  function  of  time,  it 
lias  been  taken  outside  of  the  difierential  and  computed  as  a  difference.  This 
can  be  illustrated  best  by  observing  the  equation  that  was  used 


Func 


and  the  following  equation,  which  would  be  more  accurate 


X  +  *  ntY 


that  is,  Gj  is  now  varying  with  time  and  should  be  inside  the  differential: 

(x?Gx)  =  d(T/Gi)/dt  (15) 

At  this  stage  it  may  well  be  that  the  older  formulation  which  seems 
logically  more  correct  could  be  utilized,  i.e., 

T  +  r^Cx'/Gp  =  ntY  (16) 

where  rit  *-s  used  in  place  of  rt0.  This  substitution  could  give  more  reasonable 
variations  of  k^  and  Gj.  In  addition,  it  could  aid  in  fitting  stress 
relaxation  data. 

These  equations  [(14)  and  (16)]  can  be  tested  by  our  simulation  with  minor 
changes,  although  the  use  of  eq.  (16)  will  introduce  a  trial  and  error  iteration 
for  the  Gj  time  function. 

These  formulations  might  improve  the  logical  nature  of  the  variation  of  the 
elastic  constants  in  the  forward  reaction  rate;  however,  they  will  not  improve 
the  already  excellent  fit  that  has  been  obtained  by  allowing  G^  just  to  be  a 
variable.  Consequently,  this  further  work  can  not  only  help  us  zero  in  on  the 
most  logical  theory,  but  could  help  in  improving  the  fit  for  stress  relaxation. 
Finally,  we  would  like  to  extend  the  best  analysis  to  the  treatment  of 
oscillatory  data. 


II. 


SLURRY  MATERIALS 


All  of  the  proposed  work  on  the  preliminary  slurry  formulations 
supplied  to  us  by  Sun  Tech  and  Exxon  (via  Wright-Patterson  AFB)  have  been 
completed.  In  future  work  we  would  repeat  these  tests  for  the  selected 
final  formulations  and  extend  the  measurements  to  oscillatory 
conditions.  We  also  want  to  further  investigate  the  specific  reasons 
for  the  non-Uewtonian  behavior  by  using  combined  visual  and  rheological 
unsteady  state  measurements.  Such  results  should  allow  us  to  improve 
the  logical  formulation  of  the  kinetic-elastic  theory  for  the 
description  of  these  materials. 


A.  Equipment  and  Procedures 


A.  1  Jarbon  Black  Samples 

T  ie  samples  used  in  this  study  were  suspensions  of  microsize  carbon 
black  nartic.es  dispersed  in  JP-10  solvent.  JP-10  is  a  pure  hydrocarbon 
winch  is  used  as  a  jet  fuel.  It  has  higher  combustion  energy  than  the 
curren;  commercial  jet  fuel,  JP-4,  but  it  is  more  expensive  due  to  its 
purity.  CarDon  black  is  added  to  JP-ll)  to  increase  its  density,  and 
hence,  tne  a  .  ount  of  energy  available  per  volune. 


Four  different  carbon  black  suspensions  were  used  in  this  study: 
Exxon  708-61,  Sun  Tech809-995,  Sun  Tech339-919,  and  Sun  Tech839-988. 

The  basic  physical  and  chemical  properties  of  JP-10  and  the  carbon  black 
slurries  are  shown  in  Table  VI.  As  can  be  seen,  the  carbon  black 
loading  is  about  50  weight  %,  or  0.33  volume  fraction  based  on  1.86  and 
0.93  specific  gravities  for  carbon  black  and  JP-10  respectively.  The 
particle  size  is  approximately  0.1  m.  The  small  particle  size  and  high 
loading  level  generally  result  in  a  yield  stress.  Since  we  were 
interested  i:  materials  with  a  yield  stress,  more  emphasis  was  put  on 
the  yield  materials:  Sun  Tet  h  839-988  a  id  Exxon  708-61. 

Due  to  the  high  vapor  pressure  of  the  solvent,  evaporation  occurs 
causing  the  slurry  to  sxin  over  at  the  edge  of  the  plate  and  cone.  The 
change  in  solvent  content  not  only  changed  the  measured  viscosity,  see 
Figure  30,  but  also  made  the  measurement  difficult.  A  Mooney  type 
system,  which  combines  a  cylinder  with  l  cone  and  plate,  was  used  to  try 
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Table  VI 

Physical  and  Chemical  Properties  of  JP-1 0 


Property 

JP-10 

Formula 

c10h16 

Molecular  Weight 

136 

Heat  of  Combustion 

142, 115 

(BTU/gallon) 

Flash  Point,  F 

130 

Freezing  Point,  F 

-1  30 

Boiling  Point,  F 

360 

Auto  ignition  Temp.,  F 

4  74 

m 


Physical  and  Chemical  Properties  of  Slurries 


Sample  Nunber 
Composition 

Carbon  Black  type 
Carbon  Black  size,  ns 
Carbon,  Wt .  Z 
Dispersant  Type 
Dispersant,  Wt .  Z 
JP-10,  Wt.  Z 
Heat  of  Combustion,  Calc 
(BTU/gallon) 
Density  gm/ral  @20  C 
Particle  Size,  p  by  FOG 
Average 
Largest 

Yield  Stress  of  25  C 


708-61 

839-988 

United  SL90 

SAF 

103 

70-80 

51 

50.4 

W 

ST-107 

2 

2.0 

47 

47.6 

168,300 

173,  700 

1.24 

1.234 

<6 

13 

10 

— 

220.0 
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Figure  30*  Viscosity  Change  Due  to  Skinning  During  Measurement 


to  solve  the  skinning  problem.  A  guard  ring  was  placed  on  the  top  of 
tiie  cylinder  to  minimize  the  surface  area  that  is  exposed  to  air.  Thus, 
it  should  have  minimized  the  vaporization  of  solvent  from  the  slurry. 
Unfortunately,  the  setup  did  not  stop  the  solvent,  JP-10,  from 
evaporating.  There  was  a  small  space  between  the  guard  ring  and  the 
outer  cylinder,  so  that  the  skinning  occurred  in  this  space.  This  not 
only  caused  difficulty  in  the  outer  cylinder's  rotation,  but  also 
affected  the  motion  of  the  torsion  head.  This  approach  was  abandoned 
and  new  methods  were  investigated.  Evaporation  could  be  avoided  by 
keeping  the  sample  in  a  saturated  vapor  pressure  environment.  Since 
this  procedure  would  be  difficult  to  use  with  the  R-16,  another  approach 
was  tried.  Attenpts  were  made  to  cover  the  sample  with  non-miscible 
solvents  to  prevent  the  JP-10  from  evaporating.  Water  was  tried  first, 
out  due  to  gravity,  the  water  mixed  with  the  slurry.  Next,  pure  JP-10 
was  used  for  th<.  same  purpose. 

The  slurry  sample  was  [oured  until  the  cylinder  was  filled  to  the 
free  surface  at  the  edge  of  the  plate.  Then,  JP-10  was  carefully  put  on 
the  top  of  t.ie  sample.  In  this  way,  evaporation  was  confined  to  the 
JP-ivJ,  which  was  not  under  shear.  The  excess  slurry  sam  >le  and  JP-1J 
nad  negligible  effect  on  the  measured  torque  which  was  tested  by  using 
the  above  setup  with  Newtonian  fluids  of  3.5  Poise  and  9.85  Poise  as  the 
primary  samples.  The  error  introduced  was  less  than  4%,  which  is  within 
the  experimental  error.  Thus,  in  all  runs  JP-10  was  carefully  placed  on 
top  of  the  slurry  as  a  coating  to  prevent  the  skinning  problem.  The 
mixing  between  JP-10  and  the  sample  took  place  by  molecular  diffusion 
only.  Mass  transfer  by  diffusion  was  negligible  during  the  measurement. 
This  was  checked  by  measuring  the  equilibriun  viscosity  for  samples 
resting  between  u.5  to  6  hours.  As  shown  in  Table  Vll,  the  equilibrium 
viscosity  at  a  shear  rate  of  27.6  sec~l  was  independent  of  the  resting 
time  with  tae  JP-10  coating  on  the  sample. 

At  low  temperatures,  condensation  of  water  on  the  sample  was  a  more 
serious  problem  than  skinning.  Especially  when  the  temperature  Is  lower 
than  the  freezing  point  of  water,  0  C,  it  is  important  to  isolate  the 
sample  from  moisture.  During  the  restin;  period,  a  plastic  cover  was 
used  to  cover  the  sample,  which  prevented  both  condensation  of  water  and 
vaporization  of  the  sample.  Inert  nitrogen  gas  was  slowly  bled  into  the 
chamber  to  Keep  the  plate  and  cone  in  an  inert  environment.  This  worked 
quite  satisfactorily,  except  at  the  top  of  the  chamber  which  was  open  to 
the  atomsphere.  However,  tills  condensation  did  not  interfere  with  the 
measurements. 


A.  2  Rheometer 


Auxiliar ies 


In  this  sti._  .  ail  the  data  were  obtained  on  a  Weissenberg 
Rheogoniometc  • ,  R-16,  and  its  supporting  electronic  devices.  The 
measuring  sys  ;em  can  be  divided  into  three  subsystems:  the  rheometer 
itself,  signal  conditioning  equipment,  and  data  processing  facilities  , 
as  shown  in  Figure  31. 


A. 2.1  Weissenberg  Rheogoniometer 


The  Weissenberg  Rueogoniometer ,  R-16,  is  a  plate  and  cone  rheometer 
which  provide.,  a  constant  shear  rate  in  the  gap.  The  R-16  has  been 
modified  and  updated.  The  current  version  has  the  capability  of 
measuring  the  shear  viscosity,  normal  stresses,  dynamic  viscosity,  and 
the  more  comp  ex  situation  ol  superposi, ion  of  rotation  and  oscillation. 
In  tnis  sectii n,  a  brief  description  of  the  R-16  is  given,  with  an 
emphasis  on  t,;e  newer  facilites.  More  stalled  descriptions  of  the  R-16 
have  uecn  pro\  idea  by  Denny(.l),  Lee(82  ,  and  Song(ll)  ani  the  operation 
manual (83).  he  rheometer  hi  s  three  ma  ir  parts:  the  drive  unit,  thi 
plate  and  con.  section,  and  the  torsion  lead . 

Tlie  driv.-  units  consist  of  two  mot.  rs  and  two  speed  controls.  A 
three  phase  synchronous  motor  and  a  geat  box  are  used  to  control  the 
rotational  speed  of  the  bottom  plate  at  60  selected  output  speeds.  A  DC 
motor,  fixed  , ear  box  reducer,  and  SCR  control  are  used  to  set  the 
frequency  of  oscillation.  The  instrunent  allows  torsional  testing  in 
oscillation  a  well  as  in  steady  rotation,  and  is  capable  of  measuring 
not  only  viscosity  but  also  elasticity,  Jilatancy,  thixotropy,  etc.  In 
addition,  by  i.sing  a  magnetic  clutch  to  start  or  stop  rotation, 
measurements  ;s  a  function  of  time  (transient  response)  can  be  obtained. 
The  response  time  of  the  instrunent  is  snail,  less  than  0.01  sec(75). 


The  plate  and  cone  are  the  heart  ol 
is  placed  between  them.  The  bottom  plat 
units  in  steady  rotation  and/ >r  in  a  sin 
mom..  .at  transferred  fron  the  bottom  pin 
sample  is  meat ured  by  the  torsion  head, 
the  snear  stress,  and  thus  viscosity 
can  cover  rates  of  shear  between  7x10 
changing  the  rotational  speeds  and  angle 
7.5  cm  plate  n  id  a  58  minute  4b  second  ( 


,-4 


the  measurement,  and  the  sample 
e  can  be  moved  by  the  drive 
e-wave  oscillation.  The 
:e  to  the  cone  through  the 
Much  can  be  used  to  determine 
ie  Weissenberg  Rheogoniometer 
Sec'1  to  9x10^  sec-1  by 
s  of  the  cones.  In  this  study  a 
luminal  1°  )  cone  were  used. 
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Figure  31.  !<—  1  (j  Rlieome lu r  and  auxiliaries 


RECORDER 


The  shear  rate  range  varied  from  3.54x10-3  to  10^  sec"-*-. 

The  torsion  head  includes  an  air  bearing,  torsion  bar,  and 
transducer.  The  air  bearing  is  used  to  provide  frictionless  rotation 
for  torque  measurement.  The  torsion  bar  and  transducer  are  used  to 
determine  the  torque  and  thus  the  shear  stress.  Torsion  bars  with 
diameters  of  1/16"  and  1/8"  were  used  depending  on  the  magnitude  of 
torque . 

Because  of  the  improved  stiffness  of  the  system,  the  piezoelectric 
crystal  load  cell  has  proven  to  be  a  better  force  sensor  than  the  LVDT 
(linear  variable  differential  transformer),  as  shown  by  Brodkey 
et  al.(6),  and  by  Crawley  and  Graessley(84) .  However,  the  torque 
generated  by  the  slurries  was  too  small  to  be  detected  by  crystals. 
Instead,  a  LVDT,  type  F51TM  made  by  Boulton  Paul  Aircraft  Ltd.,  was  used 
to  measure  the  torque.  It  can  produce  an  electrical  output  proportional 
to  the  displacement  of  a  separate  movable  core.  This  transducer  has  a 
high  sensitivity  for  the  measurement  ot  small  displacements  up  to 
O.lxlO-^  inch  with  reasonable  accuracy. 


A. 2 . 2  Signal  Conditioning  Equipment 

The  signal  conditioning  subsystem  consists  of  a  transducer 
conditioning  meter  (type  EP597,  B.P.A.  Electronics)  and  an  analog 
filter.  The  output  signal  from  the  transducer  is  fed  to  the 
conditioning  meter.  The  proper  range  of  the  meter  is  selected  to 
amplify  and  transform  the  signal  so  as  to  have  an  output  voltage  within 
^10V  to  meet  the  requirements  of  the  A/D  converter  for  data  acquisition. 
The  output  voltage  from  the  meter  then  passes  through  a  filter  to  remove 
the  carrier  frequency.  A  simple  RC  filter  with  a  cutoff  frequency  of 
2000Hz  is  used. 
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A. 2. 3  Data  Processing  Facilities 


The  data  processing  system  includes  a  microprocessor  (Digital 
Equipment  Co.,  LSI  11/02  with  an  RTP-A/D  converter),  a  VT100  graphics 
terminal(Digital  Engineering),  chart  recorder(Of fner  Electronics  Inc.), 
and  minicomputer  (Digital  Equipment  Co.  VAX  11/780)  with  all  of  its 
supporting  devices.  The  new  data  acquisition  system  was  a  major  change 
from  previous  equipment  (6,75),  and  has  greatly  improved  the  data 
acquisition  and  data  processing.  With  the  on-line  minicomputer,  data 
can  be  taken  fast  and  accurately.  Data  analysis  also  becomes  easier 
because  of  the  ready  access  to  a  fast  computer  with  a  high  level  of 
interactive  graphics.  The  present  system  replaces  the  PDP-15  dual 
processor  installed  about  10  years  ago. 

The  microprocessor  has  32  K  memory,  of  which  24  K  was  used  for  data 
storage,  and  8  K  was  used  for  the  data  acquisition  program.  The 
microprocessor  contains  a  12-bit  analog-to-digital  (A/D)  converter  which 
accepts  signals  between  -10  and  +10  V  and  converts  them  to  decimal 
values  between  -2000  and  2000  counts  (or  -2047  to  2047).  The  A/D 
converter  is  a  zero  order  hold  and  is  able  to  digitize  eight  signals 
simultaneously.  A  crystal  oscillator  generates  the  pulse  which 
subsequently  samples  and  acquires  data.  The  LSI  11/02  can  sample 
signals  up  to  300K  Hz.  The  sampling  frequency  is  determined 
interactively  by  a  table  which  can  have  at  most  20  elements.  These 
elements  are  the  sequence .rate , and  amount  of  data  to  be  acquired.  In 
general,  the  LSI  is  a  dummy  device  without  an  operating  system.  It  has 
to  be  initialized  by  the  VAX  11/780  minicomputer  first  by  downloading  an 
absolute  loader  and  the  data  acquisition  program  (object  code  of  LSIPGM 
program).  It  can  only  start  to  acquire  data  after  the  downloading  is 
completed.  The  microcomputer  takes  the  data  at  a  rate  based  on  the 
frequencies  set  by  the  operator,  stores  the  data  in  Its  own  memory,  and 
later  transfers  the  data  to  the  VAX,  wh.ch  stores  it  on  disk. 

The  VAX  11/780  is  a  minicomputer  with  a  semiconductor  memory  of  1.5 
megabytes  and  a  virtual  memory  of  200  megabytes.  All  computations  were 
done  with  the  aid  of  this  unit.  Programs  were  stored  on  the  VAX  disk 
and  were  then  used  for  computation  or  downloading  to  the  LSI  for  data 
acquisition  through  a  VT100  terminal.  "he  VT100  graphics  terminal  was 
used  to  communicate  with  the  rainicomput  r,  preview  data,  and  evaluate 
the  results. 


A  Dynograph  recorder  was  used  to  simultaneously  record  the  signal 
as  it  was  being  digitized.  The  recorder  contains  type  461  preamplifiers 
and  a  type  462  dual  channel  amplifier.  It  was  used  to  monitor  the 
tranducer  output  to  help  decide  if  the  data  was  satisfactory  for 
transfering  to  the  VAX.  A  digital  voltmeter  (Fluke  8020A  Multimeter, 
John  Fluke  MFG.  Co.,  Inc.)  was  also  used  to  monitor  the  magnitude  of 
signal  to  ensure  that  the  signal  was  in  the  range  of  ±10V. 


A. 2. 4  Thermostat 


A  cryogenic  unit,  Ultra  Kryostat  UK-60-SD(made  in  Germany),  was 
used  to  control  the  temperature.  Freon  R-13-B1  which  could  cool  to 
-60  C  in  an  hour,  was  used  in  this  uni:.  Methanol  was  used  as  the 
cooling  medium  for  the  rheogoniometer  and  was  circulated  through  the 
chamber  which  surrounds  the  plate  and  cone.  The  temperature  and  the 
flow  rate  of  methanol  could  be  set  at  the  cryogenic  unit.  Water-free, 
pure  nitrogen  was  slowly  blown  into  the  chamber  to  enhance  the 
convective  heat  exchange  and  prevent  it,  >isture  condensation.  The 
temperature  was  measured  on  the  cone  b/  means  of  a  digital  thermister 
thermometer  (Keithley  870  Digital  Thet  lometer).  In  a  typical 
experiment,  the  temperature  remained  c  mstant  to  within  +0.5  C. 


A. 3  Calibrati m 

Although  the  Weissenberg  Rheogoniometer  is  widely  used  to  study  the 
viscoelastic  properties  of  materials,  0avis(87)  recommended  that  the 
user  calibrate  the  instrument  with  Newtonian  fluids  of  known  viscosity. 
He  found  a  30%  difference  between  the  experimental  and  manufacturer's 
values  of  the  torsion  bar  constants.  He  also  found  that  the  constant 
used  to  convert  the  ratio  of  transducer  meter  output  voltages  to 
amplitude  ratio  was  15%  smaller  than  quoted  in  the  manual. 

In  order  to  obtain  reliable  data  from  our  Weissenberg  R-16,  each 
piece  of  equipment  was  checked  and  calibrated  carefully.  The  major 
components  included  the  transducers,  transducer  meters,  recorder,  and 
computers . 
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A. 3.1  Drive  Unit 


The  rotational  speed  of  the  bottom  plate  is  set  by  the  gear  box 
setting.  A  stroboscope  was  used  to  try  to  measure  the  rotational  speed 
of  the  plate.  Because  the  rotational  speed  of  the  plate  was  too  small 
to  be  detected  by  the  stroboscope,  a  stopwatch  and  optical  observations 
were  used.  Two  speeds  were  selected  for  comparison  with  the  values 
given  in  the  manual.  The  results  are  shown  in  Table  VIII. 


A. 3. 2  Transducers  and  Transducer  Meters 

There  are  three  Type  EP597  transducer  meters(TM).  Each  transducer 
meter  was  checke  1  with  the  same  transducer  (B.P.A.  Electronics,  Type 
F13B,  series  N  4*0)  which  is  the  type  used  for  gap  setting.  The 
amplification  ratio  of  the  first  meter  was  1.31  times  larger  than  that 
of  the  second.  7or  the  100  range, 

V  =  22.3  d  for  TM  it  t 

V  -  17. s  d  for  TM  it  > 

where  V  is  voltage  output  and  d  is  the  transducer  displacement. 

Since  it  was  more  convenient  when  setting  the  gap  to  use  the  meter 
reading  as  a  displacement  measurement,  the  calibration  was  done  through 
the  meter  reading.  Most  gap  sizes  (depending  on  the  cone  angles)  fall 
in  the  meter  ranges  of  2.5  and  10,  i.e.,  between  2.5  and  10.0 
milliinches.  Therefore,  the  range  10  was  tested  first  to  establish  the 
relationship  between  meter  reading  and  actual  displacement.  The  meters 
were  zero  set  first.  Then,  on  range  10,  the  gain  of  the  gap  setting 
transducer  (type  F13,  Series  No.  435)  was  adjusted  so  that  the  meter 
would  have  a  full-scale  deflection  when  r.he  input  displacement  was  0.01 
inch.  Unfortunately,  no  matter  what  the  gain,  the  second  meter  did  not 
have  full-scale  deflection  on  range  10  when  the  input  displacement  was 
0.01  inch.  There  was  a  linear  relationship  between  deflection  and 
displacement.  Thus,  the  second  meter  wa  -  selected  for  the  torque 
measurement  while  the  first  meter  was  us  ?d  for  the  gap  setting. 
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All  Che  transducers  were  calibrated  using  a  micrometer  to  establish 
the  constant  k,  where  V=k  X  d.  V  is  the  voltage  output  from  the 
transducer  meter,  and  d  is  the  displacement  of  the  transducer  in 
milliinches.  The  gap  setting  transducer  was  calibrated  by  mounting  the 
transducer  in  a  micrometer  calibration  unit.  The  transducer  meter  was 
zero  set  first.  Next,  the  displacement  of  the  micrometer  (in  both 
directions)  and  the  meter  reading  or  voltage  output  were  recorded 
simultaneously.  The  calibration  curve  of  meter  at  range  10  is  shown  in 
Figure  32  and  the  data  are  shown  in  Table  IX.  The  torsion  head 
transducer  differs  from  the  gap  setting  transducer  in  that  its  stator 
and  armature  are  not  connected  mechanically  and  thus  do  not  contribute 
any  friction  during  measurement.  Since  this  tranducer  cannot  be 
calibrated  in  the  micrometer  unit,  the  micrometer  was  mounted  directly 
on  the  rheometer  to  measure  the  displacement.  Both  1/3"  and  1/16" 
torsion  bars  were  used  to  calibrate  for  different  meter  ranges  (1.0-25). 
It  was  found  that  the  correlation  consta  its  are  independent  of  the 
torsion  bar  as  expected(see  Figure  33  an  I  Table  X).  The  percent  error 
of  calibration  is  3%  or  less. 

For  data  icquisition,  the  LSI  11/02  microprocessor  was  used  to 
record  the  data  rather  than  the  transduc-r  meter.  One  calibration  (on 
’EM  // 1 ,  range  2.5)  was  done  by  using  the  transducer  meter  and  the 
computer,  simultaneously.  The  results  c necked  well  as  shown  in  Table 
XI.  The  A/D  converter  provided  200  counts/volt.  Tests  on  the  second 
meter  showed  the  same  results.  In  Table  XI,  D  is  the  reading  taken  from 
the  micrometer,  V  is  output  voltage  ,  and  C  is  computer  output. 

The  entire  system,  plate  and  cone,  transducer,  and  transducer 
meter,  was  calibrated  by  standard  Newtonian  Brookfield  solutions  of 
9.85P  and  118. A  P.  Several  torsion  bars  (1/16",  1/8",  and  1/A")  were 
used  to  measure  the  viscosities,  as  shown  in  Figure  3A.  A  118. A  P 
Brookfield  solution  and  a  9.85P  Brookfield  solution  were  tested  using  a 
1/16"  torsion  bar  as  shown  in  Figure  35.  The  others  are  similiar.  The 
results  are  summarized  in  Table  XII.  The  average  error  was  less  than 
5%,  which  means  that  the  torsion  bar  constants  provided  are  reasonably 
accurate. 

In  contradiction  to  the  Davis(87)  observation,  the  torsion  bar 
constants  of  our  rheometer  agreed  to  within  5%  of  the  values  stated  by 
the  supplier.  The  same  results  were  found  by  Enthoven  and  Jalink(86). 


Typical  Displacement 


Table  IX  Calibration  of  Transducer  Meter 
at  Rang  e  10 


displacement 
*10**3  inch 

meter  reading 
(full  scale  100) 

10 

100.0 

8 

79.8 

6 

60.0 

4 

40.2 

2 

20.0 

0 

0.1 

-2 

-20.4 

-4 

-40.0 

—6 

-60.2 

-8 

-80.4 

-10 

-100.0 

Table  X  Calibration  Constants  k  from  V«k*d 


Calibration 

constants  k 

RANGE 

1/8" 

1/16" 

A  VS  k 

1.0 

1.5663 

1.6241 

1.5952 

2.5 

0.6531 

0.6570 

0.6551 

10. 0 

0.1617 

0.1642 

0. 1630 

25.0 

0.06562 

0.06654 

0.0661 

33.  Calibration  for  Two  Torsion  Bars. 


Brookfield  Solution  118.4  Poise 
Files:  ANEW4.  ANEW3,  ANEWB 


Three  Torsion  Bars  Used  to  Measure  Standard 
Solution  Viscosity  of  118. A  Poise. 


During  a  low  temperature  measurement,  the  armature  o£  the  torsion 
head  transducer  was  accidentally  broken.  The  armature  was  repaired  and 
the  transducer  recalibrated.  The  2.5  and  10.0  ranges,  with  the  1/8" 
torsion  bar,  were  checked.  The  results  were  close  to  the  previous 
values  (Table  XIII),  so  the  original  calibration  constants  were  used. 

On  April  20,1982,  TM#2  started  to  show  nonlinearity  in  its  output; 
it  was  replaced  by  TM#3.  The  gap  setting  transducer  was  used  first  to 
compare  TM#1  and  TM#3,  and  it  was  found  that  the  gains  of  the  two 
transducer  meters  were  almost  the  same.  The  constants  are  shown  in 
Table  XIV. 

TM#3  was  then  used  for  gap  setting  and  TM//1  was  used  for  the 
torsion  head  measurements.  TM#1  was  again  calibrated  with  the  torsion 
head  transducer.  The  calibration  constants  are  shown  in  Table  XIV. 

A  quick  test  was  made  to  measure  the  viscosity  of  water  at  25  C. 
The  viscosity  of  water  was  0.9671  cp  which  is  close  to  the  literature 
value  of  0.95  cp.  The  measurement  of  the  water  Is  in  the  most  critical 
range  for  the  instrument  since  high  shear  rates  and  the  most  sensitive 
meter  range  must  be  used.  The  result  was  better  than  expected. 

The  transducer  meters  are  quite  old  (20  plus  years)  and  cannot 
stand  working  too  long.  The  transducer  meter  #1  failed  during  a  low 
temperature  measurement  and  it  was  found  that  some  transistors  wet  out 
of  range  in  both  TMiPl  and  TM#2.  They  have  been  repaired  and 
recalibrated,  but  should  be  replaced  when  funds  become  available.  The 
new  calibration  constants  are  given  in  Table  XV. 


A. 4  Noise  Analysis  and  Filter  Technl  }ue 

Hie  noise  of  the  system  was  analysed  and  the  sources  of  the  noise 
were  identified.  The  output  of  the  transducer  meter  was  found  to 
contain  an  appreciable  amount  of  carrier  frequency  "noise".  The  output 
of  the  transducer  meter  is  d.c.,  but  there  is  a  high  frequency  carrier 
ripple  in  the  output  which  is  significant  on  the  more  sensitive  ranges. 
The  amplitude  of  the  noise  depends  on  the  d.c.  signal.  The  other 
predominant  noise  is  from  mechanical  sources  and  is  due  mainly  to  the 
natural  frequency  of  the  torsion  bars.  The  natural  frequencies  of  the 
bars  are  low  (around  10Hz  for  1/16"  and  60  Hz  for  1/8")  which  prevent 
using  analog  filters  to  remove  the  noise.  Filtering  of  the  noise  was 


Table  XIII 

Conver  sion 

Constants 

Range 

Original 

New 

Constants 

Constants 

2.5 

0.6551 

0.6424 

10 

0.1630 

0.1597 

Table  XIV 

Torsion  Head 

and  TMfJl 

>  8.7464 

8.7083 

2.1687 

2.1414 

0.8758 

0.8590 

0.2169 

0.2152 

0.0868 

0.0870 

0.0220 

0.0216 

Table  XV. 

Calibration  Constants 

k 

7.8547 

1.9396 

0.78066 

0.19443 

0.07821 

0.01963 


done  by  using  a  digital  filter  technique  (fast  Fourier  transform). 

The  carrier  noise  from  the  transducer  meters  was  measured  and 
filtered.  The  carrier  noise  as  estimated  by  using  an  oscilloscope  was 
2000  Hz.  Since  the  maximun  sampling  frequency  used  during  data 
acquisition  was  1000  Hz.  a  low-pass  RC  filter  with  a  cut-off  frequency 
of  2000Hz  was  used  at  the  output  of  the  transducer  meter  to  eliminate 
the  carrier  frequency.  At  zero  output,  the  output  with  and  without  the 
filter  is  shown  in  Figure  36.  The  filter  reduced  the  signal  by  only 
4.5%. 


Further  tests  were  conducted  to  identify  other  sources  of  noise. 
First,  noise  was  measured  with  the  motor  and  air  bearing  air  supply  shut 
off.  Noise  of  23  Hz  and  36  Hz  was  present,  which  must  have  come  from 
the  transducer  meter  (see  Figure  37).  Next,  an  amplifier  was  connected 
to  the  output  of  the  meter  to  see  if  it  would  contribute  noise.  The 
result  showed  that  no  additional  noise  was  coming  from  the  amplifier. 
Then  the  air  was  turned  on  and  the  torsion  bar  could  freely  vibrate  at 
its  natural  frequency.  The  noise  from  the  torsion  bar's  natural 
frequency  was  the  major  source  of  noise,  as  shown  in  Figure  38. 

Finally,  different  sampling  frequencies,  100  Hz,  500Hz,  1000  Hz,  were 
used  to  see  if  there  was  any  difference  in  the  frequency  domain 
analysis.  The  results  showed  no  difference  in  the  noise  frequency 
analysis, but  the  faster  the  sampling  period,  the  more  accurate  was  the 
frequency  analysis. 

A  filtering  technique  which  couples  curve  fitting  and  fast  Fourier 
transform  (FFT)  was  used  to  analyze  and  to  filter  the  noise.  First,  a 
curve  fitting  technique  was  applied  to  the  raw  data.  The  deviatory 
data,  the  difference  between  the  raw  da  a  and  fitted  value,  was  then 
analyzed.  The  deviatory  data  was  transformed  from  the  time  domain  to 
the  frequency  domain  by  a  fast  Fourier  transform  (FFT).  The  frequencies 
were  then  analyzed  to  determine  which  frequencies  were  noise. 

The  noise  was  filtered  out  by  a  suitable  filter:  low-pass, 
high-pass,  band-pass, or  band-reject  filter.  Next,  the  filtered 
deviatory  data  was  readded  to  the  fitted  curve  to  reconstruct  the 
filtered  data.  In  this  case,  the  filtered  data  retained  the  original 
trend  and  minimized  the  "leakage"  during  FFT. 
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Figure  36.  Transducer  Meter  Output  Showing  C 
and  with  Carrier  Frequency  Filter 


ng  Pr'prv'^^y:  100  Hz 
NOISES 


Figure  38.  FFT  of  Noise  with  Air  on.  Major  Source  of  Noise  is 
Natural  Frequency  of  Torsion  Bar  at  100  Hz  Sampling 
Frequency. 


The  filtering  technique  can  be  better  explained  by  using  an 
example,  a  second-order  response.  A  second-order  response  , shown  in 
Figure  39,  with  a  time  constant  of  1  sec  and  *  0.4  was  generated.  A 
sampling  frequency  of  100  HZ  was  used  and  then  20  Hz  and  30  Hz  noises 

were  superimposed.  The  data  was  fitted  with  a  5th-order  polynomial  as 

shown  in  Figure  40.  The  deviatory  data  (Figure  41),  which  is  the 
difference  between  the  raw  data  and  the  fitted  data,  was  transformed  to 
the  frequency  domain  by  the  FFT.  In  the  frequency  domain,  the  spectral 
density  was  evaluated  and  analyzed.  Figure  42  shows  that  there  are  3 

major  frequencies,  one  is  at  30  Hz,  another  is  at  20  Hz,  and  the  last 

one  is  the  low  frequency  contribution.  The  low  frequency  includes  the 
natural  frequency  of  the  system  and  the  frequency  generated  by  the  curve 
fitting.  A  band-reject  filter  with  rejected  frequencies  of  20  Hz  and  30 
Hz  was  then  used  to  filter  out  the  noise.  The  filtered  deviatory  data 
in  the  time  domain  is  shown  in  Figure  43.  There  are  some  undesirable 
vibrations  at  both  ends  of  the  filtered  data.  This  is  due  to  errors  in 
spectral  computation  Introduced  by  truncation  and  is  known  as  "leakage”. 
The  filtered  deviatory  data  is  then  readded  to  the  polynomial  fit  to 
construct  the  final  curve,  it  was  found  that  the  filtered  data  and  the 
original  data  essentially  coincided  except  at  the  very  ends  (Figure  44). 
The  deviation  at  both  ends  is  due  to  the  leakage  of  discrete  FFT  as 
mentioned  earlier. 

Since  thi  Initial  response  of  the  curve  is  important  to  us,  further 
studies  were  lone  to  eliminate  the  leakage.  Unfortunately  leakage  is  a 
natural  consequence  of  using  a  FFT  and  cannot  be  avoided  even  if  a 
Hanning  window  is  applied.  The  only  possihle  way  to  avoid  the  leakage 
in  the  initial  response  was  to  move  the  initial  response  away  from  the 
starting  point,  l.e.,  sampling  at  t<0. 

The  same  set  of  second-order  response  data  was  used  again;  but 
this  time  was  shifted  to  the  right  by  starting  the  response  at  t-1.0 
min.  The  same  technique  was  used  with  the  same  curve  fitting  and 
filter.  Clearly,  we  can  reconstruct  the  original  second-order  response 
without  the  leakage  problem  by  pushing  the  origin  away  from  the  desired 
point  (Figure  45).  This  technique  was  very  helpful  in  analyzing  the 
initial  response  in  our  data. 

Curve  fitting  techniques  were  further  studied  to  find  the  best  way 
of  obtaining  the  deviatory  data.  Step  changes  and  polynomials  of 
different  degrees  were  used  to  fit  the  data.^A  set  of  transient  data  on 
the  Exxon  708-61  at  a  shear  rate  of  0.55  sec'was  used  for  the 
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comp ar s ion.  Polynomials  with  degrees  of  5,  10,  and  15  were  used.  It 
was  found  that  they  give  different  devlatory  results  and  different 
spectral  density  functions.  As  expected,  they  have  the  same  noise 
frequencies:  30  Hz  and  48  Hz  (Figures  46  to  48).  The  major  difference 
between  them  is  the  low  frequency  response:  the  frequency  that  is 
generated  by  the  curve  fitting.  Therefore,  it  is  desirable  to  fit  the 
data  with  polynomials  of  the  lowest  possible  degree.  The  only 
limitation  to  using  the  low  order  polynomial  is  in  minimizing  the 
leakage  during  Fourier  Transform.  For  example,  for  the  same  set  of 
data,  the  "leakage"  was  more  serious  on  both  ends  when  no  curve  fitting 
was  done  on  the  data,  as  shown  in  Figure  49.  With  the  low-pass  filter 
with  a  10  Hz  cut-off  frequency,  the  filtered  data  using  different 
fitting  polynomials  are  shown  in  Figures  50  to  52.  There  is  little 
difference  between  them.  For  a  better  comparison,  a  set  of  raw  data, 
fitted  data,  and  filtered  data  are  plotted  in  Figure  53. 


A. 5  Experimental  Procedures 


All  measurements  were  conducted  as  Follows:  precool  the  apparatus 
to  the  desired  temperature,  set  the  zero  gap,  load  the  sample,  then  wait 
at  least  3  hours  at  the  set  temperature  before  data  acquisition. 


A. 5 . 1  Temperature  Setting 

In  this  program,  a  wide  range  of  temperatures  was  studied,  -25  C  to 
25  C.  Measurements  at  room  temperature,  25  C,  were  controlled  by 
air-conditioned  room  temperature.  Measurements  other  than  room 
temperature  were  controlled  by  the  cooling  medium  circulating  in  the 
chamber  which  surrounds  the  plate  and  cone. 

For  the  low  temperature  measurements  the  following  procedures  were 
followed.  The  cryogenic  unit  was  turned  on  an  hour  in  advance  by 
setting  the  main  switch  to  "<-0  C".  It  stabilized  at  about  -60  C  in  an 
hour  with  the  external  circulation  valve  closed.  The  circulation  valve 
was  then  opened  to  allow  the  cooling  medium  to  circulate  around  the 
plate  and  cone.  At  the  beginning,  a  large  temperature  difference  was 
used  to  force  the  system  to  cool  faster.  When  the  temperature  of  the 
plate  and  cone  approached  the  desired  temperature,  the  temperature  of 
the  refrigerator  was  reset  to  allow  the  heat  transfer  in  the  chamber  to 
just  balance  the  heat  loss.  The  temperature  difference  at  this  stage 
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was  10  to  20  C  depending  on  the  experimental  temperature.  The  nitrogen 
valve  was  also  opened  to  let  the  gas  flow  at  a  rate  of  around  0.5  SCFM. 
Minor  adjustment  could  be  made  by  changing  the  flow  rate  of  the  cooling 
agent  and  the  fLow  rate  of  the  nitrogen  gas  in  the  chamber  to  keep  the 
temperature  of  the  sample  constant.  The  nitrogen  helped  circulation  and 
prevented  the  condensation  of  water  on  the  plate  and  cone. 

To  further  reduce  the  interaction  between  the  chamber  and  room,  a 
plastLc  wrap  was  used  to  cover  the  chamber  to  minimize  the  air  entering 
from  the  bottom. 


A. 5. 2  Loading  Sample 

After  the  plate  and  cone  reached  the  desired  temperature,  the  gap 
setting  reference  measurement  was  determined  as  describe!  in  the 
manual (83) . 

Gap  setting  was  much  more  difficult  at  tie  low  temperatures  than  at 
room  teraperatur>  .  Studies  ware  undertaken  to  correlate  the  change  in 
the  zero  gap  with  temperature.  Repeat  runs  were  measured  by  decreasing 
then  increasing  the  temperature.  The  re.- ults  were  not  reproducible  at  a 
fixed  temperature.  There  w as  always  a  jump  at  zero  temperature.  By 
disregarding  the  zero  temper  it ure  jump  and  bracketing  the  gap  size 
change  into  two  regions, ^it  vas  found  that  the  change  of  gap  size  with 
temper  it  ure  was  0.05x10“  i  ich  per  degri e  ceitrlgrade.  For  better 
reliability,  the  zero  gap  wa ;  determined  individually  at  each 
temper  it  ure. 

Tie  sanple  was  loaded  a>  quickly  as  possible  and  al ' owed  to  res:  at 
least  )  hours.  The  sample  hid  to  sir  a  .‘uffl-tent  amount  of  time  to 
recover  to  its  ground  state.  Different  test  tours  were  i ested  at  a 
shear  -ate  )f  27.6  sec-l  and  ic  was  found  tha  at  least  3  hours  was 
needed  to  obtain  a  reproducible  result  (v>ithin  the  noise  limits  as 
caused  by  t  ie  n;  tural  frequency  of  the  measur  ing  torsion  bar)  in 
transient  response  as  shown  in  Figure  54. 


103 


A. 5 . 3  Data  Acquisition 

A  software  package  written  by  Kibler  and  Mohler  was  used  to  acquire 
data.  The  programs  were  written  in  an  interactive  mode  which  made  the 
data  sampling  flexible.  The  flow  chart  of  the  operation  is  shown  in 
Figure  55. 

After  the  sample  reached  and  rested  at  the  desired  temperature,  the 
following  steps  were  used  to  sample  the  data. 

(1)  Run  "PASSTHRU"  to  initiate  the  LSI  11/02. 

(2)  Load  the  "absolute  loader"  and  "LSIPGM.OBJ"  to  the 
LSI  11/02. 

(3)  Set  up  a  frequency  table  on  the  LSI  11/02,  to  ready  the  LSI 
to  acquire  data. 

(4)  Select  the  required  shear  rate  by  the  gear  box  setting. 

(5)  Turn  on  the  motor  and  set  the  electromagnetic  clutch  to 
"brake"  position. 

(6)  Switch  clutch  from  "brake"  to  "drive"  position,  simultaneously 
turning  on  the  data  acquisition  switch  of  the  LSI. 

The  data  are  taken  by  LSI  automatically  based  on  the 
table  setup,  and  ter ninatatlon  occurs  either  at  the  end 
of  select  id  time  or  when  the  switch  is  turned  off. 

(7)  Transfer  the  data  fro,i  the  LSI  to  the  VAX. 

(8)  Repent  (3)  to  (7)  until  finished  collecting  data. 

For  most  measurements,  1000  Hz  was  used  to  sample  thi  first  part  of 
stress  growth  or  stress  relaxation,  and  100  Hz  was  used  for  steady  state 
measure  lent . 
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.SIG  l.OBJ 
'ASS  NRU.EXE 


Figure  55.  Flow  Chart  of  Da  a  Acquisition 


B.  Results 


B.  1  Structural  Viscosity 

A  simple  test  was  conducted  on  each  material  to  see  if  there  was  a 
structural  viscosity.  Each  sample  was  rested  for  2  hours  before 
testing.  Then,  the  sample  was  sheared  under  a  constant  shear  rate,  an 
increasing  shear  rate  and  then  a  decreasing  shear  rate.  The  shear  rate 
was  changed  fait  enough  so  that  the  material  did  not  have  time  to  relax 
to  its  ground  state.  It  was  found  that  there  were  hysteresis  loops  for 
Sun  Tech809-991 ,  Sun  Tech  839-988,  and  Exxon  708-61  materials,  which  show 
that  these  materials  do  have  structural  viscosities.  Clearly,  the 
materials  take  time  to  recover  their  original,  states;  i.e.,  the 
relaxation  tim>  of  materials  is  not  neg  igible.  Suntec  i  839-919  had  no 
hysteresis  loo;  .  The  viscosity  of  Tun  T?ch  839-919  decreased  with  snear 
rate  out  was  ir dependent  o"  the  shearing  history.  The  -esults  are  show  i 
in  Figure  36. 

\n  easy  aid  general  technique  (88)  was  ised  to  dem  .nstrate  the 
struc  ural  vlsi  isity  and  thixotropy  of  S jn  Tech839-938 .  Shear  rates 
were  acre  lsed  stepwise  from  0.11  t>  341  sec  ^  and  tie  sample  was 
sheared  at  each  shear  rate  for  abou:  15  seconds.  The  shear  stress  was 
r  ico c  led  by  the  LSI  11  microprocessor  at  a  r.te  of  13  H  .  The  mate-ial 
showed  tlm  ;-dep cadency  over  the  ent  .re  shear  rate  range.  At  high  shear 
rates,  thicotrooic  behavior  was  observed;  i.e.,  the  viscosity  decreased 
with  i  Line.  At  Low  shear  rates,  upon  sudden  start  up,  the  viscosity 
increased  to  a  laximun  then  decreased  to  an  asymptotic  value.  After 
cessation  of  shear,  the  viscosity  did  not  return  to  xero  because  of  the 
yield  stress  of  the  material.  The  response  curves  at  different  shear 
rates  are  shown  in  Figure  57. 

After  the  sample  was  sheared  a:  the  highest  rate  of  348  sec  ,  the  speed 
was  reduced,  step  by  step.  The  sample  was  then  rested  for  a  short  time, 
about  5  minutes,  and  the  shearing  was  repeated  ascending  from  the  lowest 
shear  rate.  As  seen  in  Figure  58,  after  each  abrupt  decrease  in  shear 
rate,  a  small  increase  in  stress  was  observed.  The  initial  peak  in  the 
repeat  shearing  was  not  as  high  as  in  th i  original  test.  Figure  58  also 
shows  that  the  longer  the  rest  time,  the  higher  the  initial  peak.  This 
type  of  study  provided  a  preliminary  picture  of  Sun  Tech  839-988:  it  is 
a  material  with  structural  viscosity  and  a  yield  stress.  The  magnitude 
of  yield  stress  can  be  determined  from  the  basic  shear  diagram. 
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B . 2  Basic  Shear  Diagram 


All  of  the  measurements  were  taken  based  on  the  new  technique 
previously  mentioned;  i.e.,  the  sample  was  covered  with  JP-10  while 
resting  to  avoid  the  skinning  effect.  At  high  shear  rates,  the  excess 
JP-10  was  removed  to  minimize  any  edge  effects.  Equilibrium  stresses 
were  obtained  by  shearing  the  material  until  reaching  steady  state, 
beginning  with  the  lowest  shear  rate.  The  procedure  was  repeated  at 
high  shear  rates  until  the  highest  shear  rate  was  obtained.  All  of  the 
equilibrium  stresses  were  digitalized  at  a  rate  of  100  Hz.  Repeated 
runs  were  conducted  on  each  material  to  obtain  reliabLe  data.  The 
results  were  quite  reproducible  and  were  used  to  construct  the  basin 
shear  diagr  uis  (BSD). 


B . 2 . 1  Roo  Temperature  Results 

The  BS  is  of  slurry  materials  (Sun  Tech809-995,  Sun  Tech  839-919, 

Sun  Tech839  988,  and  Exxon  708-61)  at  room  temperatur  ■  are  s  lown  in 
Figures  59  o  62.  The  equilibrium  str  sses  are  givei  a  Table  XVI.  Toe 
plots  of  th  apparent  viscosity  as  a  function  of  the  s  ' mar  rate  for  the 
four  slurri  s  are  shown  In  Figure  63. 

f)ver  t>e  range  of  shear  cates  studied,  the  four  slurries  are 
n m-Newtonl  n  with  shear-th i nning  behavior.  Except  for  Sun  lech839-919, 
the  slurrie  are  yield  materials.  Though  Sun  Tech839-J.9  showed  a 
t  indeucy  to  increase  viscosity  at  high  shear  rates,  it  is  di  ficult  to 
say  i:  the  : ocrease  was  due  to  dilantancy  or  instability. 


B.2.2  Low  Temperature  Results 

Since  | art  of  the  temperature  effect  on  the  viscosity  of  slurri  s 
comes  from  the  solvent,  the  viscosity  of  JP-10  was  st uiied  from  25  C  ro 
-20  C,  as  shown  in  Figure  64  and  Table  XVII.  The  viscosity  of  JP-10  did 
not  change  much  with  temperature.  The  viscosity  increased  from  2.76  cp 
to  9.25  cp  as  the  temperature  dropped  from  25  C  to  -20  C.  The  change  of 
viscosity  of  JP-10  with  temperature  can  be  described  !>/  an  Arrhenius 
type  equation: 


n  -  0.0023545  *  [exp( 2098.348/T) ] 
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Table  XVI.  Equilibrium  Stresses  of  Slurries 


rate  stress 

919  995  988 


0.00348 

0.27 

0.00438 

0.32 

0.0055 

0.37 

0.00692 

0.42 

0.0087 

0.49 

0.011 

0.57 

0.0138 

0.65 

0.0174 

0.75 

0.02  76 

0.97 

0.0438 

1.25 

0.0692 

1.60 

0.11 

0.174 

2.53 

0.276 

3.16 

0.438 

0.692 

4.92 

1.1 

6.19 

1.74 

7.83 

2.76 

10.02 

4.38 

13.03 

6.92 

17.17 

11.0 

23.13 

17.4 

31.69 

27.6 

44.44 

34.8 

53.09 

43.8 

63.68 

69.2 

92.98 

110.0 

138.6 

174.0 

209.52 

2  76.0 

321.61 

348.0 

400.22 

438.0 

497.82 

550.0 

618.14 

692.0 

768.6 

870.0 

953.6 

6.47 


7.44 

8.96 

9.98  116.3 

11.21  116.6 

12.71  116.8 
117.2 

16.71  117.9 
118.8 

22.73  120.9 

124.4 
32.07  129.6 

136.7 
47.17  146.3 

58.24  158.8 

73.04  175.3 

185.4 

196.8 

120.52  224.6 

159.69  261.8 

215.67  310.6 

299.04  376.8 

418.7 

426.22  467.4 

513.15  524.1 

623.26  591.8 

763.56 


708-61 


2.33 
2.61 
2.90 
3.20 
3.51 
3.83 
4.14 
4.46 
5.10 
5.76 
6.45 
7.19 
8.01 
8.96 
10.08 
11.43 
13. II 
15.22 
17.91 
21.40 

25.90 

31.90 
39.76 
50.29 

64.39 

83.07 

108.4 

141.8 

186.1 

213.3 

244.1 
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Table  XVII.  Viscosity  of  JP-10  at  Different  Temperatures 


Temp(#C)  Viscoslty(cp)  Viscoslty(cp) 

Measured  Calculated 


24.5 

2.7564 

2.7135 

10. 0 

3.8954 

3.8932 

0.0 

4.9886 

5.1067 

-10. 0 

6.72  56 

6.83  79 

-11.0 

7.0249 

7.0490 

-19.5 

9.4395 

9.2177 

-19.5 

9.2486 

9.2177 

Standard 

Deviation  ■ 

1.82  *  10**  (-2) 

These  results  check,  well  with  the  resuLts  recenlLy  provided  by 
Wrighc-Patterson  AFB. 

Only  selected  siiirries  were  studiec  at  low  temperatures.  Sun  Tech 
309-595  was  studied  because  of  its  ease  of  handling.  Sun  Tech839-988 
and  Exxon  708-61  were  measured  at  low  temperatures  due  to  their 
prominent  yield  stresses. 

The  BSPs  of  Sun  Tech809-995  at  temperatures  of  25  C,  0  C,  -5  C 
-if)  Cj  and  -20  C  are  shown  in  Figure  br.  Flow  curves  at  different 
temperatures  have  the  same  shape  and  soift  to  the  right  as  temperature 
decreases;  that  is:  the  lower  the  temperature,  the  higher  the  viscosity. 
Exxon  708-61  was  also  studied  at  1  C  and  ~tJ  C  as  seen  in  Figure  66.  As 
expected,  the  viscosity  increased  with,  a  decrease  in  temperature  and  the 
BSD  shifted  to  the  right  as  the  temper -t ure  decreased.  4  strange 
phenomenon  was  observed;  the  viscosit.-  M  :o.  change  much  with 
temperature  at  low  shear  rates.  This  ;  !  !  discussed  later. 

More  interesting  results  were  obr  lina’  '  •_ ur.  Tech  839-988  at  low 

temperatures.  The  BSDs  of  f  ur.  Tech.339  8  -t  temperatures  of  23,  9,  0, 
-1C,  -20,  end  -25  C  are  shown  in  Figur  67.  One  interesting  result  was 
that  the  viscosity  increased  or  decree  i  v.i.o  temperature  decreased 

depending  upon  the  shear  rate.  At  Lov  sh-t  t  rates,  the  viscosity 
decreased  as  cemperature  decreased,  v.T. : -  at  high  shear  rates,  the 
viscosity  increased  as  temperature  dec  The  shear  rate  of  15  sec 

was  the  break  point  for  all  temperatur  .  This  will  be  discussed  in  the 
next  section. 


3.2.3  Yield  Stress 

BSDs  in  Figures  59  to  62  show  that  Suri  Tech809-995,  Sun  Tech 
839-988,  and  Exxon  708-61  are  yield  raa  erials.  The  existence  of  a  yield 
stress  complicated  the  measurement  at  ow  shear  rates.  In  turn,  the 
yield  stress  was  difficult  to  measure  because  very  low  shear  rates  were 
not  attainable  and  b  cause  the  instrument  lacked  enough  sensitivity  to 
measure  the  shear  stress  at  low  shear  rates. 

The  yield  stress  can  be  determine.:  from  the  vertical  section  in  the 
log-log  plot  of  the  BSD  or  by  extrapol  itlon  of  the  BSD  to  zero  shear 
rate.  The  extrapolation  is  accomplish.  1  by  fitting  the  data  at  '.ow 
shear  rates  with  a  polynomial  with  a  minimum  deviation,  then 
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extrapolating  the  polynomial  to  zero  shear  rate.  Yield  stresses 
obtained  by  both  methods  are  listed  in  Table  XVIII.  Sun  Tech  809-995  qpd 
Exxon  708-61  both  had  a  small  yield  stress,  on  the  order  of  5  dyne/cra^. 
The  yield  stress  of  Sun  Tech  839-988  was  appreciable,  around  110 
dyne/cm  \  The  yield  stresses  of  Sun  Tech  839-988  at  different 
temperatures  are  given  in  Table  XIX. 


B  2.4  n0  and  n^ 

%  and11”  are  lower  and  upper  Newtonian  viscosities.  For 
materials  without  a  yield  stress,  such  as  Sun  Tech839-919,  %  and  Too 
were  determined  from  the  apparent  viscosity  at  low  and  high  shear  rates, 
respectively.  For  yield  materials,  D0  is  infinite  and  n°°  %  as 

}  -+  °°  ,  where  np  is  the  plastic  viscosity  or  Bingham  viscosity. 
Therefore,  was  taken  as  the  plastic  viscosity  which  was  determined 
from  the  slope  of  the  linear  portion  of  the  BSD  plots  at  high  shear 
rates.  Table  XX  gives  n0  and  ^  of  different  materials  at  room 
temperature.  The  change  of  h®  with  the  temperature  of  Sun  Tech839-988 
is  given  in  Table  XXI. 


B.3  Transient  Response 

Before  model  constants  were  evaluated,  transient  data  were  all 
pretreated  in  the  following  sequence:  filtered  to  remove  noise, 
adjusted  for  top  plate  movement,  and  normalized  on  the  equilibrium 
value. 

First,  a  Newtonian  fluid  was  studied.  As  expected,  the  transient 
response  of  the  viscosity  was  a  step  change  after  a  step  change  was 
imposed  in  the  shear  rate  (see  Figure  68).  This  response  was 
independent  of  shear  rate.  The  study  of  the  Newtonian  fluid  also 
confirmed  that  the  response  time  of  the  instrument  is  fast  enough  to 
study  the  transient  response  of  slurry  materials. 

The  transient  response  of  Sun  Tech  839-919  and  Sun  Tech  809-995  were 
similiar  to  Newtonian  fluids.  They  showed  little  time  dependency  except 
at  high  shear  rates.  Typical  transient  behavior  of  these  two  materials 
is  shown  in  Figures  69  and  70.  Since  these  materials  showed  fast 
response,  they  were  not  studied  further  by  this  technique. 


Table  XVIII.  Yield  Stresses  of  Slurries 


Yield  Stress  (dyne/cm**2 ) 

material  log-log  plot  Cur/e  fitting 


Sun  Tech90 9-995  8.0 

Sun  Tech 839-988  110.0 

Exxon  708-61  6.0 


5.0 


1 10.0 
2.0 


Table  XIX.  Yield  Stresses  of  Sun  Tech839-988 
at  Different  Temperatures 

Temper ature(C)  Yield  Stress ( dyne/cm**2 ) 


Figure  69.  Transient  Response  for  Sun  Tech  839-919  at  Low  Shear. 


Figure  70.  Transient  Response  for  Sun  Tech  809-995  at  Low  Shear. 


« 


The  transient  response  of  Exxon  708-61  is  shown  in  Figure  71.  The 
material  showed  a  small  peak  at  low  shear  rates  and  behaved  as  a 
thixotropic  material. 

The  transient  response  curves  for  Sun  Tech  839-988  at  shear  rates 
from  0.055  to  2.76  sec-1  are  shown  in  Figure  72.  The  curves  show  solid 
property  characteristics  and  structural  viscosity  at  low  shear  rates. 
Even  at  low  shear  rates,  overshoot  was  observed. 

The  original  and  corrected  data  for  top  plate  movement  at  two 
different  shear  rates  are  shown  in  Figure  73  and  74.  The  correction  was 
minor  and  did  not  adversely  affect  the  initial  slope.  The  peak  values 
were  slightly  shifted. 


8.4  Kinetic-Elastic  Model 


The  kinetic-elastic  model  was  modified  to  describe  the  flow 
behavior  of  slurry  materials  for  both  steady  state  and  transient 
conditions.  The  yield  stress,  Ty  ,  was  Incorporated  into  the  equation. 
Different  approaches  were  tried  and  the  best  result  was  obtained  by 
replacing  values  of  the  T  in  the  equation  by  Ty  -  T  ,  which  is 
equivalent  to  a  shift  in  the  axis  by  Ty>  as  shown  in  Figure  75.  For 
convenience,  t'  is  used  to  represent  t  -  ty  in  the  following  section. 
Actually,  this  is  a  generalised  definition  of  the  kinetic  model.  For  a 
naterial  with  no  yield,  t  =  0,  the  equation  reduces  to  the  original 
form.  The  modified  equation  can  be  applied  to  any  system  with  or 
without  yield  stress. 


B.4.1  Definition  of  Structural  Parameters 

The  structural  parameter,  F,  is  a  function  of  v.scosity.  Different 
functions  have  been  tested  to  define  F.  Though  different  definitions  of 
F  can  be  used  with  the  kinetic  theory  to  correlate  the  steady  state 
data,  some  of  them  fail  to  describe  the  solid  behavior  at  the  low  shear 
rates  in  the  transient  response.  The  inverse  level  rule  and  the 
Rlchardson-Zaki  equation  were  used  depending  on  the  nature  of  the 
slurry. 
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Transient  Response  for  Exxon  708-61  at  Low  Shear. 
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Figure  73.  Data  Correction  for  Top  Cone  Movement. 


ORIGINAL  £  ADJUSTED  DATA 


rection  for  Top  Cone  Movement. 


For  a  slurry  without  a  yield  stress  and  having  n0  and  the 
inverse  level  rule  is  used  to  define  F.  F  is  defined  as 


(17) 


where  a=l. 


The  Richardson-Zaki  equation  is  used  to  define  F  for  the  yield 
materials:  Sun  Tech  839-988,  839-995,  and  Exxon  708-61.  In  this  case,  F 
is  defined  as 


F 


1  -  (n  /n) 

oo 


1/2.65 


(18) 


B.4.2  Steady  State 


The  kinetic  equation  alone  is  used  to  correlate  the  steady  state 
viscosity.  At  steady  state,  the  kinetic  equation  can  be  rearranged  to 


Kip  -  xn-Pz  -  c"-» 

*  ( 19) 

where  C  is  the  concentration  and  m=l,  n=2  are  used.  K  and  P  were  found 
by  a  least  squares  fit  of  log  [(1-F)2/F]  versus  log  T  .  The  optimized 
equilibrium  constants  (P's  and  K's)  for  different  slurries  are  tabulated 
in  Table  XXII.  The  experimental  data  and  the  predicted  values  are  shown 
in  Figures  76  to  79.  The  K's  and  P's  of  Suntech  839-988  at  different 
temperatures  are  summarized  in  Table  XXIII. 
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Table  XXII.  K  and  P  of  Slurries  by  Richardson-Zaki 
Material  K  P 


Suntech  809-995  0.010670  0.9991 

Suntech  839-919  0.110107  0.89086 

Suntech  839-988  0.014745  0.78287 

Exxon  708-61  0.003988  0.9935 


Table  XXIII.  K  and  P  of  Suntech  839-988  at 
Different  Temperatures 


Temp  (C) 

K 

P 

23.0 

0.014745 

0. 78287 

9.0 

0.02270 

0.7184 

0.0 

0.03633 

0.6719 

-10.0 

0.047614 

0. 73443 

-20.0 

0. 10424 

0.66912 

-26.0 

0.01899 

1. 1031 

Sun  Tech8D9-995  F.qu  i  1  ihr  lum  Data 


Figure  77.  Kinetic  Theory  Fit  for  Sun  Tech  839-919  Equilibrium  Data. 


Fit  for  Slin  Tech  839-988  Equilibrium  Data. 


B . 4 . 3  Transient  Response 


Transient  response  is  a  combinative  result  of  thixotropy  and 
elasticity.  The  thixotropic  constants,  kj ,  p^  ,  k^,  p2  and  elastic 
modulus  Gj  were  evaluated  from  transient  data.  Sun  Tech839-988  was  used 
for  this  study. 

For  convenience,  the  kinetic  and  elastic  equations  with  structural 
parameter  F  are  listed  again. 

P1  p2  2 
-dF/dt  =  kjT  F  -  k2t  (1  -  F) 


(nt/Gi^  +  T  =  n  y 


(20) 


The  parameter  G  can  be  evaluated  from  equation  (2)) 

V 

G.  =  - 1 - 

1  ntY  -  t 

at  the  initial  time,  Y  nt  >>  t  so 

T 

G-  ~~m 
1  Y 


(21) 


(22) 


The  values  of  G^  at  different  shear  rates  are  tabulated  in  the  Table 

XXIV.  The  values  of  k^  and  p^  were  calculated  from  the  initial  stress 
change,  l.e.,  the  stress  growth  near  /.  tco  time.  Once  kj  and  p^  are 
determined  ,  k2  and  P2  can  be  obtained  from  the  equilibrium  constants  K 
and  P,  i.e.,  k2=*  k]/K,  and  p2=  P2“P*  A  set  oC  these  constants  obtained 
for  Sun  Tech839-988  at  a  shear  rate  oE  2.76  sec-lare  summarized  in  Table 

XXV.  All  the  constants  were  obtained  by  the  best  fit  of  the  data  using 
an  Advanced  Continuous  Simulation  Language  (ACSL)  simulation  package. 

The  single  Gi  obtained  from  the  initial  response  cannot  describe 
the  whole  transient  response,  as  shown  in  Figure  80.  G]  is  thus  defined 
as  a  function  of  time  as  already  proposed  for  polymers.  With  a  variable 
Gi,  the  experimental  data  and  the  predicted  transient  response  are 
compared  in  Figure  81. 


K  inet  ic  -F.last  ic  Theory  Fit  Using  a 
V.Jlue  of  Gl  at  a  Shear  Rate  of  2.7S 


Figure  81.  Kinetic-Elastic  Theory  Fit  Using  Variable  Cl  with 
Model  (II)  at  a  Shear  Rate  of  2.76  sec"  . 


C.  Discussion 

The  main  purpose  of  this  study  was  to  obtain  a  better  understanding 
of  the  rheological  properties  of  carbon  black  slurries  and  to  extend  the 
kinet ic -el as  tic  model  to  predict  their  flow  behavior.  The  discussion 
focuses  on  these  two  aspects. 


C . 1  Shear  Dependenc e  of  Carbon  Black  Structure 

The  results  showed  that  all  of  the  carbon  black  slurries  tested  are 
shear-thinning  materials.  This  suggests  that,  as  in  other  suspension 
systems,  at  reasonably  high  concentrations  carbon  black  has  a  structure 
which  can  be  broken  down  under  shear.  The  structure  could  be  aggregate 
or  a  network  which  is  groups  of  particles  held  together  by  attractive 
forces.  The  formation  of  structure  results  in  a  larger  effective  volume 
fraction  than  the  actual  volume  frac  ion  of  particles.  The  higher 
effective  volume  fraction  and  the  ini eraction  between  particles  increase 
the  resistance  to  flow,  and  this  causes  the  much  larger  viscosity  of  the 
suspension  than  of  the  solvent. 

In  a  shear  field,  the  structure  is  broken  down  by  the  shear  forces. 
The  degree  to  which  the  structure  is  broken  down  depends  on  the 
magnitude  of  shear  force  and  the  strength  of  the  bonds  between  the 
particles.  As  the  shear  force  increases,  more  of  the  structure  will  be 
broken  down.  At  each  level  of  shear  rate,  an  equilibrium  structure  can 
exist  corresponding  to  an  equilibrium  viscosity.  The  effective  volume 
fraction  decreases  with  an  increase  in  shear  rate  as  does  the  viscosity. 
This  accounts  for  the  shear-thinning  behavior  of  these  slurries. 

Moreover,  at  high  enough  shear  rates,  the  structures  can  be  broken 
into  non-interactive  units  (not  necessarily  individual  particles  but  any 
basic  flow  unit)  and  the  fluid  behaves  as  a  Newtonian  fluid.  This  is 
commonly  called  the  upper  Newtonian  region.  In  this  study,  even  at  the 
highest  shear  rate,  nearly  1000  sec-T,  the  upper  Newtonian  region  was 
barely  achieved. 

Hysteresis  loops( Figure  56)  and  thixotropic  behavior  in  the  time 
domain  confirmed  the  existence  of  a  shear  rate  dependent  structure.  In 
the  hysteresis  loop  study,  with  an  i  icreasing  shear  rate( upcurve) , 
structure  is  broken  down.  As  the  shear  rate  decreases(downcurve)  the 
structure  builds  up.  The  structure  that  exists  at  each  shear  rate(in 


Che  downcurve)  is  ciose  Co  Chac  ac  Che  highesc  shear  raCe,  unless  Che 
scrucCure  build-up  is  very  fasc.  Therefore,  differenc  sCrucCures  can 
exisc  aC  the  same  shear  rate  as  a  result  of  different  history.  This  has 
been  shown  by  Jones  and  Brodkey(4).  They  constructed  flow  curves  of  a 
thixotropic  material  at  different  structural  levels.  The  structural 
curve  would  be  identical  to  the  downcurve  if  the  time  scale  of 
measurement  is  fast  enough.  They  found  there  was  an  equilibrium 
structure  associated  with  each  shear  rate  and  that  different  structures 
could  exist  at  the  same  shear  rate  depending  on  the  shear  history.  As  a 
result  different  viscosities  are  measured,  which  results  in  the 
hysteresis  loop.  During  transient  response  at  constant  shear  rate,  the 
viscosity  decreases  with  time  as  cause!  by  the  breaking  of  the 
structure.  The  structure  can  recover  its  ground  state  at  rest.  The 
rebuilding  of  the  structure  can  be  demonstrated  by  repeating  the 
transient  experiment.  After  a  sample  is  sheared  to  its  steady  state 
structure,  stop  the  shearing,  pause  for  a  few  minutes,  and  repeat  the 
test.  The  time  dependent  behavior  is  still  present,  except  that  the 
initial  peak  in  the  stress  is  smaller,  as  illustrated  in  Figure  82. 

The  hysteresis  loop  shown  in  Figure  56  is  often  used  as  a 
measurement  of  thixotropy  and  is  somet  mes  called  the  "thixotropic 
loop".  The  degree  of  thixotropy  is  of' en  related  to  the  area  between 
the  upward  and  downward  curves.  Although  this  measurement  can  describe 
a  few  systems,  it  is  usually  regarded  suspiciously  because  of  possible 
errors  arising  from  the  method  of  measurement.  The  area  in  the 
thixotropic  loop  is  dependent  on  the  procedure  used  for  the 
measurements.  It  different  paths  are  used  to  change  the  shear  rates, 
different  areas  Jill  be  obtained.  Even  if  the  same  shear  rate  path  is 
used,  different  shearing  times  will  result  in  different  areas  within  the 
loop.  The  instalment  can  also  contribute  errors  in  thixotropic  loop 
measurement,  for  example,  heating  effect,  change  of  apparent  rate  due  to 
the  movement  of  the  torque  spring,  and  inertia  of  the  moving  parts.  The 
area  within  the  loop  is  not  a  satisfactory  measurement  of  the  degree  of 
thixotropy  shown  by  a  material.  Therefore,  the  thixotropic  loop  studies 
made  here  were  only  used  as  an  indication  of  thixotropy  or  structural 
/iscosity.  Ihe  hest  way  to  describe  thixotropy  should  be  based  on  a) 

;he  equilibrium  curve  and  b)  the  rate  o  approach  to  equilibrium  after  a 
change  of  shear  rate. 
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C.  2 


Flow  Carves 


Umega(89)  suggested  that  the  rheological  beh  ivior  of  suspensions 
could  be  described  by  an  extended  Ostwald  flow  curve(Figurc  83)  if  a 
sufficiently  wide  ran,;e  01  shear  rates  was  studied.  The  flow  beha.'ior 
of  materials  progressively  passes  from  a  lower  Newtonian  region(A-b)  to 
a  pseudoplasr.ic  flow  region(B-C),  followed  by  an  upper  Newtonian 
region(C-D),  a  dilatant  region(D-E),  and  then  another  Newtonian 
region( E-F) .  In  this  study,  the  shear  rate  range  covered  was  not  high 
enough  to  include  the  whole  flow  curve;  only  the  lower  parL  of  the  curve 
was  observed. 


C.2.L  Measurements  and  ( bar  ac  ter  is;; ics  of  BSPs 

Some  ini  ere;- ting  rest  Its  were  observed  at  both  low  and  high  si  :ar 
rates.  They  will  be  disci-;:;  id  individually. 

Sun  Tech 309-995  .s  th  >  only  material  t  lat  ;how  d  little  skinni  g. 
The  BSD  of  Sun  Tech 80u-995  was  studied  with  and  without  JP-10.  The  e 
was  little  differ >nce  between  the  tw >  cases,  ex<  ept  at  high  shear  r  ites 
(see  Figure  59).  The  sample  wit  o  JP-10  lad  a  Inwer  viscosity  than  hai 
without  JP-iC  at  shear  rat  ;s  larger  than  276  see-1.  This,  nj  doubt,  was 
due  to  the  dilution  of  the  sample  by  JP-10  at  the  higher  shear  rate  A 
test  was  made  wit  lout  JP-1  )  starting  at  a  shear  rate  of  69.r?sec  a  d 
increasing  the  rate  to  692  sec_1and  was  carried  out  fast  enough  to  void 
the  skinning.  The  results  were  the  same  as  those  started  from  a  mu  h 

lower  shear  rate  without  J  '-10.  Therefore,  the  dilution  is  more  likely 

to  exist  at  the  high  shear  rate.  This  problem  was  avoided  by  removing 
the  excess  JP-10  for  the  h  gher  shear  rale  measurements. 

The  BSD  of  Sun  Tech  83')-919  was  studied  with  the  JP-10  covering  over 
tiie  shear  rates  range  of  0.00348  to  870  sec-1  (see  Figure  60).  One  .et 
of  earLier  data  with  the  s  .inning  problem  was  also  plotted  for 
comparison.  Three  runs  we' e  conducted  with  different  starting  shear 
rates.  There  was  a  difference  in  the  low  shear  rate  range  which  mi.,ht 
be  due  to  a  weak  structure  in  the  material.  However,  the  signal  at  low 

shear  rates  was  so  small  that  the  difference  could  be  the  result  of 

measurement  errors.  At  bo .h  low  and  ugh  shear  rates,  the  slope  of  the 
curve  was  close  to  1,  which  implies  tnat  the  material  was  Newtonian  at 
both  extremes,  i.e.,  it  had  high  and  low  Newtonian  viscosities. 
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Sun  Tech  839-919  was  the  material  which  showed  no  hysteresis  between 
shear  rates  of  1.1  and  276  sec- 'with  or  without  the  JP-10  covering. 

But,  at  shear  rates  above  276,  it  did  show  hysteresis,  as  shown  in 
Figure  84.  Apparently,  something  different  happened  with  this  material 
at  shear  rates  higher  than  276  sec-l  The  explanation  is  that  either  the 
material  has  some  strong  structure  that  is  broken  down  at  the  higher 
shear  rates  or  that  an  instability  generated  at  the  higher  shear  rate 
changes  the  nature  of  the  material.  There  is  little  possibility  that 
this  is  due  to  the  mixing  of  JP-10  with  the  sample  since  the  hysteresis 
loop  was  also  observed  for  the  measurement  without  JP-10  covering. 

The  BSD  of  Exxon  708-61  was  studied  at  room  temperature  from  •  = 

1 

0.00348  to  438  sec  .  The  data  branches  out  at  low  shear  rates,  as  seen 
in  Figure  62.  It  was  found  that  the  lower  the  initial  shear  rate  the 
higher  the  viscosity  obtained  at  any  specified  shear  rate.  There  are 
two  possible  reasons:  the  structural  property  or  rheopetic  property  of 
the  material.  There  could  be  a  weak  structure  which  was  broken  down 
during  the  loading  of  the  sample.  Such  weak  structures  could  build  up 
during  resting  or  under  mild  shear  rate.  If  the  structures  build  up 
under  a  mild  shear  rate,  the  material  is  rheopetic.  Otherwise,  the 
structure  builds  up  during  the  resting  period  between  runs.  Studies 
were  made  to  check  if  weak  structures  were  broken  down  during  the 
loading  of  the  sample  and  it  the  equilibrium  viscosity  at  the  low  shear 
rate  depended  on  the  resting  time.  Several  runs  were  carried  out  from 
v  =  0.0174  sec~^to  174  sec-iby  letting  the  sample  rest  from  )  to  3 
hours.  As  shown  in  Figure  85,  the  viscosity  at  the  low  shear  rate  does 
depend  on  the  resting  time.  A  large  difference  was ^ found  at  the  lower 
shear  rates  and  the  curves  coincide  at  •  >  10  sec  .  This  showed  that 

there  is  a  weak  structure  which  was  broken  easily  and  took  time  to 

rebuild.  Moreover,  as  mentioned  in  the  previous  chapter,  the  viscosity 
at  27.6  sec-was  independent  of  the  resting  time  which  is  consistent 
with  this  result.  In  order  to  obtain  the  "true"  viscosity  ai  low  shear 

rates,  we  would  have  to  let  the  sample  sit  for  a  very  long  period  to 

rebuild  its  weak  structure.  From  Figure  85,  the  raa  eriaL  would  n^ed  3 
hours  or  more  to  rebuild  such  a  weak  s.ructure  at  the  0.0174  sec-. 
Allowing  the  sample  to  rest  for  such  a  long  time  induced  skinning,  thus, 
we  wilL  always  have  difficulty  in  obtaining  viscosity  at  lew  shear 
rates.  In  addition,  from  Figure  85,  toe  shifting  of  the  curves  to  the 
right  as  the  rest  time  increased  implies  that  the  material  had  a  yield 
stress.  The  yield  stress  cannot  be  over  6  dyne/cnr,  because  any 
reasonable  extrapolation  of  the  180  mil  curve  would  give  a  value  less 
than  this.  The  transient  response  of  the  material  was  also  studied  to 


determine  the  rest  time.  The  transient  response  at  a  shear  rate  of  27.6 
sec-1with  different  resting  times  ranging  from  0.5  hour  to  6  hours  is 
shown  in  Figure  54.  Though  the  equilibrium  viscosity  at  y  =  21 .6  see  1 
was  independent  of  the  resting  time,  the  transient  response  did  depend 
on  the  resting  time.  Therefore,  letting  the  sample  rest  to  its  ground 
state  is  important  in  both  transient  and  steady  state  studies. 

The  viscosity  of  Sun  Tech  839-988  was  difficult  to  measure  due  to 
the  nature  of  the  material.  This  material  was  thicker  than  the  other 
slurry  materials,  especially  at  low  shear  rates,  which  complicated  the 
measurements.  The  same  technique,  covering  the  sample  with  JP-10,  was 
used  to  study  the  viscosity.  It  was  found  that  above  a  critical  shear 
rate  (around  1  sec  )  the  mixing  of  JP-10  and  the  slurry  became  serious 
and  the  viscosity  decreased  as  time  proceeded.  The  mixing  was  caused  by 
the  edge  effect.  In  Figure  86,  we  can  see  that  the  stress  was  constant 
>ver  a  wide  range  of  shear  rates.  In  this  region,  there  was  almost  no 
flow  inside  the  gap  because  of  the  yie  d  stress.  Above  a  shear  rate  of 
1  sec-1 >  the  structure  of  the  material  was  broken  and  the  material 
started  to  flow.  The  curve  bends  towad  lower  stresses  as  the  shear 
rate  increases  above  1  sec-l  . 

To  check  the  mixing  effect,  measurements  wer  *  made  without  JP-lo  to 
see  if  the  curve  had  the  same  character Lstics .  As  can  be  seen  in  Figure 
16,  they  did  not.  Addi'.iona!  runs  with  different  resting  time  were  ^ 
carried  out  to  see  if  the  curve  infleci ed  above  a  shear  rate  of  1  sec  . 
No  bending  was  observed  and  the  curves  shifted  to  the  right  as  the 
resting  time  increased  (Figure  87).  We  are  convinced  that  tie  decrease 
in  the  shear  stress  was  due  to  dilution  with  JP-10. 

A  systematic  experiment  was  then  made  to  study  the  resting  time 
effect  on  Sun  Tech  839-988  at  different  shear  rates.  The  sample  was 
shaken  vigorously  before  loading  to  insure  homogeneity.  Then,  the  shoar 
stress-time  data  was  obtained  every  hour. 

At  a  shear  rate  of  0.11  sec  ,  it  was  found  that  the  apparent 
;quilibrium  viscosity  was  a  function  o)  time.  However,  if  the  initial 
stress  value  was  used  rather  than  zero  as  the  ground  value,  the 
equilibrium  viscosity  was  independent  of  the  resting  time  as  seen  in 
Figure  88.  Nevertheless,  the  transient  curves  did  depend  on  the  resting 
time.  The  peak  values  of  the  transient  response  were  also  plotted  as  a 
function  of  time  (Figure  89).  They  increased  initially  and  t  len  leveled 
off  after  3  hours. 
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Shear  Stress  (dyne/crrr) 


Measurement  for  Sun  TechR39-98S  Shear  in 
Effect  of  Mix ine  of  JP-1 0  into  Somnlc. 


Figure  88.  Shifted  Equilibrium  Viscoeity  (Using  Initinl  Stress 


Figure  89.  Transient  Curves  ami  Shifting  of  Peak  Values 
as  a  Function  of  Resting  Time. 


Shear  rates  of  2.19,  5.5,  and  11)  sec  were  studied  in  the  same 
manner.  The  equilibrium  stress  as  a  function  of  time  is  shown  in  Figure 
90.  Except  for  the  shear  rate  of  2.19  sec-1  ,  they  all  showed  the 
tendency  of  decreasing  viscosity  with  time.  This  is  further  evidence 
that  dilution  by  JP-10  took  place  at  high  shear  rates.  Though  the 
equilibrium  viscosity  at  y  =2.19  sec-1  approached  an  asymptote  as  the 
time  increased,  its  peak  values  from  the  transient  responses  showed  a 
tendency  to  decrease  with  time  as  seen  in  Figure  91.  This  means  that 
the  mixing  of  JP-10  and  the  slurry  took  place  under  this  shear  rate  but 
not  seriously  enough  to  affect  the  equilibrium  value. 


Reproducibility  was  very  poor  when  the  data  was  taken  with  the 
sample  full  to  the  rim  of  the  plate  (Figure  92).  Due  to  the  high 
viscosity  of  this  material,  the  shear  rate  difference  between  the  sample 
in  the  gap  and  that  in  the  reservoir  ;  aused  mixing,  which  resulted  in 
erroneous  measurements .  To  avoid  the  edge  effect,  the  viscosity  of 
Sun  Tech839-983  was  measured  without  .  I  P—10  covering.  Nevertheless, 

JP-10  was  used  during  the  resting  period  and  was  removed  for  the  actual 
measurement.  Excess  sample  was  also  removed  from  the  reservoir  (Figure 
92).  The  results  were  reproducible  as  shown  in  Figure  62.  Due  to  the 
small  signals  and  low  rotational  rates,  the  data  scat  tered  at  low  shear 
rates.  In  spite  of  the  scattering,  the  plot  indicates  a  yield  stress  of 
around  110  dyne/cm2.  In  addition,  by  extrapolating  the  data  at  high 
shear  rates  to  zero  shear  cate,  a  Bingham  yield  stress  of  220  dyne/cm2 
is  obtained.  The  results  are  also  plotted  in  Figure  62  and  coincide 
with  the  measurement  using  JP-10  at  the  low  shear  rate.  In  order  to 
obtain  accurate  results  for  this  material,  the  edge  effect  must  be 
minimized . 


C . 2 . 2  Comparision  between  Slurries 


Though  all  the  slurries  contained  almost  the  same  loading  of  carbon 
black  (50  wt%) ,  they  possessed  different  rheological  jehavior, 
especially  at  low  shear  rates.  At  high  shear  rates,  the  viscosities  of 
the  different  slurries  were. not  much  different  and  ranged  from  50  cp  to 
100  cp.  The  difference  between  them  became  wider  as  the  shear  rate 
decreased  (Figure  63).  There  is  more  than  100  times  difference  between 
Sun  Tech839-988  and  the  others.  This  difference  was  due  to  the  particle 
size  and  surfactant  used.  These  two  factors  affected  the  formation  of 
structures  which  were  important  at  low  shear  rates.  Apparently,  Sun  Tech 
839-988  had  a  smaller  particle  size,  which  resulted  in  a  large  yield 
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Sure  <10.  Equilibrium  Viscosity  as  a  Function 

of  Rest  Time  at  Different  Shear  Rates. 


Viscosity  at  Di 


stress  at  low  shear  rates.  The  surfactant  was  used  to  stabilize  and  to 
increase  dispersity.  A  well  dispersed  system  will  favor 
particle-solvent  interaction  and  prevent  strong  structure  formation  . 
Therefore,  a  well  dispersed  system  will  have  weaker  structure  and  lower 
viscosity. 

C.3  Temperature  Effects 

According  to  Che  traditional  time-r ate-teraperat ure  superposition,  a 
decrease  in  temperature  corresponds  to  a  shift  to  higher  effective  shear 
rates  for  viscosity  versus  shear  rate  plots.  The  data  did  follow  this 
trend;  flow  curves  shifted  to  higher  shear  stress,  as  temperature 
dropped  (Figures  65  to  67).  Results  ;ere  also  consistent  with  the 
Ostwald  curve;  at  lower  temperatures  and  higher  shear  rates,  flow 
curves  shifted  to  the  right,  falling  nto  the  dilatant  region. 

The  phenomenon  observed  at  low  shear  rates  was  surprising  and  ha  I 
not  been  noted  in  the  literature.  At  temperatures  below  zero,  the  lov 
shear  viscosity  showed  a  systematic  d.  crease  as  temperature  decreased 
(Figure  67).  Though  the  result  was  u  usual,  it  can  be  explained  on  the 

basis  of  structural  viscosity.  Due  t  particle  interactions,  there  was 

a  structure  which  resulted  in  both  a  >  ield  stress  and  thixotropy.  The 
linkage  that  formed  the  structure  depended  on  the  kinetic  energy  and  the 
Brownian  motion  of  particles.  As  the  temperature  decreased,  the  linkage 
became  weaker,  which  was  then  more  easily  broken  by  shear.  Though  the 
viscosity  of  the  solvent( JP-10)  increased  as  the  temperature  decreased, 
the  structure  was  the  dominating  factor  which  determinated  the  viscositv 
at  low  shear  rates.  When  the  structure  became  weaker  at  low  shear 
rates,  the  corresponding  viscosity  became  smaller,  too. 

Structures  are  not  important  at  1  igh  shear  rates  where  nost  of  the 

structure  is  broken  down  under  such  ci  iditions.  Therefore,  the 

viscosity  change  with  temperature  foil >wed  the  traditional  trend,  i.e., 
the  viscosity  increased  as  the  temper. lture  decreased. 

The  temperature  effect  on  the  low  shear  rate  viscosity  was  a  result 
of  the  yield  stress.  The  yield  stress  was  almost  independent  of 
temperature,  rioo  was  highly  dependent  on  temperature.  The  change  of 
nrj3  with  temperature  is  plotted  in  Figure  93. 


C.4  Transient  Behavior 


Transient  response  was  used  to  evaluate  the  viscoelastic  and 
thixotropic  properties  of  the  slurries.  At  low  shear  rates,  the  stress 
growth  after  sudden  imposition  of  a  step  change  in  shear,  showed  an 
overshoot,  and  approached  a  steady  state  value  (Figure  72),  even  at  a 
very  low  shear  rate.  The  "overshoot"  has  been  reported  by  many 
researchers  in  elastic  liquids  like  polymer  solutions.  Slurries  possess 
viscoelastic  properties  at  low  shears,  but  the  elastic  property  differs 
from  that  of  a  polymer.  The  elasticity  in  slurries  comes  from  particles 
and  their  interactions  or  structure.  Since  the  elasticity  in  the  slurry 
was  due  to  the  structure,  an  effort  was  made  to  find  an  elastic  modulus, 
G,  that  could  be  related  to  the  strength  of  structure  in  the  slurries. 
The  transient  curves  were  replotted  as  shear  stress  against  shear  strain 
for  various  shear  rates.  The  shear  strain  was  obtained  by  the  product 
of  shear  rate  and  time,  i.e.,  Y  =  Y  t.  If  the  material  can  be 
described  by  a  =  GY  at  low  shear  strains,  it  would  mean  that  the 
material  acts  as  a  solid  body  with  a  onstant  elasticity  G. 
Unfortunately,  the  initial  slope  for  the  x  vs  y  plot  was  not  constant  as 
shown  in  Figure  94. 

Amari  and  Watanbe(90)  suggested  that  t maxA  was  a  good 
approximation  of  elasticity  G  of  the  structure.  They  obtained  a 
constant  G  of  5.5X103  dyne,  cm2  for  linseed  oil-carbon  black  suspens  on. 
traax/y  for  carbon  black  at  different  shear  rates  is  plotted  in  Figui e 
95.  The  data  are  not  too  far  from  a  constant  and  can  be  considered  as 
an  elasticity  constant,  G,  related  to  the  strength  of  the  structure. 

Stress  relaxation  after  cessation  of  steady  flow  was  also  examined 
to  see  if  it  could  be  described  by  a  single  relaxation  time.  None  could 
be  described  by  a  single  constant.  The  elasticity  might  change  or 
disappear  under  steady  shear.  Other  types  of  study  such  as  creep  flow 
or  dynamic  study  would  be  needed  to  measure  how  the  G  changes. 
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Figure  94.  Shear  Stress  Versus  Shear  Strain  for  Sen  Tech899-U88. 


III.  LUBRICANT  DATA 


The  viscosity  of  lubricant  5P4E  measured  with  Haake  rotor  MVI  was  19.4 
poise  and  with  the  rotor  MVIII  was  19.0  poise. 

Further  tests  on  lubricant  5P4E  were  done  on  the  Weissenberg  R-16 
Rheogoniometer .  The  value  of  the  viscosity,  determined  in  steady  state  flow,  is 
constant  with  shear  rate  varying  from  5.4  to  340  sec-'  and  is  18.6  poise,  a 
value  which  can  be  considered  in  agreement  with  that  obtained  from  the  Haake, 
i.e.,  19.2  poise.  Table  XXVI  summarizes  all  the  lubricant  data  obtained. 

Figure  96  at  1J°C  and  a  high  shear  rate  (348  sec-^)  shows  overshoot, 
which  is  an  indication  of  viscoelasticity.  The  equilibrium  viscosity  for  this 
run  was  233  poise.  Figure  97  is  an  example  of  a  lower  temperature  (2.4°C)  run 
at  a  low  shear  rate  (11  sec-^).  The  equilibrium  viscosity  increased  to  5450 
poise.  No  overshoot  is  observed  at  this  low  shear  rate.  Figure  98  at  -0.7°C 
and  34.8  sec"^  shows  overshoot  and  an  equilibrium  viscosity  of  9790  poise. 
Figures  99  and  100  it  4.9°C  and  2.76  sec-1  ire  examples  of  both  the  growth 
and  relaxation  experiments.  The  equilibrium  viscosity  is  229,700  poise.  The? 
sheer  rate  was  too  low  to  show  overshoot.  Finally,  Figure  101  at  -9°C  and 
onl  ■  2.76  sec-:  is  i  dramatic  example  of  viscoelastic  overshoot. 

The  2.76  sec”!  shear  rate  of  Figure  91  is  equivalent  to  133  sec/rev. 

This  translates  inti  a  movement  of  4.72  mi; /sec  at  the  transducer.  The 
equilibrium  transdu  er  displace  cent  was  calculated  to  be  3.38  mm  which,  if  there 
were  solid  body  rotation,  would  be  accomplished  in  0.72  sec  as  shown  in  Figure 
99.  The  instrument  response  of  0.007  s  ca  mot  even  he  seen  on  this  plot.  The 
material  response  is  noted  as  being  3.15  sec.,  which  is  over  4  times  the  time  of 
a  Newtonian  material.  Clearly  we  are  dealing  with  a  viscoelastic  material.  The 
initial  shear  rate  is  not  really  2.76  sec-1  since  the  top  cone  is  in  motion. 

The  shear  rate  as  determined  from  the  initial  slope  is  only  1.0  sec-1  an  I 
increases  to  2.76  s  ;c-^  by  3.15  sec.  As  a  further  example,  the  shear  rate  had 
increased  to  2.3  sec~^  from  the  measured  slope  at  2  sec.  The  difficulty  in 
measuring  accurate  stress  growth  curves  is  clearly  illustrated  by  this  example. 
Actually,  for  a  material  of  this  viscosity,  we  could  have  used  the  piezoelectric 
crystal  measuring  system  and  obtained  an  accurate  transient.  We  are  currently 
working  on  the  best  means  of  analysis  of  this  type  of  data  so  as  to  extract  the 
vis>  oelastic  characteristics. 

One  approximate  correction,  which  should  be  valid  at  t  =  0,  is  to  correct 
the  observed  stress  by  the  ratio  of  the  d'-sired  shear  rate  to  the  actual  rate. 
Such  a  curve  is  plotted  in  Figure  99.  The  initial  slope  of  this  curve  will  be  a 
measure  of  Gj.  Note  that  the  corrected  cu  ves  show  some  overshoot,  but  his 
may  not  be  real.  In  any  event  we  do  not  plan  to  try  to  analyze  the  entice 
curve;  we  only  want  the  initial  value  as  a  measure  ol  Gj. 

We  do  not  want  to  put  much  importance  on  the  following  values  at  this  time, 

but  from  the  corrected  curve  in  Figure  99,  the  value  of  Gj  is  1.9  x  lU^gm 
i-m/sec^.  From  the  curve  as  measured  and  u  ing  the  G^  value  and 

H0  =  310,000  P  gives  G2  as  -9.1  x  10^  gm  c  i/sec^.  The  value  of  the  time 

derivative  of  the  shear  rate  wa  ;  1.56  sec"  . 
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TABLE  XXVI.  Lubricant  5P4E 


Temp 

Shear  rate 

App.  Viscosity 

File 

(°C) 

(sec-l) 

(poise) 

name 

23°C 

5.5 

18.5 

LUB5D4 

11.0 

18.5 

LUB10D8 

110.  3 

17.8 

LUB108 

174.0 

17.9 

LUB171 

348.0 

17.8 

LUB340 

550.  7 

17.6 

LUB540 

11.4 

11.0 

218.1 

LB10D8SA 

10.6 

110.3 

254.9 

LB108SA 

10.  3 

174.0 

268.4 

LB171SA 

10.2 

348.0 

233.2 

LB340SA 

10.4 

1!  .0 

316.3 

LB10D8SB 

-0.2 

1.103 

24.6  x  103 

LB1D8II1A 

-0.8 

2.757 

28.8  x  103 

LB2D7IIIA 

-0.7 

5.507 

26.7  x  103 

LB5D4IIIA 

-0.0 

11.03 

18.1  x  103 

LB10D8FA 

+0.7 

11.03 

13.1  x  103 

LB10D8FB 

-0.8 

21.93 

20.3  x  103 

LB21D4FA 

-0.4 

0.551 

12.5  x  103 

LBD54FA 

-0.9 

2.757 

14.6  x  103 

LB2D7FB 

-1.2 

11.03 

16.8  x  103 

LB10D8FC 

-1.3 

21.93 

14.3  x  103 

LB21D4FB 

-1.1 

27.57 

12.1  x  103 

LN27FA 

-0.  7 

34.8 

9.79  x  103 

LB34FA 

Temp 

(°C) 


App.  Viscosity 
(poise) 


File 

name 


Shear  rate 
(sec-l) 


-0.2 

0.551 

8.23  x  10 

LUBD54I I 

+0.6 

1.103 

9.47  x  103 

LUB1D7II 

1.4 

2.757 

9.26  x  103 

LUB27II 

2.1 

2.757 

4.63  x  103 

LUB27IIB 

3.0 

5.507 

2.99  x  103 

LUB5D4I I 

2.4 

11.03 

5.45  x  103 

LUB10D8II 

-5.0 

0.551 

288.3  x  103 

LUBD57 

-5.2 

1.103 

290.5  x  103 

LUB1D7C 

-5.0 

1.739 

259.9  x  103 

LUB171 

3 

-4.9 

2.757 

229.7  x  10 

LUB27A 

sponse  curve 
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Figure  90.  5P4E  lubricant 
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APPENDIX  B 


appendix  c 


NOMENCLATURE 

a  constant  In  conversion  term  in  kinetic  theory 

A,B,C  reaction  components 

C  degree  Centigrade,  concentration 

CSR  constant  shear  rate 

d  displacement 

F»F*»Fo  fraction  converted 

A  D 

G,G^,G2  elastic  constants 

I  t 

k^,k2»k^,k2» • • •  specific  rate  constants  in  kinetic  theory 
K,Ka,Kb  constant  in  kinetic  theory 

constants  in  Eq.  (2) 

m,ia.,  .n^.n.n^,^  order  constants  in  the  kinetic  theory 

M  critical  molecular  weight 

c 

N  number 

p,p^,Pg  constant  in  the  kinetic  theory 

Pl*^2’‘”  sensitivity  to  stress  constant  in  kinetic  theory 

S  parameter  in  two-step  kinetic  theory 

t  time 

t  „  time  when  G1  starts  to  decay 

cut 

TAU  stress 

T 


V 


temperature 

voltage 


Greek 


6 

Y 

n 

n. 


constant  for  non-constant  shear  rate  experiment 

shear  rate 

viscosity 

lower  Newtonion  viscosity  limit 
upper  Newtonion  viscosity  limit 
thixotropic  viscosity 
stress 


Subscripts 

calc 

eq 

exp 

n 

0 

t 


calculated 

equilibrium 

experimental 

normalized 

initial 

thixotropic 


# 


ft 
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SUMMARY  OF  THE  FIELD  AMD  OUR  RESEARCH  ENDEAVORS 


The  various  approaches  towards  the  elucidation  of  the  rheological 
characteristics  of  materials  generally  fall  into  one  of  five 
classifications:  empirical,  phenomenological  (mostly  rate  processes), 
linear  viscoelastic,  nonlinear  viscoelastic,  and  microrheological  analyses. 
The  empirical  methods  correlate  data  by  curve-fitting  techniques.  The 
rate  theories  have  as  their  basis  the  assumption  that  the  nonlinear  char¬ 
acteristics  can  be  associated  with  some  structural  change  of  the  material 
whether  it  involves  particle  associations  and  dissociations,  link  formations 
and  ruptures,  or  molecular  entanglements  and  dis entanglements.  The  linear 
viscoelastic  models  are  based  upon  linear  combinations  of  Hooke’s  law  of 
elasticity  and  Newton's  law  of  viscosity.  The  nonlinear  viscoelastic 
models  are  based  on  continuum  mechanics  ;  nd  nonlinear  combination  mechanical 
models.  Finally,  microrheological  analysis  starts  with  the  basic  molecular 
or  microscopic  variables  such  as  particle  sizes,  molecular  interactions, 
and  chain  lengths.  A  simplified  mechanism  is  proposed,  and  then  it  is 
mathematically  represented  and  solved. 

Each  approach  has  made  contributions  to  the  field,  but  with  the  present 
state  of  the  art,  we  cannot  completely  describe  the  non-Newtonian  behavior 
even  in  simple  geometries.  Steady-state  flow  behavior  has  been  extensively 
investigated,  and  the  representations  available  are  adequate.  However, 
little  is  understood  about  the  unsteady- state  flow,  such  as  shear  stress 
and  normal  stress  growth  after  the  onset  of  a  sudden  shear  rate.  The 
theories  in  the  literature  are  not  able  to  always  predict  for  example  stress 
growth  at  a  given  constant  shear  rate  from  data  obtained  on  other  experi¬ 
ments. 

The  kinetic  rate  theory  we  are  developing  can  be  used  to  represent 
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transient  (unsteady-state),  equilibrium,  and  normal  stress  difference 
information.  Qhe  theory  accounts  for  elasticity  and  thixotropy  by  combining 
a  kinetic  rate  approach  with  a  modified  form  of  Oldroyd's  model.  The  basic 
approach  and  preliminary  confirmation  of  the  kinetic  theory  have  been 
described  in  the  papers  by  Brodkey,  Denny  and  Kim  (2,  3)-  Kim  and  Brodkey 
(3)  provided  the  test  of  the  steady-state  theory  by  correlating  polymer 
solution  data  over  a  wide  concentration  range.  Lewis  and  Brodkey  (5) 
offered  further  support  of  the  theory  by  describing  the  phenomena  of 
shear  stress  growth  over  a  range  of  shear  rates.  However,  thus  far*  the 
tests  of  the  theory  were  still  empirical  since  the  necessary  constants  and 
parameters  had  been  obtained  from  the  best  fit  to  the  data.  Lee  and 
Brodkey  (7)  attempted  to  provide  an  independent  test  of  the  kinetic 
approach  by  using  data  obtained  under  constant  shear  stress  conditions 
to  evaluate  the  necessary  constants  and  parameters  in  an  idealized  manner 
as  suggested  by  the  theory,  and  not  by  simple  curve-fitting  methods.  The 
results  were  then  used  for  predicting  time-dependent  constant  shear  rate 
measurements  as  well  as  the  constant  stress  data.  Jachimiak  et  al.  (8) 
and  Jachimiak  (9)  extended  the  approach  to  allow  the  prediction  of  the 
first  normal  stress  difference;  but,  because  of  the  difficulty  of  obtain¬ 
ing  the  elastic  parameter  under  constant  stress  conditions,  the  prediction 
was  unsatisfactory.  Song  (10)  did  a  complete  analysis  of  the  theory  and 
showed  that  the  previous  formulations  could  not  be  used  to  predict  all  the 
various  possible  types  of  experiments.  He  was  able  to  formulate  one  that 
could  be  used  which  will  be  presented  in  the  theory  section.  Song  (11) 
then  obtained  a  comprehensive  set  of  data  on  stress  growth,  stress 
relaxation  and  transient  shear  rate  experiments.  He  attempted  to  fit  all 
of  these  to  the  theory  with  some,  but  not  complete,  success.  Pan  1  ala i  (12) 
obtained  data  for  very  large  amplitude  oscillation  and  extended  the  theory 


for  this.  Finally,  Landers  (13)  obtained  new  and  improved  oscillatory  data 

on  the  same  system.  The  most  recent  is  by  Park  (75)  and  is  summarized  in  App. 

The  present  summary  is  slanted  so  far  toward  our  own  work,  and  the 
background  literature  has  been  covered  in  several  of  our  previous  publi¬ 
cations.  However,  it  might  be  well  to  put  our  approach  into  proper 
perspective  with  other  means  of  treating  rheological  data,  since  the  area 
of  research  is  large  and  all  readers  may  not  be  equally  familiar  with  our 
efforts.  First,  there  are  the  empirical  methods  such  as  the  poy-r  or 
Ostwald-deWaele  (l4,  15),  Ellis  (1 6,  17),  Sisko  (l8),  Reiner-Philippoff 
(19),  Powell-Eyring  (20),  and  Rodriguez  (21)  models.  Detailed  explanation 
of  these  c^n  be  found  in  most  textbooks  (1,  22,  23).  The  semiempirical 
or  phenomenological  models  form  a  second  class  of  which  the  Eyring  work 
(20,  24-26)  is  representative.  It  is  in  this  class  that  the  present 
effort  resides  i.s  well  as  a  number  of  other  kinetic  approaches  (27-30). 

These  latter  (27-30)  have  been  less  general  kinetic  approaches  that  have 
assumed  some  specific  mechanism  and  then  suggested  a  rate  equation  to 
describe  the  change;  however,  since  these  are  specific  to  a  particular 
suggested  mechanism  and  not  general  to  aid  in  establishing  the  mechanisr , 
they  will  not  be  discussed  further.  The  basic  linear  viscoelastic  model 
is  due  to  Clerk  Maxwell  (31) •  The  equation  is  a  linear  combination  of 
Newton's  law  of  viscosity  of  an  ideal  liquid  and  Hook's  law  for  an  ideal 
solid.  Even  though  it  is  a  linear  equation,  it  has  been  a  starting  point 
of  almost  all  the  nonlinear  rheological  models  proposed  in  the  last  two 
decades.  Most  of  the  models  for  nonlinear  viscoelasticity  are  a  generali¬ 
zation  of  Maxwell's  equation  using  tensor  calculus.  Typical  equatior s  of 
this  category  are  discussed  in  references  32  through  55 •  Most  are 
formulated  in  a  convected  (co-deforr  m: )  coordinate  system  as  was  first 


suggested  by  Oldroyd  (3*0*  Some  are  expressed  in  a  co- rotational  coordinate 
system,  which  has  been  reviewed  by  Bird  et  al.  (53)*  In  spite  of  all  the 
work,  the  ability  to  fit  transient  data  is  far  from  satisfactory  as  can  be 
seen  in  Chen  and  Bogue  and  Christiansen  and  Leppard  and  others  (51,54,55, 

65  -  67).  Figures  1  and  2  show  some  comparisons  from  the  latter  work. 
Finally,  there  are  the  microrheological  approaches  which  assume  a  mechanism 
and  are  directed  toward  elucidating  the  basic  reason  why  specific  systems 
are  non- Newtonian.  In  this  class  are  such  analyses  as  Rouse  (56),  Bueche 
(57),  Michaels  et  al.  (5S-6O),  and  Mason  (6l,  62). 

In  the  empirical  models,  nothing  is  learned  about  the  basic  mechanism 
because  little  specific  meaning  can  be  ;  iven  to  the  constants.  In  the 
phenomenological  approaches  as  little  is  assumed  about  the  specific 
mechanism  a.:  possible,  and  cae  hopes  to  -ain  some  insight  into  this  by 
evaluation  of  the  constants  and  their  variation  with  the  parameters  of 
the  system  (such  as  raolecultr  weight  iis tribution } .  The  viscoelastic 
theories  have  been  the  recent  mainstay  f  rheological  representation  and 
as  many  believe  will  eventually  be  ale  puate  far  the  task.  The  micro- 
rheological  analyses  depend  >n  a  knowl  tge  ,f  the  m*  cnanism  and  are  only 
valid  fir  systems  that  fulfill  the  di;  •  issompt.i  -n.  in:  ij.prox Lmations  ,f 
the  theories.  In  general,  a  mechanism  is  suggested,  .no  then  >ne  attempts 
to  confirm  the  hypothesis  by  comparison  to  experiment  1;  lat j..  •  ,pe fully, 

in  time  the  methods  will  complement  eacr.  otrw  r ,  :  .om  &  ri  lv:  :1  if.:  .ides 
providing  the  necessary  insight  to  allow  a  more  r-  as  nit  i-  me  ;:,  if. is::.  • 
proposed  than  has  heretofore  been  possible,  and  micr  •••..  •  .1  . 

allowing  estimation  and  •  liigent  variation,  at  ♦:.•  m/- s*  .  .♦ 
f  the  characteristic  pa  ,ers  of  mate  ials. 
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Typicol  shear  stress  decoy  with  the  predictions  of  the 
Sprifps,  Corrcou,  and  Bogus-Chen  models. 


Both  figures  reproduced,  from  reference  55 
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The  possible  application  of  our  kinetic-elastic  model  (as  well  as 
alternate  models)  has  been  mainly  to  polymer  rheology  problems  where  time 
dependency  is  recognized  but  is  not  as  extreme  as  it  might  be  in  elasto- 

O 

hydrodynamic  lubrication  (^ls  versus  ^  10  s).  Temperature  dependency 

is  important  in  both,  but  pressure  effects  are  much  more  extreme  in  the 
lubrication  application.  Johnson  and  Tevaarwerk  (69)  have  recently 
reviewed  the  shear  behavior  of  EHD  oil  films,  but  it  should  be  emphasized 
that  the  models  considered  by  them  can  at  best  only  describe  a  part  of 
the  known  response  characteristics  of  non-Newtonian  materials.  Never¬ 
theless,  the  reasonable  success  they  demonstrated  in  matching  traction 
curves  leads  one  to  hope  that  a  better  model,  more  accurately  descriptive 
of  the  process,  would  allow  nearly  a  complete  description  of  the  dynamics 
of  the  problem. 

The  model  used  by  Johnson  and  Tevaarwerk  was  the  simplest  form  of  the 
sinh  law  for  the  non-Newtonian  representation  (see  references  20  and  ’4-26) 
and  the  Maxwell  model  (31)  for  the  viscoelastic  contribution.  In  both 
cases,  the  representation  of  known  response  characteristics  is  limited  and 
does  not  compare  with  the  kinetic  -elastic  model  to  be  described.  There 
are  of  course  a  number  of  earlier  models  that  have  been  used  to  describe 
traction  curves.  Crook  (70)  suggested  a  Newtonian  fluid  and  temperature 
dependence  to  explain  the  nonlinear  region.  Johnson  and  Cameron  (71) 
showed  that  non-Newtonian  characteristics  were  indeed  important.  Barlow 
et  al.  (72)  showed  that  lubricants  could  be  viscoelastic.  Byson  (73) 
suggested  that  viscoelasticity  could  be  important  in  the  linear  region 
under  large  deformations  as  well  as  in  the  nonlinear  region. 

Ihe  knowledge  of  the  characteristics  of  lubricants  is  clearly  limited 
and  the  best  means  to  describe  the  fluid  during  a  concentrated  contact  is 
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still  to  be  determined  thus  there  exists  a  need  for  looking  into  this  applied 
problem  in  more  detail  using  the  Theologically  sound  kinetic-elastic  approach. 

The  application  of  the  theory  to  the  field  of  slurry  fuels  is  still  another 
possibility.  These  materials  will  usually  exhibit  a  yield  stress.  Unfortunately 
data  on  the  time  response  of  such  materials  does  not  appear  to  have  been 
measured  and  clearly  work  is  needed  here. 


THEORY 


Polymer  rheology  is  primarily  concerned  with  three  fundamental 
properties  of  fluids:  the  time  rate  of  change  of  viscosity  or  structure, 
called  here  the  thixotropic  property;  the  steady-state  level  obtained, 
which  is  a  function  of  the  shear  rate;  and  the  lag  of  the  stress  behind 
what  one  would  expect  from  the  thixotropic  behavior,  called  here  the 
elastic  property.  The  change  in  viscosity  can  be  attributed  to  some 
reversible  internal  mechanism  such  as  entanglement,  association,  or  winding 
of  long-chain  polymer  molecules.  The  lag  of  the  stress  is  due  to  the 
property  of  the  material  which  causes  it  to  resist  deformation  and  thereby 
recover  its  original  shape  and  size  when  deforming  forces  are 
removed.  The  nonlinear  time-dependent  flow  behavior  for  such  polymeric 
materials  can  best  be  described  by  a  combination  of  elasticity  and 
thixotropy.  The  thixotropic  and  steady-state  behavior  is  described  by  the 
kinetic  theory  and  the  elastic  property  is  described  by  an  elastic  model. 


A  KINETIC  MODEL  FOR  THIXOTROPIC  CHANGE 


In  the  previous  work  (1-13),  we  developed  a  kinetic  model  for  thixo¬ 
tropic  materials  which  is  descriptive  of  the  viscosity  change  that  occurs 


in  polymer  solutions.  The  model  is  different  from  previous  rate  models 


in  that  molecular  phenomena  such  as  links,  particle  spheres,  dipole 
moments,  or  hydrogen  bonding  are  not  used  directly  as  a  basis  of  the  model. 
Instead,  let  us  define  a  new  variable,  F^,  as  the  fraction  of  the  thixotropic 
fluid  structure  unchanged.  This  new  variable  F^  is  a  lumped  parameter 
which  contains  the  overall  effects  of  molecular  phenomena  such  as  disentangle 


ment  and  entanglement  of  chains,  link  rupture  and  association  and  others. 


At  zero  time  and  no  stress,  the  fluid  thixotropic  structure  would  be 
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unchanged  and  =  1.  As  time  progresses  (unsteady  state),  the  thixo¬ 
tropic  fluid  structure  undergoes  change  and  F^  decreases  and  continues  to 
decrease  as  the  structure  changes.  Finally,  at  steady  state,  an  equili¬ 
brium  is  reached  between  the  unchanged  and  changed  thixotropic  fluid 
structure  and  F^  =  Fg^,  the  equilibrium  value  of  the  variable.  The 
situation  being  considered  reminds  one  of  a  simple  kinetic  reaction  of  Ajffi. 
If  we  apply  the  homogeneous  kinetic  reaction  concept,  the  simplest  equation 
we  could  get  is 


d(CF. ) 

-af  -  4  <cpt>"  -4  rc(i-Ft)]n  (u 


where  k|  and  k^  are  forward  and  reverse  reaction  constants,  m  and  n 
are  forward  and  reverse  orders  of  reaction,  and  C  is  the  concentration  of 
polymer  solution.  For  steady  state  (equilibrium)  conditions,  Eq.  (1) 
reduces  to 


4  <CV“  -4  " 0 


(1-F  )  *• 

y*  -  4 

where  the  subscript  eq  means  that  F^  is  at  the  equilibrium  state. 

What  we  are  eventually  interested  in  is  not  the  thixotropic  structure 
unchanged,  but  rather  the  real  fluid  property  viscosity;  therefore,  we 
need  some  functional  form  of  F^  expressed  in  terms  of  viscosity.  The 
limiting  viscosities  (zero  shear  rate  viscosity,  7,  ,  and  upper  Newtonian 
viscosity,  Tj  )  can  be  utilized  to  define  the  thixotropic  structure  unchanged, 

CO 

F^,  since  the  apparent  viscosity  of  a  polymer  solution  falls  within  these 
two  limiting  values.  The  structural  variable  F^  is  defined  by  an  inverse- 
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lever  arm  principle  (1,3): 


f(Tit)  -  t(\) 


where  T|  can  be  called  the  thixotropic  structural  viscosity.  The 
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functional  form  taken  for  f(T|)  is  &f]  ,  which  when  combined  with  Eq.  (4) 


gives : 


\  ^ 
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where  a  is  constant  to  be  determined.  The  constant  a  can  be  selected  to 
be  (1/3-4)  from  experimental  evidence  (63)  that  the  zero  shear  rate  viscosity 
T)o  is  proportional  to  the  3-4  power  of  the  molecular  weight.  But,  for 
specific  cases  or  other  materials  any  value  which  gives  a  good  fit  to  the 
experimental  data  would  be  correct.  For  slurries  a  value  of  1  can  be  used. 

To  complete  the  kinetic  model,  we  need  one  more  basic  functional 
assumption  about  the  reaction  constants,  k|  and  kg,  in  Eq.  (1).  In 
general,  these  reaction  constants  are  considered  to  be  functions  of  shear 
stress  and  temperature.  For  isothermal  condition,  the  following  functional 
forms  are  suggested  from  empiricism  and  logic  (see  references  1  to  II ) 


ki  -  h T 


where  k^  and  kg,  and  the  p's  can  be  said  to  be  rates  and  susceptibilities 
to  stress  level,  respectively.  The  term  k^  is  a  constant  in  this  case 
since  it  is  considered  to  be  a  function  of  temperature  only.  It  accounts 
for  the  physical  fact  that  a  material  recovers  its  original  structure 


from  molecular  or  Brownian  motion  without  other  external  help.  It  might 
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be  considered  small  and  thus  negligible  compared  with  other  terms,  except 
when  shear  stress  is  zero  as  a  result  of  the  material  simply  standing  for 
long  periods  of  time  in  order  to  allow  structural  recovery.  Equation  (1) 
is  a  scalar  equation  and  the  stress  in  Eqs.  (6)  and  (7)  must  be  scalar. 

If  the  structure  is  only  affected  by  the  shear  stress  then  t  in  Eqs.  (6) 
and  (7)  must  reflect  this.  If  the  structure  depends  on  the  total  stress 
in  the  system  then  the  double  dot  product  of  the  stress  (t:t)  might  be 
more  appropriate. 

We  define  a  thixotropic  structual  stress,  t^,  as 

r t  =  'Tit  Y  (8) 

This  stress,  T  ,  follows  the  dashed  line  in  Fig.  3-  However,  the  real 
stress,  t,  follows  the  solid  line  in  the  same  figure.  There  is  a  gap 
between  the  real  and  newly  defined  thixotropic  structural  stress.  To 
account  for  this  difference  one  needs  to  consider  the  elastic  property 
of  the  material. 

B  ELASTIC  MODEL  FOR  VISCOELASTICITY 


The  kinetic  theory  of  thixotropy  just  presented  assumed  that  there 
was  no  elasticity  in  the  material;  however,  most  polymeric  solutions  are 
viscoelastic  and  the  thixotropic  presentation  alone  is  not  adequate. 

In  order  to  incorporate  the  elastic  property  into  the  theory,  a 
three  constant  model  which  is  similar*  to  Maxwell's  equation: 


T  + 
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G  dt 
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(9) 


and  originally  proposed  by  Jeffreys  (64)  and  used  by  Oldroyd  (34,35)  is 
used.  The  equation  is 


T  +  XJL(dT/dt)  =  -t1q  fy  +  X2(dy/dt)] 

(10) 

where  X^  and  are  relaxation  times  and  T)0  is  zero  shear  rate  viscosity. 
In  our  earlier  work  (5j  7)  an  integrated  form  of  Eq.  (10)  was  used.  The 
modified  version  of  Eqs.  (9)  and  (10)  currently  used  is 


T  +  (H^O,*  -  -Ht  tv  +  (yo2>?5 


where  G^  and  are  modulus  type  parameters  to  be  determined.  Connected 
derivatives  are  used  so  that  normal  stresses  can  be  accounted  for.  This 
equation  will  be  called  the  elastic  model  since  it  is  developed  to  account 
for  the  lagging  or  elastic  part  of  viscoelasticity.  Equation  (l)  was 

called  the  kinetic  model  since  it  is  formulated  with  the  kinetic  concept. 
Thus,  the  rate  theory  is  called  the  kinetic -elastic  model. 

C  LIMITING  CONDITIONS  FOR  THE  KINETIC-ELASTIC  MODEL 

C.l  STEADY-STATE  CONDITIONS 

Introducing  Eqs.  (6)  and  (7)  into  Eq.  (3)  gives 


(1  -  p  )n 

— sa'  =  k  c  m_n  tp 
f1 

eq 

(12) 
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K  "  k2  +  (k3/TP27 

(13) 

P  =  P1  '  P2 

(14) 

The  subscript  eq  means  the  equilibrium  value. 

It  has  been  found  tint  K  is  essentially  constant  (3)*  This  would  imply 
that  under  equilibrium  conditions  where  t  is  finite  k^/TP2  is  small  compared 


to  kg  and/or  that  Pg  is  very  close  to  zero,  so  that  K  =  k^/(kg  +  k^). 
Apparently  the  latter  is  true  since  Pg  has  been  found  to  be  very  close  to 
zero  (7,  9»  10,  11).  In  any  event  K  and  p  can  be  taken  as  constants  over 
the  entire  range  of  shear  rates  or  stresses.  It  has  been  found  that  m  -  1 
and  n  =  2  is  satisfactory;  thus,  if  so,  Eq.  (12)  becomes 

Feq  =  1  +  (  K/2C)tP  -/( 1  +  (K/2C)tP)2-1 

(15) 

If  is  negligible  compared  withT|0  the  non- Newtonian  (steady- state) 
viscosity  expression  from  Eqs.  (5)  and  (15)  is 

I)  =  T]o  [1  4-  (K/2C)tP  -  /(1-!  (k/2C)tP)2-1] 

(16) 


C.2  SHEAR  STRESS  GROWTH  AT  CONSTANT  SHEAR  RATE 


At  constant  shear  rate  and  considering  only  the  shear  stress,  the 
elastic  Eq.  (11)  can  be  simplified  to 

dTn/dt  =  tn.  (Gj/7|t)(Ttn  -  Tn)  (|7) 

"here  Tn  m  Ttn  -  Tt/Teq. 

The  other  equations  are: 


d(CF.) 
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with  k£  and  kg  given  by  Eqs.  (6)  and  (7) 


1/a 


(14) 
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The  limiting  conditions  are  at  time  =  0, 


mi 
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and  at  time  =  «  or  equilibrium 


Ft  "  1,0»  Tn  "  °»  \  ~  \ 


F.  =  constant  <  1,  t  =  1.0,  7).  =  T1 
t  =  n  't  1 


eq 


C.3  SHEAR  STRESS  RELAXATION  AFTER  CESSATION  OF  SHEAR 


The  elastic  Eq.  (11)  takes  the  simple  form  (shear  rate  equal  to  zero) 

d  rn/dt  =  -  (G1/T]t)  Tn  (19) 


Equations  (l)  and  (l8)  in  the  previous  section  apply  here.  The  limiting 
conditions  are  at  t  =  0, 

F.  -  F  <1.0,  T|  =  Tl 
t  eq  ’  't  ‘eq 

and  at  t  ~  , 


Ft  =  1.0,  T)t  =  'lo 

C.4  TRANSIENT  SHEAR  RATE  AT  CONSTANT  SHEAR  STRESS 

For  this  case,  a  little  manipulation  of  the  elastic  Eq.  (11)  at 
constant  stress  generates  Eq.  (20): 


^2  .  y'  _  QJ_ 
dt  Yn  TL 
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where  ^  is  (-  y)/(-  ;gq). 


and 


The  other  pertinent  equations  are, 
d(CF.  ) 
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The  boundary  conditions  are  at  t  =  0, 


F.  =  1.0  and  y  is  guessed  by  data 

O  'll  ‘ 


and  t  = 

F.  =  F  and  v  =  1 
t  eq  *n 

The  special  feature  of  this  case  is  that  Eq.  (1)  is  not  coupled  with 
Eq.  (20).  Thus,  Eq.  (1)  can  be  numerically  integrated,  independently, 
'"his  information  is  then  utilized  for  the  integration  of  Eq.  (20). 

C. 5  NORMAL  STRESS 

In  a  straight  forward  manner,  as  shown  by  Oldroyd  (see  1),  one  can 
obtain  from  the  model  in  the  form  of  convective  derivatives,  the  first 
normal  stress  difference  as 


where  the  latter  form  is  obtained  using  Eq.  (8).  In  references  (8)  and 
(9),  was  not  satisfactorily  obtained,  thus  an  adequate  test  of  the 
prediction  was  not  made. 

C.6  OSCILLATORY  CONDITIONS 

The  kinetic  equation  (1)  along  with  the  equation  to  account  for  the 
elastic  part  of  the  deformation  (11)  (using  a  sinusoidal  input  expression 
for  shear  rate)  can  be  solved  simultaneously.  For  the  oscillatory  shear 
rate  and  its  derivative  one  has 


Y  =  Vp  +  v0  sin  (cut) 

(22) 

Y  =  y0co  cos  (out) 

(23) 

where  Vp  is  a  fixed  shear  rate  level.  Neglecting  Tj^  compared  to  T|0  and  7)^ 
in  Equation  (5),  we  have 


Therefore, 


Ft  -  <K 

(24) 

dVat  -  \  \ 

(25) 

Substituting  the  above  equation  in  Eq.  (1),  using  Eqs.  (6)  and  (7)  with 
=  0  and  rearranging  gives 


Rearranging  Eq.  (11)  gives 


(26) 


(27) 


Equations  (26)  and  (27)  can  be  simultaneously  solved  numerically  fusing 
Eqs.  (22)  and  (23)]  to  get  time  vs.  sheer  stress  and  7)  and  the  results 
compared  to  values  calculated  from  experimental  data. 


D  SOME  SPECIFIC  RESULTS 


Figures  4  and  5  show  typical  predictions  obtained  from  the  theory  for 
stress  growth  and  relaxation.  Typical  constants  used  for  these  are: 

From  equilibrium  data,  assumed  m=l,  n=2,  a=l  and  7]^  <  7) 

h 

evaluated  experimentally  7]  =  3«1  x  10  poise 

for  35$  PMMA  (mw  =  1.15  x  ^J5)  in  DEP 
evaluated  p  =  1.621  and  log  K  =  -8.936 


From  constant  shear  rate  transient  results 

P2  =  -O.Ol  k2  =  20.0 

=  f(y)  and.  is  about  1.0  x  10^ 

From  constant  shear  stress  transient  results 

G2  =  F(t)  and  is  about  0.8  x  10^ 

As  can  be  seen  from  the  figures  the  trend  of  the  predictions  are  correct. 

These  specific  results  are  not  optimized  to  fit  our  full  range  of  lata 

and  thus  any  comparison  would  not  be  meaningful. 

Some  limited  and  specific  comparisons  can  be  cited  from  our  previous 

work  (3,  5,  7,  8).  Figure  6  is  an  example  of  the  basic  shear  diagram 

compared  to  the  theory  for  a  range  of  concentrations  (3).  Figure  7  gives 

similar  results  for  the  polymethylmethacrylate  solution  (7)*  Figure  8 

(0.44  sec-1)  and  Figure  9  (8*68  sec"1)  show  the  empirical  fit  obtained  by 

2 

Lewis  (5),  while  Figure  10  (496  dynes/om  )  shows  the  shear  rate  build-up 
and  the  fit  obtained  by  Lee  (7)*  These,  as  noted,  are  the  earlier  empirical 
fits  based  on  an  integrated  form  of  Eq.  (10)  and  are  not  based  on  the 
latest  equations  given  in  this  proposal.  Again,  evaluation  of  the  present 
equations  with  our  data  is  the  subject  matter  of  this  proposal. 

Pandalai  (12)  was  able  to  reproduce  large  amplitude  oscillatory  data 
superimposed  on  a  steady  shear  rate  by  using  the  theory  as  presented  in 
Eqs.  (26)  and  (27)  and  a  set  of  constants  selected  to  fit  the  oscillatory 
data.  The  set  of  constants  was  not  necessarily  unique  since  no  overall 
optimization  survey  was  made.  Pandalai  was  not  able  to  reproduce  stress 
growth  and  relaxation  experiments  with  exactly  the  same  constants  but  was 
able  to  do  reasonably  well  by  adjusting  only  G^  and  k2» 


BASIC  SHEAR  DIAGRAM 
FOR  POLYACRYLIC  ACID 


5 


Fig.  ? 


BASIC  SHEAR  DIAGRAM 
FOR  POLYMETHYL¬ 
METHACRYLATE  IN 
DIETHYLPHTALATE 
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COMPARISON  OF 
DATA  AND  RESULTS 
FOR  O.L36  sec-1 


Fig.  9 

COMPARISON  OF 
DATA  AND  RESULTS 
FOR  8.68  sec 


Fig.  10 

COMPARISON  OF 
DATA  AND  RESULTS 
FOR  CONSTANT 
STRESS  EXPERIMENT 


The  kinetic-elastic  theory  is  quite  general  and  can  be  used  as  a 
basis  for  describing  lubricating  materials.  As  long  as  the  mechanism  of 
the  flow,  i.e.,  the  reasons  for  the  nonlinearity  of  the  fluid,  does  not 
change  the  model  should  be  valid.  Ideally,  one  would  like  to  make  measure¬ 
ments  of  the  rheological  characteristics  under  isothermal  and  atmospheric 
pressure  conditions  and  to  be  able  to  use  this  information  under  conditions 
of  temperature  variation  and  extreme  pressure  levels.  The  degree  that 
the  model  can  be  used  to  describe  such  changes  is  unknown  and  must  be 
determined. 

Certainly  the  time  response  is  realistically  described  by  the  kinetic- 
elastic  model.  Our  past  work  hrs  demonstrated  this  for  polymers  and  their 
solutions  and  there  is  no  reason  to  expect  lubrication  materials  to  be 


different. 


It  should  be  noted  that  all  alternate  means  of  analysis  have  been  for 
a  Maxwell  type  model  which  means  the  derivative  of  the  shear  rate  in 
Eq.  (11)  has  been  neglected  or  effectively  G2  was  taken  as  infinite.  Our 
experiments  show  that  the  latter  is  not  true  and  certainly  time  rate  of 
changes  of  shear  rate  exists  in  lubrication  systems,  so  that  one  wonders 
what  is  the  effect  of  only  using  an  equation  like  Eq.  (9)  rather  than  (10) 
or  its  generalization,  Eq.  (11). 

Temperature  changes  are  introduced  naturally  through  the  values  for 
the  lower  and  upper  Newtonian  viscosities  and  through  the  rate  constants 


(k^)  and  the  equilibrium  constant  (K)  (see  reference  1).  The  rate 
constants  and  the  limiting  viscosities  depend  on  temperature  through 
equations  of  the  form 


-  v  'E%1 


T1  .  A’e  -B/RT 
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where  E*  is  an  activation  energy  of  the  non-Newtonian  reaction  and 
k  ,  A’,  B  are  all  constants.  A  similar  equation  is  used  for  T]  and  K 
obeys  a  van  Hoff  type  equation  of  thermodynamics  (i.e.,  from  kj/l^). 

The  elastic  parameters  are  not  expected  to  depend  on  temperature,  but 
this  will  have  to  be  verified  by  experimentation.  We  have  done  some 
work  on  the  effect  of  temperature  on  steady-state  results,  but  have  not 
as  yet  tried  to  extend  this  type  of  analysis  to  predict  results  at  one 
temperature  from,  another. 

The  effect  of  polymer  concentration  and  the  adequacy  of  the  model 
to  describe  steady  state  conditions  was  demonstrated  with  Fig.  6  where 
data  over  a  wide  range  of  concentrations  was  represented.  It  would  be 
interesting  to  establish  if  the  effect  of  pressure  on  the  time  response 
problem  is  simply  a  change  in  concentration  (i.e.,  density  in  gm/cc)  in 
the  kinetic  equation  (1)  or  does  it  also  involve  a  coupling  through  the 
normal  stress  given  in  Eq.  (21).  Of  course,  changes  in  the  limiting 
viscosities  with  pressure  would  have  to  be  established  by  experiments 
using  for  example  capillary  viscometers. 


F.  APPLICATION  TO  SLURRY  FUELS 


We  have  not  been  able  to  locate  as  yet  any  existing  data  that  is  adquate 
for  the  time  response  analysis  of  slurry  materials  with  a  yield.  However,  some 
steady  state  (equilibrium)  data  is  available  and  this  has  been  analysed.  As 
far  as  we  can  tell,  slurry  materials  with  a  yield  can  be  treated  by  the  theory 
by  using  a  minor  modification  to  Eq.  (5)  for  the  definition  of  the  fraction 
converted.  This  is  necessary  because  n0  for  such  materials  is  infinite  (i.e., 
solid-like).  The  modifications  attempted  for  the  steady-state  data  that  proved 
promising  were:  1)  a  simple  shifting  of  the  stress  axis  to  make  the  yield  point 
(to)  the  origin,  2)  using  a  in  Eq.  (5)  as  negative  (so  that  the  n0  term  would  be 
zero;  i.e.  l/nua),  and  3)  using  ’i  defined  as  the  slope  of  the  curve  rather  than 
by  Newton’s  law. 


